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This new edition is a self-contained, concise but rigorous book 
introducing the reader to the basics of the subject. It has been 
brought completely up to date and reorganized to improve the 
quality and flow of the material. 

The introductory chapters provide definitions and useful 
mathematical and physical notes to help readers understand the 
basics. The book then describes the topics relevant to atmospheric 
processes, including the properties of moist air and atmospheric 
stability. It concludes with a brief introduction to the problem of 
weather forecasting and the relevance of thermodynamics. Each 
chapter contains worked examples to complement the theory, as 
well as a set of student exercises. Solutions to these are available 
to instructors on a password protected website at 
www.cambridge.org/9780521696289. 

The author has taught atmospheric thermodynamics at 
undergraduate level for over 20 years and is a highly respected 
researcher in his field. This book provides an ideal text for short 
undergraduate courses taken as part of an atmospheric science, 
meteorology, physics, or natural science program. 
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and their application in climate, climate variability, predictability, 
and nonlinear time series analysis. He is a member of the 
American Geophysical Union and the European Geosciences 
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We have explained that the causes of the elements are four: 
the hot, the cold, the dry, and the moist. In every case, 

heat and cold determine, conjoin, and change things. Thus, 
hot and cold we describe as active, for combining is a sort 
of activity. Things dry and moist, on the other hand, are the 
subjects of that determination. In virtue of their being acted 
upon, they are thus passive. 


Aristotle, Meteorology, Book IV 
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PREFACE 





This book is intended for a semester undergraduate course in 
atmospheric thermodynamics. Writing it has been in my mind 
for a while. The main reason for wanting to write a book like 
this was that, simply, no such text in atmospheric thermodynam- 
ics exists. Do not get me wrong here. Excellent books treating 
the subject do exist and I have been positively influenced and 
guided by them in writing this one. However, in the past, atmo- 
spheric thermodynamics was either treated at graduate level or 
at undergraduate level in a partial way (using part of a general 
book in atmospheric physics) or too fully (thus making it diffi- 
cult to fit it into a semester course). Starting from this point, my 
idea was to write a self-contained, short, but rigorous book that 
provides the basics in atmospheric thermodynamics and prepares 
undergraduates for the next level. Since atmospheric thermody- 
namics is established material, the originality of this book lies in 
its concise style and, I hope, in the effectiveness with which the 
material is presented. The first two chapters provide basic defini- 
tions and some useful mathematical and physical notes that we 
employ throughout the book. The next three chapters deal with 
more or less classical thermodynamical issues such as basic gas 
laws and the first and second laws of thermodynamics. In Chapter 
6 we introduce the thermodynamics of water, and in Chapter 7 
we discuss in detail the properties of moist air and its role in 
atmospheric processes. In Chapter 8 we discuss atmospheric sta- 
bility, and in Chapter 9 we introduce thermodynamic diagrams 
as tools to visualize thermodynamic processes in the atmosphere 
and to forecast storm development. Chapter 10 serves as an epi- 
logue and briefly discusses how thermodynamics blends into the 
weather prediction problem. At the end of each chapter solved 
examples are supplied. These examples were chosen to comple- 
ment the theory and provide some direction for the unsolved 
problems. 


PREFACE 


Finally, I would like to extend my sincere thanks to Ms Gail 
Boviall for typing this book and to Ms Donna Genzmer for draft- 
ing the figures. 


Anastasios A. Tsonis 
Milwaukee 


CHAPTER ONE 





Basic definitions 


e Thermodynamics is defined as the study of equilibrium states of a 
system which has been subjected to some energy transformation. 
More specifically, thermodynamics is concerned with transforma- 
tions of heat into mechanical work and of mechanical work into 
heat. 


e A system is a specific sample of matter. In the atmosphere a par- 
cel of air is a system. A system is called open when it exchanges 
matter and energy with its surroundings (Figure 1.1). In the 
atmosphere all systems are more or less open. A closed system is 
a system that does not exchange matter with its surroundings. 
In this case, the system is always composed of the same point- 
masses (a point-mass refers to a very small object, for example 
a molecule). Obviously, the mathematical treatment of closed 
systems is not as involved as the one for open systems, which 
are extremely hard to handle. Because of that, in atmospheric 
thermodynamics, we assume that most systems are closed. This 
assumption is justified when the interactions associated with 
open systems can be neglected. This is approximately true in the 
following cases. (a) The system is large enough to ignore mixing 
with its surroundings at the boundaries. For example, a large 
cumulonimbus cloud may be considered as a closed system but a 
small cumulus may not. (b) The system is part of a larger homo- 
geneous system. In this case mixing does not significantly change 
its composition. A system is called isolated when it exchanges 
neither matter nor energy with its surroundings. 


e The state of a system (in classical mechanics) is completely speci- 
fied at a given time if the position and velocity of each point-mass 
is known. Thus, in a three-dimensional world, for a system of N 
point-masses, 6N variables need to be known at any time. When 


Figure I.1 

In an open system mass 
and energy can be 
exchanged with its 
environment. A system is 
defined as closed when it 
exchanges energy but not 
matter with its 
environment, and as 
isolated if it exchanges 
neither mass nor energy. 
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Open system 





N is very large (like in any parcel of air) this dynamical defini- 
tion of state is not practical. As such, in thermodynamics we are 
dealing with the average properties of the system. 


If the system is a homogeneous fluid consisting of just one com- 
ponent, then its thermodynamic state can be defined by its 
geometry, by its temperature, T, and pressure, p. The geometry 
of a system is defined by its volume, V, and its shape. How- 
ever, most thermodynamic properties do not depend on shape. 
As such, volume is the only variable needed to characterize geom- 
etry. Since p, V, and T determine the state of the system, they 
must be connected. Their functional relationship f(p,V,T) = 0 
is called the equation of state. Accordingly, any one of these vari- 
ables can be expressed as a function of the other two. It follows 
that the state of a one-component homogeneous system can be 
completely defined by any two of the three state variables. This 
provides an easy way to visualize the evolution of such a sys- 
tem by simply plotting V against p in a rectangular coordinate 
system. In such a system, states of equal temperature define an 
isotherm. Any other thermodynamic variables that depend on the 
state defined by the two independent state variables are called 
state functions. State functions are thus dependent variables and 
state variables are independent variables; the two do not differ 
in other respects. That is why in the literature there is hardly 
any distinction between state variables and state functions. State 
variables and state functions have the property that their changes 
depend only on the initial and final states, not on the particular 
way by which the change happened. If the system is composed of 
a homogeneous mixture of several components, then in order to 
define the state of the system we need, in addition to p, V, 7, the 
concentrations of the different components. If the system is non- 
homogeneous, we must divide it into a number of homogeneous 
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parts. In this case, p, V, and T of a given homogeneous part are 
connected via an equation of state. 

For a closed system, the chemical composition and its mass 
describe the system itself. Its volume, pressure, and tempera- 
ture describe the state of system. Properties of the system are 
referred to as extensive if they depend on the size of the sys- 
tem and as intensive if they are independent of the size of the 
system. An extensive variable can be converted into an intensive 
one by dividing by the mass. In the literature it is common to 
use capital letters to describe quantities that depend on mass 
(work, W, entropy, S) and lower case letters to describe inten- 
sive variables (specific work, w, specific heat, q). The mass, m, 
and temperature, T, will be exceptions to this rule. 


An equilibrium state is defined as a state in which the system’s 
properties, so long as the external conditions (surroundings) 
remain unchanged, do not change in time. For example, a gas 
enclosed in a container of constant volume is in equilibrium if 
its pressure is constant throughout and its temperature is equal 
to that of the surroundings. An equilibrium state can be stable, 
unstable, or metastable. It is stable when small variations about 
the equilibrium state do not take the system away from the equi- 
librium state, and it is unstable if they do. An equilibrium state 
is called metastable if the system is stable with respect to small 
variations in certain properties and unstable with respect to small 
changes in other properties. 


A transformation takes a system from an initial state 7 to a final 
state f. In a (p,V) diagram such a transformation will be rep- 
resented by a curve J connecting i and f. We will denote this 
as i = f. A transformation can be reversible or irreversible. 
Formally, a reversible transformation is one in which the succes- 
sive states (those between i and f) differ by infinitesimals from 
equilibrium states. Accordingly, a reversible transformation can 
only connect those 7 and f states which are equilibrium states. 
It follows that a reversible process is one which can be reversed 
anywhere along its path in such a way that both the system 
and its surroundings return to their initial states. In practice a 
reversible transformation is realized only when the external con- 
ditions change very slowly so that the system has time to adjust 
to the new conditions. For example, assume that our system is 
a gas enclosed in a container with a movable piston. As long as 
the piston moves from 2 to f very slowly the system adjusts and 
all intermediate states are equilibrium states. If the piston does 
not move slowly, then currents will be created in the expanding 
gas and the intermediate states will not be equilibrium states. 
From this example, it follows that turbulent mixing in the atmo- 
sphere is a source of irreversibility. If a system goes from i to 


Figure 1.2 

(a) A reversible cyclic 
process; (b) an 
irreversible cyclic process. 
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f reversibly, then it could go from f to i again reversibly if the 
same steps were followed backwards. If the same steps cannot 
be followed exactly, then this transformation is represented by 
another curve I’ in the (p, V) diagram (i.e. f a i) and may or 
may not be reversible. In other words the system may return to 
its initial state but the surroundings may not. Any transforma- 
tion i — f —> i is called a cyclic transformation. Given the 
discussion above we can have cyclic tansformations which are 
reversible or irreversible (Figure 1.2). A transformation i — f 
is called isothermal if J is an isotherm, isochoric if I is a constant 
volume line, isobaric if J is a constant pressure line, and adiabatic 
if during the transformation the system does not exchange heat 
with its surroundings (environment). Note and keep in mind for 
later that adiabatic transformations are not isothermal. 


Energy is something that can be defined formally (we have to 
wait a bit for this), but its concept is not easily understood by 
defining it. We all feel we understand what is meant by energy, 
but if we verbally attempt to explain what energy is we will get 
upset with ourselves. At this point, let us just recall that for 
a point-mass with a mass m, moving with speed v in a uni- 
form gravitational field g, Newton’s second law takes the form 


l= reversible 
l’ = reversible 








M 





(a) y 
. f 
p | 
r 
i 
l= reversible 
I’ = irreversible 





(b) v 
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d(K + P)/dt = 0 where K = m,v?/2, P = m,gz, tis the 
time, and z is the height. K is called the kinetic energy and P is 
called the potential energy. The total energy of the point-mass 
E = K + P is, therefore, conserved. If we consider a system of 
N interacting point-masses that may be subjected to external 
forces (other than gravity), then the total energy of the system 
is the sum of the kinetic energy about the centre of gravity of 
all point-masses (internal kinetic energy), the kinetic energy of 
the centre of gravity, the potential energy due to interactions 
between the point-masses (internal potential energy), and the 
potential energy due to external forces. The sum of the internal 
kinetic and internal potential energy is called the internal energy 
of the system, U. A system is called conservative if dU/dt = 0 
and dissipative otherwise. For systems considered here we are 
interested in the internal energies only. 





CHAPTER TWO 





Some useful mathematical 
and physical topics 


2.1 Exact differentials 


If z is a function of the variables x and y, then by definition 


Oz Oz 
dz = (a) a+ (Z) a (2.1) 


where dz is an exact differential. Now let us assume that a quan- 
tity dz can be expressed according to the following differential 
relationship 

ôz = Mdx + Ndy (2.2) 


where x and y are independent variables and M and N are functions 
of x and y. If we integrate equation (2.2) we have that 


fsz= [Macs | Nay. 


Since M and N are functions of x and y, the above integration 
cannot be done unless a functional relationship f(x, y) = 0 between 
x and y is chosen. This relationship defines a path in the (x,y) 
domain along which the integration will be performed. This is called 
a line integral and its result depends entirely on the prescribed path 
in the (x,y) domain. If it is that 


Oz Oz 
M= F N= a (2.3) 
then equation (2.2) becomes 
Oz Oz 


The right-hand side of the above equation is the exact or total 
differential dz. In this case dz is an exact differential. If we 
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now integrate dz from some initial state i to a final state f we 
obtain 


a ôz = f dz = 2(x5, Ys) — 2(is Yi). (2.4) 


Clearly, if 6z is an exact differential its net change along a path 
i —> f depends only on points 7 and f and not on the particular 
path from 7 to f. We say that in this case z is a point function. 
All three state variables are exact differentials (i.e. dp = dp, ôT = 
dT, V = dV). It follows that all quantities that are a function of 
any two state variables will be exact differentials. 

If the final and initial states coincide (i.e. we go back to the initial 
state via a cyclic process), then from equation (2.4) we have that 


f z =0. (2.5) 


The above alternative condition indicates that 6z is an exact differ- 
ential if its integral along any closed path is zero. At this point we 
have to clarify a point so that we do not get confused later. When we 
deal with pure mathematical functions our ability to evaluate $ ôz 
does not depend on the direction of the closed path, i.e. whether 
we go from i back to i via i = f b ior viai D fi 
(see Figure 1.2). When we deal with natural systems we must view 
the condition $ ôz = 0 in relation to reversible and irreversible 
processes. If, somehow, it is possible to go from 7 back to 7 via 
jeij 5 i but impossible via 7 ahs f i (for example, 
when J’ is an irreversible transformation), then computation of 6z 
depends on the direction and as such it is not unique. Therefore, 
for physical systems, the condition $ 6z = 0 when 6z is an exact 
differential applies only to reversible processes. 

Note that since 


ð Oz o?z Or z ð Oz 
OyOx OyOx OxOy Ox Oy 


it follows that an equivalent condition for dz to be an exact 
differential is that 











OM ON 

Oy Ox” 
Equations (2.3)—(2.6) are equivalent conditions that define z as a 
point function and subsequently 6z as an exact differential. If a 
thermodynamic quantity is not a point function or an exact differ- 
ential then its change along a path depends on the path. Moreover, 
its change along a closed path is not zero. Such quantities are path 
functions. For a path function the thermodynamic processes must 
be specified completely in order to define the quantity. For the rest 
of this book an exact differential will be denoted by dz whereas a 
non-exact differential will be denoted by 6z. Finally, note that if z 


(2.6) 
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is not an exact differential and if only two variables are involved, a 
factor A (called the integration factor) may exist such that Adz is 
an exact differential. 


2.2 Kinetic theory of heat 


Let us consider a system at a temperature T, consisting of N 
point-masses (molecules). According to the kinetic theory of heat, 
these molecules move randomly at all directions traversing rectilin- 
ear lines. This motion is called Brownian motion. Because of the 
complete randomness of this motion, the internal energies of the 
point-masses not only are not equal to each other, but they change 
in time. If, however, we calculate the mean internal energy, we will 
find that it remains constant in time. The kinetic theory of heat 
accepts that the mean internal energy of each point-mass, U, is 
proportional to the absolute temperature of the system (a formal 
definition of absolute temperature will come later; for now let us 
denote it by T), 

U = constant x T (2.7) 
Let us for a minute assume that N = 1. Then the point has 
only three degrees of freedom which here are called thermodynamic 
degrees of freedom and are equal to the number of independent vari- 
ables needed to completely define the energy of the point (unlike 
the degrees of freedom in Hamiltonian dynamics which are defined 
as the least number of independent variables that completely define 
the position of the point in state space). The velocity v of the point 
can be written as 

v? =v +o + o. 


Because we only assumed one point, then the total internal energy 
is equal to its kinetic energy. Thus, 








y= Mpv? 
2 
or > ; 3 
U = eie. Uae Ü= e 
2 2 2 


where to each component corresponds one degree of freedom. 
According to the equal distribution of energy theorem, the mean 
kinetic energy of the point, U, is distributed equally to the three 
degrees of freedom i.e. U, = U, = U,. Accordingly, from equation 
(2.7) we can write that 


VU, AT. 4a 2s 


where the constant A is a universal constant (i.e. it does not depend 
on the degrees of freedom or the type of the gas). We denote this 
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constant as k/2 where k is Boltzmann’s constant (k = 1.38 x 107? 
J K~") (for a review of units, see Table A1 in the Appendix). There- 
fore, the mean kinetic energy of a point with three degrees of 
freedom is equal to 





= kT 
U = ont (2.8) 
2 
or 
Mpv? _ 3kT 
ZB 


The theorem of equal energy distribution can be extended to N 
points. In this case, the degrees of freedom are 3N and 








D3 Mpv? _ a 


or 





S 





PE 


(25y = : (2.9) 


where (mp) is the average kinetic energy of all N points. Note that 
the above is true only if the points are considered as monatomic. 
If they are not, extra degrees of freedom are present that corre- 
spond to other motions such as rotation about the center of gravity, 
oscillation about the equilibrium positions, etc. 

The kinetic theory of heat has found many applications in the 
kinetic theory of ideal gases. An ideal gas is one for which the 
following apply: 


vl 


or 








(a) the molecules move randomly in all directions and in such a 
way that the same number of molecules move in any direction; 

(b) during the motion the molecules do not exert forces except 
when they collide with each other or with the walls of the 
container. As such the motion of each molecule between two 
collisions is linear and of uniform speed; 

(c) the collisions between molecules are considered elastic. This 
is necessary because otherwise with each collision the kinetic 
energy of the molecules will be reduced thereby resulting in a 
temperature decrease. Also, a collision obeys the law of specular 
reflection (the angle of incidence equals the angle of reflection); 

(d) the sum of the volumes of the molecules is negligible compared 
with the volume of the container. 


Now let us consider a molecule of mass m, whose velocity is v and 
which is moving in a direction perpendicular to a wall (Figure 2.1). 
The molecule has a momentum P = m,v. Since we accept that the 
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Figure 2.1 
A molecule of mass mp V 
moving with a velocity v 
and hitting a surface S. If 
this collision is assumed 
elastic and specular, then 
the change in momentum 
is 2mpv. 














collision is elastic and specular, the magnitude of the momentum 
after the collision is —m,v. Thus the total change in momentum is 
mpu — (mpv) = 2m,v. According to Newton’s second law, F = 
dP/dt. If we consider N molecules occupying a volume V we can 
calculate the change in momentum dP of all molecules in the time 
interval dt by multiplying the change in momentum of one molecule 
(2m,v) by the number of molecules dN that hit a given area S' on 
the wall, i.e. 


dP = 2m,vdN (2.10) 


Note that here we have assumed that all molecules have the same 
speed. The number dN of molecules hitting area S during dt is equal 
to the number of molecules which move to the right and which are 
included in a box with base S and length vdt. Since the motion 
is completely random, we can assume that i of the molecules will 
be moving to the right, will be moving to the left, and 4 will be 
moving along the directions of the other two coordinates. Since the 
volume of the box is Svdt and the number of molecules per unit 


volume is N/V then the number of molecules inside the box is 


N 


Accordingly, the number of molecules moving to the right and 
colliding with S is 


N 
dN = —Svdt 2.11 
6v" ae, 
From equations (2.10) and (2.11) it follows that 


N 
dP = 2mp_v" gy ot 
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Recalling the definition of pressure, p (pressure = force/area), and 
Newton’s second law we obtain 


1N a 
p= 3 pie’ ; 
The above formula resulted by assuming that all molecules move 
with the same speed. This is not true, and because of that mpv? in 
the above equation should be replaced by the average of all points, 
(m,v?). We thus arrive at 


1N 

ae, 

It can easily be shown that by combining equations (2.9) and (2.12), 
we can derive an equation that includes all three state variables: 








(m,v?). (2.12) 


pV = NKT. (2.13) 


Equation (2.13) provides the functional relationship of the equation 
of state f(p,V,T) and it is called the ideal gas law. More details 
follow in the next chapter. 


CHAPTER THREE 





Early experiments and laws 


At the end of Chapter 2 we derived theoretically the equation of 
state or the ideal gas law. This law was first derived experimentally. 
The relevant experiments provide many interesting insights about 
the properties of ideal gases and confirm the theory. As such a little 
discussion is necessary. 


3.1 The first law of Gay-Lussac 


Through experiments Gay-Lussac was able to show that, when pres- 
sure is constant, the increase in volume of an ideal gas, dV, is 
proportional to the volume Vo that it has at a temperature (mea- 
sured in the Celsius scale) of 6 = 0°C and proportional to the 
temperature increase, d0: 


dV = aVodd. (3.1) 


The coefficient a is called the volume coefficient of thermal expan- 
sion at a constant pressure and it has the value of 1/273 deg~* for 
all gases. The physical meaning of œ can be understood if we solve 
equation (3.1) for a: 
1 dV 
a= a Te 
From the above equation it follows that if we increase the temper- 
ature of an ideal gas by 1°C, while we keep the pressure constant, 
the volume will increase by 1/273 of the volume the gas occupies 
at 0°C. By integrating equation (3.1) we obtain the relationship 


between V and 8: P f 
f dV = J aVodé 
Vo 0 


V — Vo = avd 


or 
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Figure 3.1 Vt 
Graphical representation 
of the first law of 
Gay-Lussac. 


p=constant 








= = (9 > 
-273 0 6 (°C) 


or 
V = V(1 + a8). (3.2) 


This is a linear relationship and its graph is shown in Figure 3.1. 


3.2 The second law of Gay-Lussac 


Again through experimentation Gay-Lussac was able to show that 
for a constant volume the increase in pressure of an ideal gas, dp, 
is proportional to the pressure pọ that it has at a temperature of 
0°C and proportional to the increase in temperature dé: 


dp = Bpodé. 


The coefficient 8 is the pressure coefficient of thermal expansion at 
a constant volume and has the value 1/273 deg for all gases: 


go”. 

Po 
The above formula indicates that an increase in temperature by 
1°C (while V is constant) results in an increase in pressure by 
1/273 of the pressure the gas had at 0°C. 
As in the first law again it follows that (Figure 3.2) 


p = po(1 + b0). (3.3) 


Application 


The second law of Gay-Lussac can easily explain why, in the winter, 
heating a house at a much greater temperature than the outside 
temperature does not increase the pressure enough to break the 
windows. A difference of 20°C between the inside and the out- 
side air increases the pressure inside the house by 7.3%. Glass can 
withstand such pressure changes easily. 


Figure 3.2 

Graphical representation 
of the second law of 
Gay-Lussac. 
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V = constant 








= = 6" > 
-273 0 6 (°C) 


3.3 Absolute temperature 


From equation (3.2) it follows that if we extrapolate to 0=—273 °C, 
then V=0. This means that if we were able to cool an ideal gas 
to —273 °C, while keeping the pressure constant, the volume would 
become zero. Similarly, from equation (3.3) it follows that if we 
were able to cool an ideal gas to the temperature of —273 °C, while 
keeping the volume constant, the pressure would become zero. This 
temperature of —273°C we call absolute zero. Up to now for the 
measurement of temperature the Celsius scale has been used, start- 
ing from the temperature “zero Celsius”. If we extend the Celsius 
scale to the absolute zero (i.e. —273 °C), then the temperature mea- 
sured from the absolute zero is called absolute temperature, T. This 
defines a new scale called the Kelvin scale (T = 273 + 0). 


3.4 Another form of the Gay-Lussac laws 


Using the absolute temperature we can present the Gay-Lussac laws 
as follows. From Figure 3.3 it follows from the similarity of triangles 
ABC and AB’C’ that his first law can be expressed as 


T 


for p = constant, yop 


(3.4) 


Similarly from Figure 3.4, his second law can be expressed as 


T 
for V = constant, ae (3.5) 
p T’ 
Verbally stated, the volumes of an ideal gas under constant pressure 
and the pressures under constant volume are proportional to the 
absolute temperatures. 
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Figure 3.3 

Expressing Gay-Lussac’s 
first law for two states C 
and C”. 


Figure 3.4 

Expressing Gay-Lussac’s 
second law for two states 
C and C”. 
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3.5 Boyle’s law 


While Gay-Lussac’s laws provide the change in volume or pressure 
as a function of temperature, Boyle’s law provides the change in 
pressure as the volume varies at a constant temperature. The law 
is expressed as follows: 

pV = pV". (3.6) 


3.6 Avogadro’s hypothesis 


The formal definition of a mole is the amount of substance which 
contains the same number of particles (atoms, molecules, ions, or 
electrons) as there are in 12 grams of !C. This number of atoms is 
equal to N = 6.022 x 10” and it is known as Avogadro’s number. 
It is obtained under the condition that one mole of carbon weighs 
12 grams. In general, a substance of a weight equal to its molecular 
weight contains one mole of the substance. Thus 25 grams of water 
contain 27/18 = 1.5 mole of water. In 1811 the Italian physicist 
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and mathematician Amedeo Avogadro proposed his hypothesis that 
the volume of a gas, V, is directly proportional to the number of 
molecules of the gas, N’, 


V =aN' 


where a is a constant. To relate to this hypothesis, simply imagine 
inflating a balloon. The more the air you pump into it, the bigger 
the volume. It follows that at constant temperature and pressure, 
equal volumes of gases contain the same number of molecules. For 
one mole N’ = N. Since N is the same for one mole of any gas, 
Avogadro’s hypothesis (actually a law by now) can be stated simply 
as: A mole of any gas at constant temperature and pressure occupies 
the same volume. This volume for the standard state Tọ = 0°C, 
Po = 1 atmosphere has the value (see equation (2.13)) 


Vrp = 22-4 liters mol™* 
= 22400cm* mol”. 


3.7 The ideal gas law 


Let us consider an ideal gas at a state p, V, T which is heated under 
constant volume to a state pı, V, T” (Figure 3.5). Then, according 
to Gay-Lussac’s second law 
T". 
Pi = Pr: 


If subsequently we keep the temperature constant and increase the 
volume to V’ the gas goes to a state p'V'T’. Then, according to 
Boyle’s law, 

pP'V' = pV. 


By combining the above two equations we obtain 


T 





pV =p 
or 
pV pV! 
-em ; 3.7 
TT (3.7) 


This is called the Boyle-Gay-Lussac law and it indicates that if we 
change the pressure, the volume, and the temperature the value of 
PV/T remains constant. Thus, equation (3.7) can be written as 


pV = AT (3.8) 


where A is, according to equation (3.7), a constant depending on 
the type and mass of the gas. The dependence on the mass follows 
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Figure 3.5 p 

The three steps used to 

derive the ideal gas law 

experimentally. P V,T’ 


yer Isotherm 





from the fact that if we consider twice the mass and keep the pres- 
sure and the temperature constant, then the volume will double 
(Avogadro’s hypothesis). Therefore, the value of pV/T (i.e. A) will 
double. Accordingly, we can write A = mR where m is the mass of 
the gas and R is a constant independent of the mass but dependent 
on the gas. This constant is called the specific gas constant. We can 
then write equation (3.8) as 


pV = mRT. (3.9) 


Since m = Nm, (where N is the total number of molecules each 
having mass m,), equation (3.9) is identical to equation (2.13) with 
R = k/m,. Equation (3.9) is thus the equation of state or ideal 
gas law which is now derived experimentally. Equation (3.9) can be 
written as p = pRT where p = m/V is the mass density. By defining 
the specific volume a = 1/p, the ideal gas law takes the form 


pa = RT. 


If M denotes the molecular weight of the gas then the number of 
moles is n = m/M. It follows that 


pV =nMRT 
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or 
pV =nR*T (3.10) 
or 
R* 
= T 
pa M 


where R* = MR. Considering that R = k/m, it follows that 
R* = Mk/m, = mk/nm, = Nk/n. Now recall that the number 
of molecules in one mole, N4, is equal to 6.022 x 107% (Avogadro’s 
number). Then, N = nN. In this case R* = Nak, which is the 
product of two constants. This new constant is called the universal 
gas constant and its value is 8.3143 JK~'mol"}. 


3.8 A little discussion on the ideal gas law 


e Because the ideal gas law relates three variables (rather than 
two), we have to be careful when we interpret changes in one of 
the variables. For example, temperature increases when pressure 
increases but only if the volume (density) increases (decreases) 
or remains constant or decreases (increases) by a smaller amount 
compared with pressure. Let us consider the ideal gas law pV = 
mRT or p = pRT and differentiate it: 


pdV + Vdp = mRdT 
dp = RTdp + RpdT. 


We can easily see that if dp > 0 and dV > 0 then dT > 0. Also, if 
dp > 0 and dV < 0 and |Vdp| > |pdV| then dT > 0. And so on! 
It follows that cold air is denser than warm air only if pressure 
remains constant or if its change does not offset the temperature 
difference. 

Similarly, if we consider the ideal gas law in the form 








R* 
= T 
pa M 
or 
= T 
P=p M 


we can reason that moist air (in which some of the dry air of 
molecular weight 28 g mol~' has been replaced by water vapor of 
molecular weight 18gmol~') is less dense than dry air at more 
or less the same pressure and temperature. This explains some 
interesting statistics in the American game of baseball where 
more home runs occur when the weather is hot and humid. In 
a warmer and more moist environment where p = constant, the 
density of air is smaller. Therefore, the ball has less resistance. 
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For a path of 400 feet or so the effect can be significant, resulting 
in higher chances for a home run. 


e Recall equations (2.9) and (3.2). According to equation (2.9) 
T = 0K when v = 0. Because of this one might interpret the 
zero absolute temperature as the temperature where all motion 
in an ideal gas ceases. On the other hand, when equation (3.2) 
is extrapolated to the temperature at which V = 0 it results 
in Tabsolute = (—+)°C = —273°C = OK. While this tempera- 
ture corresponds to vanishing volumes, it is not necessarily the 
temperature at which motion ceases. If both V = 0 and v = 0 
are true when T = OK, then pressure can be anything, but this 
does not invalidate the ideal gas law, both sides of which will 
equal zero. To avoid confusion, we must keep in mind that these 
issues emerge when these equations are extrapolated to states 
where one cannot assume that a gas behaves as an ideal gas. 
At those states, equations (2.9) and (3.2) simply do not mean 
much. 


3.9 Mixture of gases — Dalton’s law 


Consider a mixture of two gases occupying a volume V and con- 
sisting of N, molecules of gas 1 and N, molecules of gas 2. The 
total pressure on the walls will be the result of all the collisions, 
i.e. the collisions by the molecules of gas 1 and the collisions by the 
molecules of gas 2. Thus we can write equation (2.13) as 

(N: + N2) 


= atr 
$ V 


or 








_MkE | NokT 
© V V ` 


The first term on the right-hand side of the above equation is 
exactly the pressure that gas 1 would create if all its N; molecules 
occupied the volume V, i.e. the partial pressure of gas 1. The same 
is valid for the second term. We thus conclude that the total pres- 
sure is the sum of the partial pressures. This expresses Dalton’s 
law which states that for a mixture of K components, each one of 
which obeys the ideal gas law, the total pressure, p, exerted by the 
mixture is equal to the sum of the partial pressures which would 
be exerted by each gas if it alone occupied the entire volume at the 
temperature of the mixture, T: 
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If the volume of the mixture is V and the mass and molecular 
weight of the ith constituent are m; and M; respectively, then for 
each constituent 

S R* m; 

=V 





Pi T. 


Applying Dalton’s law, we then can write that for a mixture 


X R*T mi 
a? EI 




















or 
RT Am, 
ao 
Us i=l M; 
Since the total mass of the mixture m = DS m,;, the above 
equation becomes 
K mi 
_ R*Tm Des Mi 
” V D nh 
or 
K mi 
_ ÆT aia m 


S (3.11) 


a DF Mi 


In order for the mixture to obey the ideal gas law, it must satisfy 
the equation 





RT RT 
SSS ee =— 3.12 
Pog ae (3.12) 





where M is the mean molecular weight of the mixture. By compar- 
ing equations (3.11) and (3.12) we see that this would be possible 
if 


K 
M = Se (3.13) 
Diz m; 
This provides the proper way to compute the mean molecular 
weight for a mixture and indicates that the mixture follows the 
ideal gas law as long as there is no condensation (if there is, then 
some m;’s may not remain constant). 


Examples 


(3.1) Determine the mean molecular weight of dry air. 


For our planet, the lowest 25km of the atmosphere is made 
up almost entirely by nitrogen (N3), oxygen (O3), argon (A), 
and carbon dioxide (CO2) (75.51%, 23.14%, 1.3%, and 0.05% 
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(3.2) 


by mass, respectively). Thus, for dry air the mean molecular 
weight is 


az _ MN T Mo, TMa T Mco 
M = 2 2 2 























mN | Oo ma MCOg 
Mya Moz Ma Mcog 
or 
s 75.51 + 23.16 + 1.3 + 0.05 i 
M = B3 1 214 13 | 20 g mol 
28.02 ' 32.0 T 39.94 ' 44.01 
or 
M = 28.97 g mol`! 
or 


M = 0.02897 kg mol™*. 
It follows that the specific gas constant for dry air is 


Ra = R*/M _= 287 Jkg K. 


Determine the mean molecular weight of a mixture of dry air 
saturated with water vapor at 0°C and | atmosphere pressure. 
The partial pressure of water vapor at 0°C is 6.1 | mb. 


The mixture consists of dry air and water vapor. If pa is the 
pressure due to dry air and p, the pressure due to water 
vapor, then the pressure of the mixture is p = pa + py. If 
we recall that the number of moles n = m/M, then we can 
write equation (3.13) as 





IF BE _ = Ea! ree ` 
M = a= = 3.14 
sae =e ee 


where n; and n are the corresponding number of moles of the 
constituents of the mixture and the total number of moles 
in the mixture. Assuming that both dry air and water vapor 
are ideal gases we have that 


pV = nR*T for the mixture 
pV =n,R*T for each of the constituents. 


From these two equations it follows that 


Combining the above relationship and equation (3.14) yields 


Poi 


<P 
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(3.3) 


In our case i = 1,2. Thus, 
M = Pi wr + PM, 
Pp p 
or changing 1 > d (dry air) and 2 > v (water vapor) 


M =P M+ EM, 
p p 


or 





_ A A 
Pp 


M M,+ 2M, 
Pp 
or 
M =28.9¢mol". 


The above estimated mean molecular weight is not very 
much different from the mean molecular weight of dry air. 
The difference becomes more apparent for higher temper- 
atures. For example for T = 35°C the partial pressure of 
water vapor is 57.6mb and M =28.3gmol' . 


Two containers A and B of volumes Va = 800cm? and Vg = 
600 cm’, respectively, are connected with a tube that closes and 
opens by a hinge. The containers are filled with a gas under 
pressures of 1000 mb and 800 mb, respectively. If we open the 
connection what will the final pressure be in each of the con- 
tainers? Assume that the temperature remains constant. 


Once the connection is open, each gas will expand to fill the 
total volume V = Va +V thereby equalizing the difference in 
pressure between the two containers. Thus the final pressure 
in each container will be the same. Let us first assume that 
container B is empty. Since the gas in container A expands 
at a constant temperature we have (Boyle’s law) 


PVF = DiVi 


where į and f stand for initial and final. It follows that 
py = 571 mb. 

Similarly if we assume that container A is empty the final 
pressure in both containers will be 





1 DiVi 
p, = Pit = 343 mb. 
f V; 


Since ps(p) is the pressure the gas in A(B) would exert if it 
occupied the total volume Va + Vg, py and p, can be consid- 
ered as partial pressures of two gases. Then from Dalton’s 
law it follows that the final pressure in each container should 
be 914 mb. 
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(3.4) 


If Boyle had observed ,/pV = constant, what would the equa- 
tion of state for an ideal gas be? In this case calculate the 
temperature of a sample of nitrogen which has a pressure of 
800 mb and a specific volume of 1200 cm? g~. How does this 
temperature differ from that estimated when the correct law 
pV = constant is observed? Can you explain the difference? 


If we follow the procedure to arrive at equation (3.7) but 
with the new Boyle’s law, we have that 


T’ 
pPr= Pr 


VP'V' = Jn. 


By combining the above equations we have that 


; [ 1 
Jp'V = VP —V 
or 


VPV’ _ BV 


= —— = constant. 


v vT 





As we know, if we keep the temperature and pressure con- 
stant and double the mass, the volume will double. As such, 
here as well the ratio \/pV/ VT will double. Then we can 
write that \/pV/ VT = mR where m is the total mass and 
R is the specific gas constant. Thus, the ideal gas law in this 
case will be 





Solving for T we have that 
VT = 1.143 


or 


T= 1.07K. 


Under normal circumstances (i.e. when pa = (R*/M)T), 
we find that T = 323.3K which makes more sense. The 
tremendous difference is due to the fact that the square root 
in ,/pa = constant introduces two corrections to the ideal gas 
law (both pressure and temperature appear under a square 
root) having the net result of significantly reducing T. 
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Problems 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


What is the mass of dry air occupying a room of dimensions 
3x5 x 4m at p=1 atmosphere and T = 20°C? (72.3kg) 
Graph the relationship V = f(T), = constant from absolute 
zero up to high temperatures, for two samples of the same 
gas which at 0°C occupy volumes of 1000 and 2000cm?, 
respectively. 

In a 2-D coordinate system with axes the absolute temper- 
ature and volume, graph (a) an isobaric (constant pressure) 
change of one mole of an ideal gas for p = 1 atmosphere, 
and (b) the same when p = 2 atmospheres. 

Determine the molecular weight of the Venusian atmosphere 
assuming that it consists of 95% CO, and 5% N, by volume. 
What is the gas constant for 1kg of such an atmosphere? 
(43.2 gmol-?, 192.5 Jkg-! K-?) 

If p = 1 atmosphere and T = 0°C how many molecules are 
there in 1cm? of dry air? (2.6884 x 10'° molecules) 

Given the two states p, V,T and p’,V’,T’, define on a (p, V) 
diagram the state p,,V,7”’ that was used to arrive at the 
ideal gas law (section 3.7). 

An ideal gas of p,V,T undergoes the following successive 
changes: (1) it is warmed under a constant pressure until 
its volume doubles, (2) it is warmed under constant volume 
until its pressure doubles, and (3) it expands isothermally 
until its pressure returns to p. Calculate in each case the 
values of p,V,T and plot all three changes in a (p,V) 
diagram. 


CHAPTER FOUR 





The first law of 
thermodynamics 


4.1 Work 


As we have already mentioned, in atmospheric thermodynamics we 
will be dealing with equilibrium states of air. If a system (par- 
cel of air, for example) is at equilibrium with its environment no 
changes take place in either of them. We can imagine the “shape” 
of the parcel remaining unchanged in time. If the pressure of the 
surroundings changes, then the force associated with the pressure 
change will disturb the parcel thereby forcing it away from equilib- 
rium. In order for the parcel to adjust to the pressure changes of 
the surroundings, the parcel will either contract or expand. If the 
parcel expands we say that the parcel performs work on the envi- 
ronment and if the parcel contracts we say that the environment 
performs work on the parcel. By definition, if the volume change is 
dV then the incremental work done, dW, is 


dW = pdV. 


Accordingly, when the system changes from an initial state 7 to a 
final state f the total work done, either by the system or on the 


system, is 
f 
W = / pdV. 


The above equation indicates that the work done is given by an area 
in a (p, V) diagram (Figure 4.1) [or a (p,a) diagram, where a is the 
specific volume]. If dV > 0 (the system expands) it follows that 
W > 0 and if dV < 0 (system contracts) it follows that W < 0. 
Thus, positive work corresponds to work done by the system on 
the environment and negative work corresponds to work done to 
the system by the environment. 

Now let us consider a situation where the system expands 
through a reversible transformation from i to f and then contracts 
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Figure 4.1 

The shaded area gives the 
work done when a system 
changes from an initial 
state i to a final state f. 


Figure 4.2 

The shaded area gives the 
work done during a cyclic 
reversible transformation. 
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PF f 
f 

i we | pdV 
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V; V; V 
pt 

f 
l W=Anfai 
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y, V; V 


from f to i along exactly the same path in the (p,V) diagram. 
Then, the total work done will be 


f i f f 
w=) pdv + f pav = f pav — f pdV = 0. 
i f i i 


If, however, the system expands from 7 to f and then contracts 
from f to i but along a different reversible transformation (different 
path), then the total work done would be (see Figure 4.2) 


[f vav] + [fra] 


= area under curve 1 — area under curve 2 
= Ay f2i # 0 


where Aj1,2; is the area enclosed by the two paths. It follows that 
$ dW = $ pdV £ 0 and thus dW is not an exact differential, which 
means that work is not a state function. As such it depends on the 
particular way the system goes from i to f. Because of this from 
now on we will denote the incremental change in work as dW not 
as dW. 


W = 
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If we write dV = dAds (where dA is an area element and ds a 
distance element) we have that DW = pdAds = Fds = Fvdt or 
that 


sw 
CRAN 
dt £ 
or 
ôW dv 
a VH 
or 
W d(l 
dt dt \2 
K SW dK 
ase (4.1) 


where v denotes the velocity of the parcel and K is the kinetic 
energy of the parcel. It follows that the work done and the kinetic 
energy are related. The last equation indicates that one way by 
which a thermodynamic system can exchange energy with its envi- 
ronment is by performing work. The other is through transfer of 
heat. More on that will follow soon. According to (4.1) the units 
for work are those of energy. Thus, the unit for work in the MKS 
system is the joule which is defined as J = Nm where the newton 
N = kgm s7?. In the cgs system, the unit is the erg which is defined 
as erg = dyncm where dyn = gcms™?. It follows that 1 joule = 
10” erg. 


4.2 Definition of energy 


The first law of thermodynamics expresses the principle of con- 
servation of energy for thermodynamical systems. The idea here 
is that energy cannot be created or destroyed. It can only change 
from one form to another. In this regard if during a transforma- 
tion the energy of the system increases by some amount, then this 
amount is equal to the amount of energy the system receives from 
its surroundings in some other form. 

Let us consider a closed system contained in an adiabatic enclosure. 
In this case the energy of the system, U, is equal to the sum of the 
potential and kinetic energy of all its molecules. The sum of the 
energies of all molecules depends on the state of the system at a 
given moment (i.e. on the values p, V, T) but obviously is indepen- 
dent of past states. It follows that the internal energy of the system 
depends on the state in which it exists but not on the way it arrived 
at that state. Thus, in a transformation i > f, AU = Up — Uj. 
Then, for a cyclic transformation, ¢ AU = 0 which means that for 
an infinitesimal process dU is an exact differential. If no external 
forces are acting upon the system (i.e. the system is at equilibrium 
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Figure 4.3 

The work done by an 
external force on a 
system that is taken 
adiabatically from a 
reference state O ora 
reference state O’ toa 
state A. In this case the 
work done depends on 
the initial and final states, 
not on the particular path 
from an initial state to a 
final state. 
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Moy. 





awe 


with its environment), then its energy remains the same (AU = 0). 
If some external force acts upon the system taking it from state i 
to state f then the conservation of energy principle dictates that 


U; — U; = -W° (4.2) 
where —W°¢4 is the work done adiabatically by the external force 
on the system (+W7?? is the work done by the system). Note that 
equation (4.2) implies that the work done in this case depends only 
on states 7 and f, not on the particular transformation i > f. This 
is contrary to what we proved previously — that the work done 
is not an exact differential and that it depends on the particular 
transformation i + f. The discrepancy is due to the fact that here 
we are dealing with adiabatic transformations. Only in this case 
does the work done not depend on the particular transformation. 
Assume now a standard state O of this system where by definition 
Uo =0 and a state A where U4 #4 0. Assume further that through 
some external influences we take the system from state O to state 
A. Then according to equation (4.2) the energy of state A is 


U, = —wet, (4.3) 


Apparently, equation (4.3) (and hence any definition of energy) 
depends on the standard state O. If instead of O we had chosen 
another standard state O’ we would have obtained U4. In this case 
(see Figure 4.3) we would have (since W*4 depends on initial and 
final states only) 


-W3 = -W31 — wad 
or 
U, = U} + Uo 
or 
U, — UL = Uo 


where Uo, is the energy of the standard state O’ and thus is a 
constant. It follows that had we chosen a different standard state, 
U'ʻ would differ by a constant (additive) amount. This constant is 
an essential feature of the concept of energy which has no effect on 
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the final results when we are working with differences (changes) of 
energy and not with actual amounts. 


4.3 Equivalence between heat and work done 


Consider a sample of water at a state i where the temperature is 
T; and at a state f where the temperature is T; > T;. In how many 
different ways can we take the system from 7 to f? Obviously we 
can simply warm up the system. In this case, and as long as Tọ 
is low enough for evaporation to be negligible, the volume of the 
system remains unchanged and therefore no work is done by the 
external forces. This is the first way. The second way is to raise 
the temperature via friction. Just imagine a set of paddles rotating 
inside the system. In this case, we will observe that as long as the 
paddles continue to rotate the temperature of the water increases. 
Since the water resists the motion of the paddles we must perform 
mechanical work in order to keep the paddles moving until T} is 
reached. 

Thus, the work performed by the system depends on whether we 
go from i to f by means of the first or the second way. If we 
assume that the principle of the conservation of energy holds, then 
we must admit that the energy transmitted to the water in the 
form of mechanical work from the rotating paddles is transmitted 
to the water in the first way in a non-mechanical form called heat. 
It follows that heat and (mechanical) work are equivalent and two 
different aspects of the same thing: energy. In the first way, where 
no work is performed, the change of energy can be expressed as 


AU= Qw-=o 


where Qw-=o is the amount of heat received. In the second way (the 
experiment can be thought of as taking place in some adiabatic 
enclosure) the change in energy can be expressed as 


AU = —w*. 


When both work and heat are allowed to be exchanged, we replace 
the above equations by the general expression 


AU+W=Q. (4.4) 


The quantity Qw-o required to raise 1g of water from 14.5°C to 
15.5 °C is by definition equal to 1 cal. In this case W*4 turns out to 
be equal to 4.185 J. This value is known as the mechanical equiv- 
alent of heat. Note that since dW is not an exact differential, 6Q 
is not either. Thus, for an infinitesimal process the above equation 
takes the form 

dU + ôW = ôQ. (4.5) 
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Note that for a cyclic transformation $ dU = 0 and thus W = Q. 
Equation (4.5) is the mathematical expression of the first law of 
thermodynamics. Since ÂW = pdV we can write the first law as: 


dU + pdV = 5Q (4.6) 


or for unit mass 





du + pda = ôq. (4.7) 


4.4 Thermal capacities 


Since dU is an exact differential we can write it as: 


w= (Fr) ar+ (F), dV (4.8) 
(Ge) ere(G) fe 8 
_ w= (22) ne(%) a any 


By combining equations (4.6) and (4.8) we get 


can) eer) |e ea 


By combining equations (4.6) and (4.9) and using the fact that 


OV 
dT 
w= (5), 0+ a), 
we arrive at 


(sr), +? (ar) e [(ap), E) 


(4.12) 
By combining (4.6) and (4.10) we get 


(= ) (5 ) 
dp 4 

Op) y OV r 

We define as thermal capacity, C, the ratio 6Q/dT and we denote 

it as Cy if the thermal capacity is measured at a constant V and as 


C, if it is measured at constant p. According to this definition if in 
a process at constant volume there is no change of the physical and 

















p| dV =6Q. (4.13) 
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chemical state of a homogeneous system, then the heat absorbed is 
proportional to the variation in temperature 


6Q = CydT 


or 
6Q = cymdT 

where m is the mass of the system and cy is the specific heat 

capacity at constant volume (i.e. cy = Cy/m). Another specific 


heat is the molar specific heat Cym = Cy /n where n is the number 
of moles. Similarly, at constant pressure 


ÔQ =C,dT 


or 


ÔQ = cpmdT. 


From equation (4.11) if we set dV = 0 we obtain 
ôQ OU 
c= a= (ar), 
and (4.14) 


_ (au 
Cy = ƏT “a 


From equation (4.12) if we set dp = 0 we obtain 


it ON see Ole 
2 NOE) AOE) 


By defining here a new state function, the enthalpy, H = U + pV 
(h = u + pa), it follows from the above equation that 


OH 
o= (Fr) 


and (4.15) 
_ (Oh 

=| zr) - 

p 
Again here we can define the molar specific heat at constant pres- 
sure AS Cpm = Cp/n. Also, since pV = nR*T, for ideal gases the 
enthalpy is a function of two state variables (U, T) and as such it 
can be expressed as an exact differential. Note that for isobaric pro- 


cesses dH = dU + pdV = 6Q, i.e. the change in enthalpy equals the 
amount of heat absorbed or given away. 
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Figure 4.4 
The experimental design 
of Joule’s proof that 


U =U(T). 
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4.5 More on the relation between U and T (Joule’s law) 


As we have stated earlier, any thermodynamical quantity can be 
expressed as a function of two state variables. As such we can 
assume that U = f(V,T). In Chapter 2 we discussed the fact that 
the kinetic theory of heat accepts that U is proportional to the 
absolute temperature. This indicates that U is a function of only 
one state variable, U = f(T). An experimental proof of this was 
provided by Joule, who performed the following experiment. He 
constructed a calorimeter and placed in it a container having two 
chambers, A and B, separated from each other by a stopcock (see 
Figure 4.4). He filled A with an ideal gas and evacuated B. After 
the whole system settled to a thermal equilibrium (as indicated 
by the thermometer in the calorimeter), Joule opened the stopcock 
thus allowing the ideal gas to flow from A to B until pressure every- 
where equalized. He observed that this did not cause any change in 
the temperature reading. From the first law we thus have Q = 0 or 
from equation (4.4) 


AU = —-W. 


Now note that since the volume of the system (the container) does 
not change, the system performs no work. As such W = 0 and 
AU = 0. But, you may ask, how about the gas inside chamber 
A?. Its volume increases as it expands to fill container B. What 
is the effect of that? Since there was no temperature change and 
no variation in energy during the process we must conclude that a 
variation in volume at a constant temperature produces no variation 
in energy. In other words U 4 f(V) which means that U is only a 
function of T: 


U =U(T). (4.16) 


Later in problem 6.2 you will be asked to prove this mathematically. 
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More on thermal capacities 


Because of equation (4.16), (0U/OT), = dU/dT. As such we can 
write equation (4.14) as 


dU 

Cy = ar 
(4.17) 

_ du 

Cy = aT 

Similarly, H =U + pV =U+nR*T = H(T) and 

dH 

= oF 
(4.18) 

_ dh 

C= op 


In the atmospheric range of temperatures, cy and c, are nearly 
constant. From equation (4.17) we then have that 


U= fevers constant 


or 
U = CyT + constant. 


Since only differences in internal energy (and thus in enthalpy) 
are physically relevant, we can set the integration constant to zero 
and obtain 


U x CyT. 


Similarly H ~ C,T,u ~ cyT, and h ~ cT. Equations (4.17) and 
(4.18) can be combined to obtain 


C, — Cy =nR* (4.19) 
which implies that 
Cp— cy = R, Com — Cvm = R*. (4.20) 


Recall now from Chapter 2 (equation (2.8)) that for a monatomic 
gas consisting of N point-masses the total internal energy is 











3 
= - NKT. 
4 2 
It follows that 
dU 3 3 3 Srah 
Cy IT zE 5k Fe ie gn ; 
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and 3 3 3 
* aoe n Peco woe 
Cum = git 3 Cy 2 ae 5k 
Similarly, 
dH dU ge D ey 
oa apt pe 
and 5 5 
Cpm = af Cp = ait 
For a diatomic gas these values change to 
5 5 5 
Cy = ze, Cvm = Jea Cy = ait 
and 
7 7 6 


Cp = gk", Cpm = gf: Cp = at 


Note that often in the literature the values of Cy and Cym or Cp 
and Cpm are assumed the same. This is correct only in the case of 
one mole. From the above it follows that 

Cy re, Cpm Cp 


Cy Cvm Cy 





We will denote these ratios as y. 


Dry air is considered a diatomic gas. According to the above 
relationships we estimate that for dry air 


ya = 1.4 

Cva = 718 Jkg K7! = 171 calkg™* K7! 

Cpa = 1005 Jkg™* K~* = 240 cal kg™* K7* 

Ra = 287 J kg K7. 
According to the above definitions we can write the first law as: 

CydT + pdV = ôQ (4.21) 
or 
cydT + pda = ôq 


Considering that Cy = C, — nR* and that pV = nR*T, equation 
(4.21) becomes 


(Cp — nR*)dT + (nR*dT — V dp) = 6Q 


or 


C,dT — Vdp = 6Q. (4.22) 


For convenience Table 4.1 summarizes all the different expressions 
of the first law. 
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Table 4.1. Expressions of the first law 








For a gas of mass m For unit mass 





dU + ôW =6Q du + dw = ôq 
dU + pdV = 6Q du + pda = ôq 
CydT + pdV =Q cvdT + pda = ôq 
C dT —Vdp=6Q cpdT — adp = ôq 














4.6 Consequences of the first law 


Let us look at the form of the first law in the following special cases. 


. me a= tant 
e Isothermal transformations: i 20 f 


During such transformations dT = 0 and as such from equation 
(4.21) it follows that 6Q = dW. Then the amount of heat exchanged 


f f 
q= | sw = | pdV 


1S 
f aV 
V 


or 


Q=nR*T 
or 
V; 
=nR*T In. 
Q=n ny 


eè Isochoric transformations: i © "3°" f 

In this case dV = 0 and from the first law it follows that 
6Q = dU = CydT 

or 


f 


e Isobaric transformations: i 25 f 
In this case it follows that 
6Q => C,dT or Q = C,(T; E T;), 


ôW =pdV or W = p(V; — V;). 
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and 
AU = C,(T; — T;) — p(V; — Vi) 


e Cyclic transformations i — f — i 


As we know, in cyclic transformations ¢ dU = 0. Then, from the 


first law it follows that 
gow z f 5Q. 


e Adiabatic transformations — Poisson’s relations 


During an adiabatic process there is no exchange of heat between 
a system and its surroundings by virtue of temperature differences 
between them. Thus, 6Q = 0. Then the first law can be expressed 
as 








dU = —ôW 
or 
CydT = —pdV 
or 
C dT = Vdp 
or 
dT V dp 
T= Ot 
or using pV = mRT 
dT  mRdp 
T ©p 
or 
dT nk* dp 
T Cp p 


or using C, — Cy = nR* 


dT ( ce) dp 








T Cp) P 
or 
dT 1 \ dp 
= {1 i 4.23 
T ( 3 p n 


Integrating equation (4.23) yields 


1 
ln T= (1 — ) In p + In(constant) 
Y 


or 
ln T = ln p7 + In(constant) 


Figure 4.5 

The various relationships 
between adiabats, 
isotherms, isobars, and 
isochores in partial and in 


the complete state spaces. 
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or 


y-1 
T = constant Xx p 7 
or 
1-7 
Tp 7 = constant. (4.24) 


Through the ideal gas law we can then derive the following 
equivalent expressions of equation (4.24): 


TV = constant (4.25) 


pV” = constant. (4.26) 


Equations (4.24)—(4.26) express Poisson’s relations for adiabatic 
processes. Since for dry air y = 1.4 > 0 equation (4.25) proves 
something that we all more or less know: that as an air parcel 
rises and expands (i.e. V increases), its temperature must decrease. 
For adiabatic expansion, since dQ is zero this decrease is due to 
expansion work that the parcel performs on the environment. 
From the ideal gas law it follows that isotherms are given by the 
equation 
pV = constant. 


Thus on a (p,V) diagram isotherms are equilateral hyperbolas 
(Figure 4.5(a); an equilateral hyperbola is one whose asymptotes 
are perpendicular to each other). From equation (4.26) we see that 
the adiabats on a (p, V) diagram will also be equilateral hyper- 
bolas but since y > 1 the adiabats will be steeper than the 


Ty 
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isotherms. Obviously in a (p,V) diagram isobaric and isochoric 
transformations are represented by straight lines on the p and V 
axis, respectively. On a (p, T) diagram the isochores are given by 
the equation 
nR*T 
Pp 


V= = constant’ 





or 
Tp * = constant 


which is the equation of a straight line. From equation (4.24) it 
also follows that in a (p, T) diagram the adiabats are again hyper- 
bolas (Figure 4.5(b)). On a (7,V) diagram (Figure 4.5(c)) the 
isobars are straight lines (as before TV~* = constant) and the adi- 
abats are hyperbolas (see equation (4.25)). The complete (p, V, T) 
space and transformations are shown in Figure 4.5(d). Since they 
require no heat exchange between a parcel and its environment (i.e. 
they require perfect isolation), adiabatic processes are idealizations 
that never happen in nature. However, for processes that happen 
quickly, energy transfer is not very significant. Such processes are 
thus adiabatic to a good approximation. 


e Polytropic transformations 
During an isothermal transformation the system can exchange 
heat with its environment but its temperature remains constant. 
This, however, is an ideal situation. In reality the equiliza- 
tion of temperatures does not happen instantaneously. Simi- 
larly, an adiabatic transformation is an ideal situation since it 
requires perfect heat insulation. The same comments can be 
made for isobaric and isochoric transformations. Thus, the rela- 
tions pV = constant, pV’ = constant, Tp™t = constant, and 
TV-t = constant are special cases of the general relationship 
pV” = constant where 7 can take on any value. Such transforma- 
tions are called polytropic transformations (in Greek polytropic 
means multi-behavioral). If we assume that 7 = y/e where e 
is a constant, we can work our way backwards from equation 


i= 
pV" = constant or Tp ™ = constant to arrive at 


7-1 =) 
T y) p` 


Considering that y = & the above equation can be expressed 


as 
dT  (C,—Cye\ dp 
ee Cy p 








or 





dT rs Cy — Cy — Cye+ Cy dp 
P C. p 


P 
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or 
dr (25) dp  Cvy(1-— €) dp 
T Cy p` Cp p 
or 
dT _ nk" dp | (1-6) dp 
T œp Y pP 
or 
dr V dp W= e) dp 
TO Cp T y P 
i 1 CLT 
C,dT = Vdp + Ty 
i Gezi yv 
T= Je 
C,dT = Vdp pe 


Then, according to the first law, the second term on the right- 
hand side is the heat exchanged, which is zero if e = 1 (i.e. 7 = y, 
adiabatic process), greater than zero (i.e. the system absorbs heat 
from the environment) if 7 < y (e > 1) and dp < 0, orn > 
(e < 1) and dp > 0, and less than zero (i.e. the system gives 
away heat to the environment) if n > y and dp < 0, or 7 < y and 
dp > 0. From the above equation we also find that for € = 0 (i.e. 
n = œ) nR*dT = Vdp, which owing to the ideal gas law implies 
dV =0. For e = y (ie. 7 = 1), it reduces to C,dT = 0, while for 
€ = œ (ie. 7 = 0) it dictates that dp = 0. As such a polytropic 
process pV” = constant reduces to an isobaric, to an isothermal, 
to an adiabatic, and to isochoric process for 7 = 0,1,y, and oo, 
respectively. 

Having said this we must add that in the atmosphere over a 
rather wide range of motions the timescale for an air parcel (our 
system) needed to adjust to environmental changes of pressure, 
and to perform work, is short compared with the corresponding 
timescale of heat transfer. For example, above the boundary layer 
and outside the clouds the timescale of heat transfer is about two 
weeks whereas the timescale of displacements that affect a parcel 
is of the order of hours to a day. It thus appears that adiabatic 
approximations are good approximations for many atmospheric 
phenomena. 


e Dry adiabatic lapse rate 
From equation (4.24) we have 


ee) 
T = constant -p 7 . 


By taking the logarithmic derivative of the above equation we 
have 


Ai hiy 
y 
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or 


or 
1dT y-—11dp 


Tdz y pdz 





(4.27) 


For large-scale motions the hydrostatic approximation states 
that pressure gradient force balances the force due to gravity. 
Therefore, 


dp 
q = PS: (4.28) 


For a parcel of air rising in the atmosphere, equation (4.27) is 
valid as long as its ascent is an adiabatic process. For the same 
parcel equation (4.28) is valid only if dp/dz experienced by the 
parcel is equal to that of the large scale (surroundings). We 
assume this is more or less true but in equation (4.28) p rep- 
resents the density of the surroundings which from the ideal gas 


law is 
p 


PEA 4.29 
RT; (4.29) 


p 


where Ts denotes the temperature of the surroundings. In this 
case by combining equations (4.27)—(4.29) we have that 

dr!  y-igT 

dz y RT; 





Recall that by definition an adiabatic ascent or descent is one in 
which there is no exchange of energy with its surroundings by 
virtue of temperature differences between them. Accordingly, if 
the parcel remains dry (i.e. no condensation takes place) we may 
define from the above equation the dry adiabatic lapse rate, Ij, as 
the atmospheric temperature profile such that the temperature 
of the parcel is always at the temperature of its surroundings. 
For such a profile T = Ts and as such for dry air 

af y-1g g 


y= = 


= = 9.8°Ckm™?. 4.30 
dz y Ra Cpa ai ( ) 





This value is larger than the observed average decrease of envi- 
ronmental temperature with altitude. The difference is mainly 
due to the fact that in the derivation of [4 we neglected the 
presence of water vapor. When vapor is present we have new 
lapse rates, which we will discuss in a later chapter. 
Recalling that h = u + pa and that h = cpaT we can write 
equation (4.30) as 

dh 

dz I 
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or 


Aei +gz)=0 (4.31) 
dz 

where h + gz is defined as the dry static energy. The above 
equations indicate (1) that enthalpy of a parcel decreases as it 
rises adiabatically because it is doing gravitational work, and (2) 
that the static energy (the sum of enthalpy and the gravitational 
potential energy) is conserved in an adiabatic motion (ascent or 
descent). Enthalpy is thus the specific internal energy plus the 
term pa that accounts for the work done by the parcel on the 
surroundings. 


e Potential temperature 
Let us assume that there is a variable 0 defined by the relation 
6 = ATp-* where A and 8 are constants. By taking the loga- 
rithm on both sides we get In 0 = In A+ 1n T — £ ln p. Then by 
differentiating we arrive at 


dln 0 =dln T — d ln p. 


Consider now the first law in the form c,dT —adp = ôq and divide 
it by T. Then with the help of the ideal gas we obtain 
tare nee 
Cc Cyl’ 


P P 


For 6 = R/c, the above two equations combine to give 


ôq 
In 0 = —. . 
dln a (4.32) 


p 


For adiabatic processes ôq = 0 and thus d ln 0 = 0. As such we 
should expect that in adiabatic processes there should exist a 
measure 0 which is conserved. This measure can be defined as 
follows. 

From equation (4.24) we have 


where k = (y—1)/y = 1—cy /cp = R/c, = 0.286 (for dry air). The 
state (To, po) can be taken as a reference state. As such we can 
choose po = 1000 mb and denote the corresponding temperature 
as Ty. The above equation is expressed as 


k 
n=r() 
p 
or (4.33) 
To = poTp™. 
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If we compare the above expression to 6 = AT'p~* we see that 
To = 6 for A = pë and 8 = k. We call 0 the potential tempera- 
ture and we regard it as the temperature a parcel will have if it is 
compressed or expanded adiabatically from any state (T, p) to the 
1000 mb level. It follows that the potential temperature remains 
invariant during an adiabatic process. Thus, under adiabatic con- 
ditions, 0 can be used as a tracer of air motion. On timescales for 
which a parcel can be considered as adiabatic, constant values of 
0 track the motion of this parcel, which takes place on a surface 
of constant 0. From equation (4.33) it is clear that the distri- 
bution of 0 in the atmosphere depends on the distribution of T 
and p. In the atmosphere dp/dz >> dT /dz (just think that 5km 
above the surface pressure decreases on the average from about 
1000 mb to 500 mb while the temperature decreases from 288 K to 
238K). Therefore, surfaces of constant 6 tend to be like isobaric 
surfaces. Before, we saw that for an ascent (descent) the parcel’s 
temperature must decrease (increase) owing to the work done by 
the parcel (by the surroundings). With the definition of potential 
temperature we can extend this statement to read as follows: dur- 
ing an adiabatic ascent (descent) the parcel’s temperature must 
decrease (increase) but in such a proportion as to preserve the 
parcel’s potential temperature. Note that equation (4.32) indi- 
cates that for non-adiabatic processes the change in the potential 
temperature is a direct measure of the heat exchanged between 
the parcel and its environment. As a result a non-adiabatic par- 
cel will drift across potential temperature surfaces in proportion 
to the net amount of heat exchanged with its environment. 


Examples 


(4.1) One mole of a gas expands from a volume of 10 liters and a 
temperature of 300 K to (a) a volume of 14 liters and a temper- 
ature of 300K and (b) a volume of 14 liters and a temperature 
of 290 K. What is the work done by the gas on the environment 
in each case? 


(a) According to the definition of work 


V2 V2 R*T 
w=] pav=] ” 


Via Vi V 


dv. 





Since the temperature remains the same the above 
equation gives 
V; 
W =nR*T In z = 839 J. 


1 


Figure 4.6 
Diagram for example 4.1. 
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(b) In this case T does not remain constant. As such it 
cannot be taken outside the integral and the above 
equation cannot be applied. 


When T does not remain constant, in order to calculate the 
work done we must have the explicit function describing the 
path from i to f in the (p, V) diagram, not just i and f. An 
approximation to this would be to define a temperature T 
or a pressure Pp that satisfied the relationship 


L TdaV T dV 
v V v V 


or the relationship 


Vo Vo 
I pdV =p 1 dV 
Vi Vi 


Since in our case no clues are given about the function 
describing the path from i to f we are free to make (rea- 
sonable) assumptions. We can thus assume that i and f 
are connected by a straight line. From the data we can 
estimate that p; = 249420Pa and py = 172200Pa. From 
Figure 4.6 it is easy to see that the work done is the area of 
the trapezoid ifV;V;. Because of the equality of the shaded 
triangles this area is equal to the area of the rectangle 
defined by the lines p= 0, p = P, V = Vi, V = V; where 
P = pı + p2/2. Thus 








W = P(V; — V;) = 843 J. 


Pi 
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Figure 4.7 pt 
Diagram for example 4.2. 


Pi} $% , 


~< T=constant 


V; V; V 


(4.2) Consider an initial, 7, and a final, f, state of a gas on an isotherm 
with V; > V;. The gas goes from į to f via an isobaric expansion 
and subsequent isochoric cooling; sketch on a (p, V) diagram 
the complete transformation and find the work done. If instead 
the gas goes from 7 to f via an isochoric cooling followed by an 
isobaric expansion, will the work done be different from before? 


In the first case the gas goes from i to m (isobaric from V; to 
V;) and then from m to f (isochoric from p; to ps) (Figure 
4.7). The total work done from i + f will be the sum of the 
work done from i —> m and from m > f 


Vf Vf 
w= | pav + f pdV 


Vi Vy 


Vf 
=p. | dV = p;i(V; — V;). 


y 


i 


In the second case the gas goes from i to m’ (isochoric from 
p; to pp) and then from m’ to f (isobaric from V; to Vp). As 
before we have that the total work from i —> f is 


Vi Vf 
W: = f pdV + / pdV 


Vi Vi 


Vy 
= v f dV = p;(V; — Vi). 


V; 


i 


Since p; # py it follows that W, # W3. 
(4.3) Calculate Q,W, and AU for the following processes 
(Figure 4.8): 


(a) isothermal reversible compression from state i = (p;, V;) 
to state f = (py, Vy); 
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Figure 4.8 p 


Diagram describing the 
processes in example 4.3. 


PF 


Pi 








(b) 


(c) 


f m m 
x Isochore 
Adiabat 
Isotherm T 
i 
> 
vV; v, v, V 


adiabatic reversible compression from state i = (p;, V;) to 
state m = (pf, Vm) and a subsequent isobaric reversible 
compression to state f = (py, Vy); 

reversible isochoric increase of the temperature from 
state i = (pi, V;) to state m’ = (p, V;) and a subse- 
quent reversible isobaric decrease in temperature to state 


f = (ps, Vp). 


Express all answers in terms of p;, pf, and T, where T is the 
temperature of the isotherm in (a). 


(a) 


During this transformation T; = T; = T, AT = 0 and 
AU = 0. Then from the first law it follows that 


f fd 
q=w- | pdv =nkT | d 


= nR*T ln t; =nR*T ln Pi. 
V; Pr 


This transformation consists of two branches (i > m 
and m — f). From (a), the initial and final temper- 
ature is the same, i.e. AU = 0. Thus, again we have 
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Q=W with 
AU = AU; pm Res, 
Q z= Qism AR Quis} 
and 
W = Wim a Wat: 





Since the branch i — m is adiabatic it follows that 
ism = 0 and 


Tm and T, are related via the equation 


1-7 1-7 


Tipi” = Tope , 


where Pm = p;. It follows that 


i-v dey, 
ron(5)”=1(5) 
Pr Pf 


Then from equation (4.34) we obtain 


(=) ro J ; (4.35) 


Branch m — f is isobaric. Thus, 


Wim = —-CyT 





f 
Wor = J pdV = pi(Vs — Vm). (4.36) 


m 


Vm and V; are related via Poisson’s equation as follows 
pi V? = Pm Va 


Since pm = ps the above equation becomes 


(FE) =F 
V; Py 


N Vy 
Vin = Vi (=) l 
Pf 


Then equation (4.36) can be written as 


A Vy 
nfe" 
Pf 
; ra 
Vi \ Pf 





or 





= pV} 
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ioe 
i= (=) _ (4.37) 


From equations (4.35) and (4.37) it follows that overall 
any 
ie 
Pf 
iai 
eto 
Pf 
iar 
i= (=) . 
Pf 


(c) Here also the transformation consists of two branches, 
i => m and m’ > f. Branch i —> m’ is isochoric, so 


= nRT 





=nR*T 





Q=W=-Cy\T 





+ nR*T 





=C,T 





Wim! = 0. 


Branch m’ —> f is isobaric, which gives 


f 
Wms = I pdV = pi(V; — Vm) = pi Vs — Vi) 


or 


V; 
Wms = ps Vr (L- — 
aes 


= nR'T (1 = 2r) 
Pi 


Here again T; = Ty = T, AT = 0 and AU = 0. 
Thus, overall 


g=w-=nrr (1-2), 


4 


(4.4) This following problem is due to Iribarne and Godson (1973). 
Students should try to solve it before looking at the solution 
here. Figure 4.9 shows an insulated box with two compartments 
A and B. Both compartments contain a monatomic ideal gas and 
they are separated by a wall that on one hand does not allow 
any heat through it but on the other hand is flexible enough to 
ensure equalization of pressure in both compartments. The initial 
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Figure 4.9 
Diagram for example 4.4. 


LA nee ee 


conditions for both compartments are the same: T; = 273K, 
V; = 1000 cm?, and p; = | atmosphere (1013 x 10?Pa). Then, 
by means of an electrical resistance, heat is supplied to the gas 
in A until its pressure becomes ten times its initial pressure. 
Estimate: (a) the final temperature of the gas in B, (b) the work 
performed on the gas in B, (c) the final temperature of the gas 
in A, and (d) the heat absorbed by the gas in A. 


(a) 


The heat that is supplied to gas in A results in an 
increase in all its state variables. As it expands, it com- 
presses the gas in B but it does not provide it with any 
of its heat because of the insulating wall between them. 
Since the wall is flexible the gas in B is compressed 
adiabatically until its pressure equals that inside com- 
partment A. As such pgs = 10 atmospheres. From 
Poisson’s equation we can then solve for Tgș 


= 1-7 


i-v Izy 
TpiPpi =TBf PBF 


Since the gas is monatomic y = 1.666. Thus, 





—0.4 
Taj = (=) Tp: ~ 686K. 
PBf 


Having Tgş we can use the other form of Poisson’s 
equation to find Vz: 


—1 —1 
Tp: Vg; = Tay Vay 


Tiy 273\ ee 
Ver = (=) Vai = (as) Vai 


= 0.25 Vg; = 250em*. 


or 
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(b) The work performed on the gas in B is given by 
W = —AU — —Cy (Tey T Tp) = —413 Cy is 


As we know, for a monatomic gas Cy = ae The 
number of moles n can be found from the ideal gas law 
and the initial conditions. It follows that 


_ DVi 


= eq, = 020446 mol. 


n 





Accordingly, Cy = 0.556 J K~* and 
W x —230J 


or 
W x —55 cal. 


(c) Since the final volume in B is 250cm? it follows that 
the final volume in A is 1750 cm?. By applying the ideal 
gas law for the initial and final conditions in A we find 
that 

V, 
Tay = AAI i = 4777 K. 
PaiVai 

(d) The heat absorbed by the gas in A is given from the 

first law 


Q=AU+W. 


Here W = 55 cal (i.e. the opposite of the work done on 
gas B), and AU = Cy(Ta; — Tai) = 25043 = 598 cal. 
Thus 

Q = 598 + 55 = 653 cal. 


Problems 


(4.1) 


(4.2) 


Commercial aircraft fly near 200mb where typically the 
outside temperature is —60°C. (a) Calculate the tempera- 
ture of air if compressed adiabatically to the cabin pressure 
of 1000mb. (b) How much specific heat must be added 
or removed (isobarically) to maintain the cabin at 25°C? 
Consider the air as dry air. (337.5K,9.5calg7! has to be 
removed) 

A sample of 100g of dry air has an initial temperature of 
270K and pressure 900 mb. During an isobaric process heat 
is added and the volume expands by 20% of its initial vol- 
ume. Estimate: (a) the final temperature of the air,(b) the 
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(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10 


m= 


amount of heat added, and (c) the work done against the 
environment. (324 K, 5427 J, 1550 J) 

An air parcel moves from p; to pr. If its initial temperature 
is T; find: (a) the specific work done by or on the parcel 
if its change happens adiabatically, (b) the specific work if 
its change happens isothermally, and (c) the change in the 
parcel’s potential temperature in (a) and (b). 

One mole of dry air has an initial state T = 273K and 
p = 1 atmosphere. It undergoes a process in which its 
volume becomes four times its initial volume at 400 mb. If 
air is considered an ideal gas and if the process obeys the 
law pV” = constant, estimate (a) the value of 7, (b) the 
final temperature, (c) the change in internal energy, the work 
done, and the heat exchanged between the air and its sur- 
roundings. Look at your results carefully and elaborate on 
what kind of process would generate results like these. Is 
this process a realistic process? (0.67, 431K, 3286 J, 3989 J, 
7275 J) 

Assume that you are on the top of a mountain at an altitude 
of 15000 feet and there are no clouds above or below. If the 
temperature is —12°C what will the temperature be at 3500 
feet altitude? (22.4°C) 

A parcel of dry air has a volume 10 liters, a temperature 
27°C, and a pressure 1 atmosphere. The parcel (a) is com- 
pressed isothermally until its volume becomes 2 liters, and 
(b) expands adiabatically until its volume becomes 10 liters. 
Describe these changes graphically in a (p,V) diagram. On 
it mark the values of p, V,T for each case. 

A sample of dry air has a temperature of 300 K, a volume 
of 3 liters, and a pressure of 4 atmospheres. The air then 
undergoes the following changes: (a) it is warmed under a 
constant pressure to 500K, (b) it is cooled under constant 
volume to 250K, (c) it is cooled under constant pressure 
to 150K, and (d) it is warmed under constant volume to 
300K. (1) Describe graphically on a (p, V) diagram each of 
the changes and mark on it the final values for pressure and 
volume after each change. (2) Calculate the total work done. 
(405.2 J) 

For an ideal gas what on a (p, V) diagram describes the lines 
for which U = constant? 

An ideal gas undergoes two transformations from state 1 to 
state 2 which follow the two curves in Figure 4.10. In which 
of the transformations is (a) the change in internal energy 
greater? (b) the heat absorbed greater? 

One kilogram of dry air is warmed from 20°C to 70°C under 
a constant pressure of 1 atmosphere. Calculate: (a) the heat 


Figure 4.10 
Diagram for Problem 4.9. 
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(4.11 


(4.12 


(4.13 


(4.14 


(4.15 


WS 


Ne 


) 
) 


— 


T= constant 








absorbed by the gas, (b) the work done, and (c) the change 
in internal energy. (50250 J, 14350 J, 35 900 J) 

Hydrogen under p = 5 atmospheres and T = 20°C is 
warmed isobarically until its volume increases from 1 to 2 
liters. Calculate the absorbed heat. (1773.4 J) 

Inside a cylinder closed by a piston there is 88g of an 
unknown diatomic gas of a temperature 0 °C. The gas is then 
compressed adiabatically until its volume becomes equal to 
1/10 of the initial volume. The change in internal energy is 
17 158 J. Identify the gas. (Carbon dioxide) 

Provide an example that proves that heat transfer is an 
irreversible process. 

At T =0°C and p = 1000 mb, 1g of dry air receives during 
an isochoric process an amount of heat. It is then observed 
that its pressure increases by 50 mb. What is the change in 
the temperature of the air and what is the amount of heat 
absorbed? (13.6°C, 2.33 cal) 

Calculate the change in specific energy of an adiabatic 
parcel whose speed changes from 10ms~! to 25ms7}. 
(—262.5 Jkg~*) 


CHAPTER FIVE 





The second law of 
thermodynamics 


The first law of thermodynamics arose from the conservation of 
energy principle. The first law, even though it implies that we can- 
not create or destroy energy, places no limits on how energy can be 
transformed from one form to another. Thus, on the basis of the 
first law, heat can be transformed into work, work into heat, work 
can be done at the expense of internal energy, and so on. However, 
if no other laws existed the first law would allow certain phenom- 
ena to happen that never happen in reality. For example, consider a 
heavy body falling on the ground. We will observe that during the 
impact the body will warm. The opposite phenomenon according 
to which a body at rest on the ground begins to rise by itself while 
it is cooling is impossible. Similarly, no engine has yet been built 
which, for instance, would receive heat from the sea, transform it 
to work, and then set a ship in motion. Both the above examples 
are not in disagreement with the first law since the work would be 
done at the expense of the internal energy of the soil or the sea. 
The impossibility of these phenomena is due to the second law of 
thermodynamics, often hailed as the supreme law of nature. We 
will start our discussion of this law with the following example. 


5.1 The Carnot cycle 


The Carnot cycle is a thermal engine. A thermal engine is one 
that receives from some source an amount of heat, part of which it 
transforms into work. In thermodynamical terms the Carnot cycle 
is a cyclic transformation consisting of the following four steps 
(Figure 5.1): 


Step 1: reversible isothermal expansion (1 — 2), T=T, = 
constant 
Step 2: reversible adiabatic expansion (2 — 3), 6 = 6, = constant 
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Figure 5.1 pa 
Steps in a Carnot cycle. 











Vi V4 V> V3 V 


Step 3: reversible isothermal compression (3 — 4), T = T, = 
constant 

Step 4: reversible adiabatic compression (4 —> 1), 0 = 0, = 
constant 


with T, < Ti, 02 < 04. 

Let us now calculate the work done and heat absorbed when an 
ideal gas is subjected to this cycle. The change in internal energy 
is zero because the overall transformation is a cyclic one. We will 
consider each step separately. 


Step | 
AT =0 
a “%2 aV Vo 
W. =] pdV = nR*T, I — = nRT, ln = 
12 P 1 ie V 1 Vi 
AU 12 = Cy AT = 0 
Qı2 = Wie 


Since V, > V, it follows that W,2,Qi2 > 0 indicating that during 
step 1 work was done by the gas and heat was absorbed by the gas 
from the heat source at temperature 7}. 


Step 2 
Qos = 0 
AU 23 = Cy (To — Tı) < 0 
Was =] —AU 43 = —Cy (Te == Tı) > 0 
Step 3 


AT =0 
Va Vi 
Wa == / pdV = nR*Tə In Vv. 


V3 3 
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AU 34 = 0 
Q34 = Wa 
Since V, < V3 it follows that W34, Q34 < 0 indicating that now work 


is done on the gas and heat is given away by the gas to the source 
at temperature Tz. 


Step 4 
Qu =0 
AU 4, = Cy(T; = T) = 0 
War = —AU 4, = -Cy (Tı baa Tə) < 0 
From the above it follows that the total work done is 
W = Wiz + Wo3 + W4 + Way 
V; V, 
= nR*T, In v +nR*T, In A (5.1) 
Because transformations 2 — 3 and 4 — 1 are adiabatic the 
following relations apply: 
HVIS = DV 
and 
nV a Ta V~. 
By dividing the above two equations we have 


Va _ Vs 


Yo Vi 
Accordingly, equation (5.1) can be written as 


V; V. 
W =nR*T, ln A — n R*T, ln v 
or 


W =nR* (1m 2) (T, — Tp). 


1 


The total amount of heat absorbed during this cycle is 


Q = Qiz et Qsa = Wia + W34 (5.2) 
or 
V2 
Q=W=nR* (1m =) (T-T) (5.3) 


For simplicity we will now denote the amount of heat absorbed by 
the gas from the heat source at T, (i.e. Q12) as Qı and the amount 
of heat given by the gas to the source at T, (i.e. Q34) as Qo. Since 
W is the area enclosed by the cycle it follows that W > 0. Then 
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from equation (5.3) it follows that Q > 0. Since W = Qı + Q2 and 
Qs is negative it follows that only part of the heat absorbed by 
the gas at the source T, (higher temperature) can be transformed 
into work. The other part is surrendered to the source at To (lower 
temperature). We define the efficiency of the Carnot cycle, 7, as 
the ratio between the work done and the heat absorbed by the gas 
at the source T}. 
n= PEE (Q > 0,02 <0), 

By considering the relations derived for Q, and Q> we can write 
the above ratio as 


nR*T, in 
ar a 
nR Tı h 
or 
To 
=1- 5.4 
n T, (5.4) 


5.2 Lessons learned from the Carnot cycle 


e According to equation (5.4) the thermodynamic efficiency of the 
cycle depends only on the temperatures of the two heat sources 
and it becomes zero if Tz = T, (i.e. if the two sources collapse into 
one source). As such we conclude that doing work via a thermal 
engine operating at only one source of heat is impossible. This is 
known as Kelvin’s postulate and is an expression of the second 
law of thermodynamics. Another way of stating this postulate 
is that it is impossible to construct an engine which transforms 
heat into work without surrendering some heat to a source at a 
lower temperature. This explains why stones on the ground do 
not rise in the air or why no machine has been discovered which 
by taking heat from the sea or the air sets in motion a ship or 
a car. According to the second law, apart from the warm heat 
source (sea or air), a second colder heat source would be required 
which we do not really have. I use the word really because some 
differences in temperature do exist in both the oceans and the 
atmosphere. For example, the temperature of the surface air is 
around 35K warmer than the air at 5km, and the temperature 
at the surface of the oceans is about 20 K warmer than at a depth 
of 1km. In both cases, however, it turns out that the Carnot effi- 
ciency is about 10%. This is a low efficiency given the tremendous 
expenses that would be required to build such engines. On the 
other hand if the atmosphere, in the horizontal, is thought of as 
a Carnot cycle operating at T} = equatorial temperatures and 
Tə = polar temperatures we can explain why in the winter the 
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circulation is stronger. Recall from above that 





T: 
(2 ape 
Qı Tı 
or 
QtQ h-r 
Qı To 


If we think of summer as the limiting case where T, — T, > 0, 
then during summer Q2 + —Q, and 7 —> 0. On the other hand if 
we think of winter as the limiting case where Ta << Tj, then 
during winter 7 —> 1 and Q2 —> 0. It follows that the eff- 
ciency of this “machine” is higher in the winter than in the 
summer. Accordingly, the amount of heat given by the atmo- 
sphere to the cold source is smaller in the winter, than in the 
summer. Since the amount of heat absorbed by the atmosphere 
from the warm source is more or less constant, this means that 
there is more heat “remaining” in the atmosphere in the win- 
ter than in the summer that can be transformed into kinetic 
energy. 

Apart from the impossibility of transformations whose only final 
result is to transform into work heat extracted from only one 
source of heat, there is another type of transformation that is 
impossible. Such transformations refer to transfer of heat from 
a cold body to a warm body without work being done. This 
phenomenon is allowed by the first law (since heat lost by the 
cold body will be exactly equal to the heat gained by the warm 
body) but has never been observed. Of course transformations 
that result in heat flow from a colder body to a warmer body 
are allowed if work is done (such as in refrigerators). This leads 
to another expression of the second law of thermodynamics. A 
transformation whose only final result is to transfer heat from a 
body at a given temperature to a body at a higher temperature 
is impossible. This is known as Clausius’s postulate. It is easy 
to show that both postulates are equivalent and two expressions 
of the same law. To do that we shall prove that if Clausius was 
wrong then Kelvin would be wrong and vice versa. If Kelvin was 
wrong and we could transform into work heat extracted from only 
one heat source, then we could transform this work (by means 
of friction) into heat and with this heat raise the temperature 
of a given body of a higher initial temperature. This would be 
a violation of Clausius’s postulate. Now if Clausius was wrong 
and we could transfer a certain amount Q, from a source at 
T, to a source at Tı (Tə < T,) in such a way that no other 
change in the state of the system occurred, then with the help 
of a Carnot cycle we could absorb this amount and produce an 
amount of work, W. Since the source T, receives and gives up 
the same amount of heat (Qı) it suffers no final change. But 
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what we have just described is a process whose only final result 
is to transform into work the heat extracted from a source at the 
same temperature T, throughout. This is contrary to Kelvin’s 
postulate. 


e From the definition of efficiency it follows that 


E Q2 To 
Tr = & 
T° Oy (Qı > 0,Q2 <0). 


Based on the above relationship an important application of 
Carnot cycle is that it offers a definition of the absolute tempera- 
ture (Kelvin) scale based on pure thermodynamic arguments. For 
this we only need to have a value for a given temperature. An 
example for such value could be the melting point of ice. As such 
if we choose for the melting point of ice the value T} = 273°C 
every other temperature Tə can be defined via a Carnot cycle 
operating between heat source T, (consisting of melting ice) and 
a heat source T>. If we then measure Q, and Q> we can apply 
the above relationship to calculate T,. This thermodynamic def- 
inition of the Kelvin scale allows us to present a thermodynamic 
definition of absolute zero: the absolute zero is the temperature 
of a heat source (assumed to be the cold source of a Carnot 
cycle) which absorbs no heat (i.e. Q2 = 0) even though the cycle 
produces work. 


e In reality natural processes are irreversible processes. However, 


if some processes happen very slowly and the heat losses are 
small enough not to affect the temperature (which supposes that 
friction is also negligible), then the results will not differ from 
those expected from perfectly reversible processes. Such processes 
are often called partially reversible. If on the other hand fric- 
tion and heat flow are important then the degree of reversibility 
decreases. In nature the degree of reversibility varies between 
almost perfectly reversible and perfectly irreversible. 

If in a thermal engine we increase the degree of irreversibility 
(for example, by increasing friction or direct leakages of heat), 
then the work from a given amount of heat will decrease. If 
the thermal losses keep on increasing we will end up with a 
perfectly irreversible engine. It follows that the efficiency of a 
thermal engine can take on values from zero (corresponding to 
the worst thermal engine operating at perfectly irreversible condi- 
tions) to some maximum value (corresponding to the best engine 
operating at perfectly reversible conditions). The Carnot cycle 
is a perfectly reversible cycle and as such it gives the maximum 
value the efficiency can take. This expresses the Carnot theorem 
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which states that it is impossible to construct an engine oper- 
ating between two heat sources which would have an efficiency 
greater than the efficiency of a Carnot cycle operating between 
the same heat sources. However, this maximum value will not 
be equal to one as this would require that T> = OK which is 
practically impossible. Recall from section 5.1 that 


7 
Qi = nR*T, In v 


and Vv 
Q2 = nR*T, In r 


3 
It follows that 


Qi V2 
nR* In 
Pe V, 
and Q 7 
2 2 
— = —nR* ln = 
2 Vi 
or Qi Q 
1 AA . 
T, + 17 0 (Qi >0,Q2 <0) (5.5) 


Now by considering equation (5.4) we can, for the reversible 
Carnot cycle, rewrite equation (5.5) as 


Q1Nrev 4 Q2 
T-b h 


Since Q2/T> is less than zero it follows that QıMev/(Tı — Tə) 
must be greater than zero. Since their sum must equal zero 
and, as we discussed above, the value of 7 is maximum for per- 
fectly reversible processes, it follows that for irreversible processes 
where 7 becomes smaller 

Qı Nirrev Qe 


S2 <o. 
T-L 7 





= 0. 





This proves that 

2 

ves 

i=l T, 
where the equality applies to reversible processes. The above 
argument can be extended to any cyclic transformation in which 
the system receives heat from or surrenders heat to a set of tem- 
peratures T,,T>,...,Ty (Fermi, 1936). If the amounts of heat 
exchanged are Q1,Qz,...,Qn (positive when heat is given to 
the system, negative otherwise) it can be shown that 

N 

D 


i=1 


alg 
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where again the equality applies to reversible processes. In the 
limit where N —> oo the above equation takes the form 

dQ 

— <0. 5.6 

7s (5.6) 
We may conclude that with the aid of the Carnot cycle we can 
introduce a quantity, which for reversible processes is an exact 
differential. This quantity, 6Q/T, is defined as dS, where S is a 
new state function called entropy. Since 6Q/T is, for reversible 
processes, an exact differential it follows that for any reversible 
transformation 1 — f 


f 
| oe a ee (5.7) 


In other words the change in entropy depends only on the initial 
and final states, not on the particular transformation. This result 
is fundamental in the formulation of the second law. 

Since S' is a state function it follows that entropy changes are 
due to changes in both temperature and volume (or pressure). 
Using the first law we have that for a reversible process 














ôQ dT dV 
dS = = Ov tP 
= dT dV 
aS = Cv -p + nk v 
or 
dS dT, nR*dV 
Cy T Vary. 
or 
dS dT dV 
Go Py 
or 
y-1 
S; = Si+ Cy ln (H) , (5.8) 


Equation (5.8) provides the change in entropy of an ideal gas 
(since the ideal gas law was used) as a function of the initial and 
final temperature and volume. This relation cannot, therefore, 
be used for liquids or solids. For liquids or solids Cy = C, = C 
and one can show (but we will not go into this proof here) that 
in this case the entropy change is dominated by the temperature 
change. As such for liquids or solids and for reversible processes 


Ty 


_~§=Cl 
$;-$.=Clng 


(5.9) 


Figure 5.2 

A cyclic transformation 
consisting of a reversible 
part from i to f(R) and 
an irreversible part from 


f toi(). 
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5.3 More on entropy 


As in the case of energy, we will choose an arbritary equilibrium 
state O as the standard state where by definition So is zero. Then 
we will define the entropy S4 of an equilibrium state A as 


4 5Q 
s= F 
A o T 


where the integral is taken over a reversible transformation. It can 
easily be shown that if instead of a standard state O we had chosen 
a standard state O’ then the entropy of state A would differ from 
the original one by an additive constant: from the above equation 


we have P 
ô 
E f ZA 
d T 
From equation (5.7) it follows that 
4, = Sa So 
or 
Sa = A = Sor. 


Since O’ is fixed, So, is some constant. The entropy is thus defined 
except for an additive constant. As in the case of energy this does 
not present a problem as long as we are dealing with differences, 
not actual values of entropy. Now let us consider the cyclic trans- 
formation shown in Figure 5.2. The path i ic f represents a 
reversible transformation and the path f eee represents an irre- 
versible transformation. As we know, for any cycle $ 5Q/T < 0. As 


such 5 
ME 
iRfIi T 


U FLAU F], 


or 
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Since i > f is reversible we have that 


al | 


Combining the last two equations leads to 


5,-5+|f 2 <0 
f T I 


ES 


By interchanging i and f we conclude that for any process 


or 


f 5Q 
or 
f dQ 
as> f = (5.11) 
or 
6Q 
> — . 
dia (5.12) 


where the equality applies to reversible processes only. Equation 
(5.12) is the general expression of the second law of thermody- 
namics and it indicates that the upper bound to the heat that 
can be absorbed by the system during a given change is ôQ = 
TdS. For a completely isolated system (Q = 0) equation (5.10) 
translates to 


ERA (5.13) 


This is an important conclusion as it indicates that for any spon- 
taneous irreversible transformation (i.e. one not related to external 
influences) occurring in an isolated system, the final entropy is 
greater than the initial entropy. It follows that when an isolated 
system has attained a state of maximum entropy it cannot undergo 
any further transformation because any change would decrease 
its entropy, and this is not allowed by equation (5.13). Thus the 
state of maximum entropy is a state of stable equilibrium. Note 
that if the isolated system consists of a number of microscopic 
subsystems it would be possible for some of the subsystems to 
reduce their entropy but only at the expense of the rest of the 
subsystems, whose entropy must increase enough that the over- 
all entropy of the system increases. A perfect isolated system is 
our universe (assuming that no other universes exist with which 
our universe can interact). It follows that the entropy of our uni- 
verse increases in time. Equations (5.10)-(5.12) as well as equation 
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(5.13) are the basic mathematical expressions of the second law of 
thermodynamics. 


5.4 Special forms of the second law 


Given the above definitions of the second law we can derive its 
expression for special cases such as the following. 


Finite isothermal transformations (AU = 0) 
From equation (5.11) it follows that in this case 


1 f 
ase a | 6Q 


or Q 
AS > = 
S2 T 
or i 
AS > —. .14 
S2 a (5.14) 
Adiabatic transformations 
Using equation (5.12) we obtain 
dS > 0. (5.15) 


Isentropic transformations 
An isentropic transformation is one during which the entropy 
does not change. In this case it is clear that (recall equation 
(5.12)) 

6Q <0. (5.16) 


Note that according to equation (5.15) a reversible adiabatic 
process is isentropic. 


Isochoric transformations 
From the first law we have that when dV = 0 then 6Q = CydT. 
As such for isochoric transformations 


dT 
dS > Cv 
or 
Ts 
Isobaric transformations 
In this case 
6Q = C,dT 
and it follows that 
dT 
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or 
T 
AS > C,ln (5.18) 


From equation (5.14) it follows that irreversible work increases 
a system’s entropy. From equation (5.17) it follows that in the 
absence of work the change in entropy depends on the relation 
between T; and T;. 


5.5 Combining the first and second laws 
Consider the following form of the first law: 
ôQ = C,dT — Vdp. 
Combining the above equation with dS > 8 yields 
TdS > C,dT — Vdp. 
Recalling that C, = dH/dT we can reduce the above equation to 
TdS > dH — Vdp 


or 


dH < TdS + Vdp. (5.19) 


Similarly, from 
dQ = CydT + pdV 


we arrive at 


TdS > dU + ôW 


or 


dU < TdS — pav. (5.20) 


It is often convenient to introduce two new functions: the Helmholtz 
function F = U — TS and the Gibbs function G = H-TS=U+ 
pV —TS. Since S = S(T,V) and U = U(T) it follows that both 
F and G are state functions (i.e. F = F(T,V),G = G(T,V)), and 
thus are exact differentials. The advantage of these functions is that 
they can be used to express equations (5.19) and (5.20) in a form 
where the pairs (T,p) and (T, V) appear as the independent vari- 
ables instead of the pairs (S,p) and (S,V). In this case equations 
(5.20) and (5.19) can be written as: 


dF < —SdT — pdV (5.21) 


and 
dG < —SdT + Vdp. (5.22) 


The interpretation of the above functions is that for isothermal 
processes dF < —ôQ or dF < —dW which makes F the energy 
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available for conversion into work. The usefulness of G becomes 
clearer for isothermal-isobaric transformations. During such trans- 
formations (which apply in phase transitions, for example water to 
vapor) dG = 0 and therefore G is conserved. Relations (5.19)—(5.22) 
are often referred to as the fundamental relations. 


5.6 Some consequences of the second law 


Thermodynamic definition of temperature 


As we know, the integral ff ÔQ depends on the path from i to f, 
not just the states 7 and f. In this chapter we showed that the 
integral fÍ 6Q/T depends only on the initial and final states, not 
on the particular path from i to f. As such 1/T is the integration 
factor that makes dQ an exact differential. It follows that in ther- 
modynamical terms temperature is the inverse of the integration 
factor of the differential of heat in reversible processes. 


The statistical nature of thermodynamics 


The change in entropy during a reversible cooling at constant 
pressure of m grams of dry air from T; to T, (Tı > T>) is 


T2 6 

AS = Q 
Tı T 
Since p = constant we have that 6Q = mcpadT where cpa is the 
specific heat capacity of dry air. It follows that 

Me dT 

AS = MCa — 

a ae 

or assuming that Cpa does not vary significantly with temperature 


AS = Sr, — Sr, = MCpa In ir 

Tı 

Since m and Cpa are greater than zero and In A < 0 it follows 
that AS < 0. Thus, the entropy decreases dune the cooling. This 
decrease does not violate the second law. Even though the dry air 
decreases its entropy the environment that provided the heat suffers 
a positive change so that the overall entropy change of the system 
is positive. We can extend the above arguments to show that as 
the temperature goes down so does the entropy (or as temperature 
goes up so does the entropy). The question is why and what is the 
physical meaning of this. Kinetic theory provides an explanation 
for the above. According to the theory, at very low temperatures 
the motion of the molecules is very slow with the molecules being 
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Figure 5.3 

The situation in (a) 
represents greater order 
and a lower number of 
complexions than in (b) 
where the number of 
complexions and disorder 
have increased. 


5 THE SECOND LAW OF THERMODYNAMICS 





| B 








(b) <`. AJ. B 














N = N/2 No=N/2 


more or less uniformly distributed in space. This is a picture of a 
high order, with complete order (motionless molecules uniformly 
spaced) presumably achieved at absolute zero. As the temperature 
increases the motion of the molecules increases and the order is soon 
destroyed. If we combine this with the above result it follows that we 
can associate entropy with order or disorder. Accordingly, decreas- 
ing entropy implies increasing order and increasing entropy implies 
increasing disorder. This association justifies the word entropy, 
which in Greek means the “inner behavior” (evtpomia). Now if 
such an association is correct we should be able to find a measure 
that defines disorder which should be mathematically related to 
entropy. 

Let us consider a volume V consisting of two compartments A and 
B, with N total particles moving around. Let us denote with N; 
the number of particles in compartment A and with N, the number 
of particles in compartment B. Statistics tells us that the number 
of different ways, P, in which N; particles will be found in A and 
N particles in B is 

N! 


P=: 
N! No! 


(5.23) 
Now consider two cases as shown in Figure 5.3. In the first case all 
particles are found in A (Figure 5.3(a)), and in the second case the 
particles are all over the box (Figure 5.3(b)). The first case (M1, = 
N, N, = 0) obviously represents a state of more order than the 
second where it can be assumed that N, = N = N/2. Applying 
equation (5.23) to both cases yields 


Case 1: 
N! 


= Noi 7! 
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Case 2: 
N! 


P = —_ 
NINJ” 


1. 

We thus see that we have a measure P (often called the number 
of complexions) that increases as disorder increases. If we consider 
the first case as the initial state and the second case as the final 
state we can conclude (since i > f is an irreversible process) that 
as the entropy increases P increases. So P is a possible candidate 
to relate disorder to entropy. But how can we find their functional 
relationship? 

The derivation of this relationship (one of the most beautiful yet 
simplest ones) is due to Boltzmann (Fermi, 1936). From the above 
we saw that as S increases, P increases and vice versa. So let us 
start with a general relationship of the form 


S= f(P). 


Next, consider a system consisting of two subsystems and let Sı 
and S> be the entropies and P, and P, the corresponding number 
of complexions. Then, Sı = f(P,) and S2 = f(P2). For the whole 
system we can then write that 


S = Sı + S2 


and 
P=P,P> 


or 
f(PiP2) = f( Pr) + f(P2) 


i.e. the function f obeys the functional equation 
f(xy) = F(a) + fy). 


Since the above equation is true for all values of x and y we may 
take y= 1 + € with e < 1. Then we can write it as 


f(x + ze) = f(x) + fl+e). (5.24) 


Expanding both sides of the above equation using Taylor’s theorem 
and neglecting terms higher than first order we get 


f(a) + zef (x) = f(@) + f0) + f0) (5.25) 


where f’ indicates the first derivative. For €e = 0 from equation 
(5.24) we have that f(1) = 0. Therefore, equation (5.25) reduces to 


we f(x) = ef'(1) 


or 
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Since f’(1) is a derivative evaluated at a value of 1 it is a constant. 
Thus 


ef'(e) =k 
or 
f=" 
or 
f(x) = T ‘ da 
or 


f(x) =k ln x+ constant 
or by changing z > P, 
S= k ln P+ constant. (5.26) 


This equation relates order and entropy and constitutes a funda- 
mental relationship connecting statistical mechanics and thermo- 
dynamics, thereby giving thermodynamics a statistical character. 


Entropy and potential temperature 
Recall that the potential temperature 0 is given by 
1000) */° 
jar (=») | 
p 
By taking the logarithms of both sides we obtain 


Indé=In T+ A ia 1000 — Hiap 
Cp Cp 
or 
cpd ln 0 = cpd In T — Rd ln p. (5.27) 
From the first law we have 


Q = C,dT — Vdp 








or 
ÔQ n dT dp 
T =C T Kr 
or A 
Q adr dp 
T = Cy T ni 


or (for reversible processes) 


dT d 
ds = e — Re (5.28) 


Figure 5.4 

A hypothetical upper-level 
flow and the trajectory 
followed by a parcel. 


5.6 SOME CONSEQUENCES OF THE SECOND LAW 71 


—zZ 


where ds is the specific entropy. Combining (5.28) and (5.27) 
leads to 


ds = cpd ln 0 (5.29) 


or 


s = cp ln 0 + constant. (5.30) 


Thus, except for an additive constant, the specific entropy of a sys- 
tem is given by the logarithm of the potential temperature. When 
0 remains constant the entropy remains constant. It follows that a 
reversible adiabatic process is isentropic. For an irreversible adia- 
batic process ds > 0,d0 = 0. In this case, the increase in entropy 
comes from irreversible work (for example, dissipation of kinetic 
energy due to friction). It follows that an isentropic process is 
adiabatic, but an adiabatic process may not be isentropic. 


Atmospheric motions 


Figure 5.4 shows a hypothetical upper-level flow close to some lat- 
itude A. Along this latitude the radiative-equilibrium temperature 
Tre(A) is assumed to be the same. Such thermal structure can 
be achieved if the motion is completely zonal (i.e. parallel to lati- 
tude circles), since air parcels then have infinite time to adjust to 
local thermal equilibrium. If at position 1 a parcel is allowed to 
move with the flow two possibilities exist. If the motion is very 
slow the parcel has time to come into equilibrium with its sur- 
roundings, so its temperature differs from Tre only infinitesimally. 
Therefore, we can assume that heat transfer along the trajectory 
of the parcel is reversible. Then for positions 1 and 2 we have 


that 
2 2 2 q 
odin = | as= f —. 
J j 1 1 T 
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If we assume that ôq depends only on the parcel’s temperature (for 
example, dg = Tf(T)) then we can write the above equation as 


[ c,d In 0 = a df(T) = f(T) — f(T). 


Since T = T; = Trp the above relationship reduces to 


2 
J cpd In 8 =0 
1 


A, = 0. 


or 


Thus, the parcel is restored to its initial thermodynamic state when 
it returns to the original latitude. Since 0 depends on p and p 
depends on height it follows that under the slow motion condi- 
tions, successive crossings of the initial latitude result in zero net 
vertical motion and zero heat transfer. 

On the other hand if the motion is fast enough that reversibility 
cannot be assumed, then 


2 
[cdo <ds 40 
1 


and therefore 62 Æ 6,. In this case when the parcel returns to the 
initial latitude it is displaced vertically and the net heat trans- 
fer is not zero. Vertical motion resulting from irreversibility in 
the atmosphere plays an important role in maintaining meridional 
circulations, heat and moisture transfer, etc. 


Examples 


(5.1) Imagine a volume V divided into two smaller volumes V; and 
V by a partition. Assume that m, grams of an ideal gas at 
temperature T, occupy volume V; and ma grams of another 
ideal gas at temperature T, occupy volume V2. Then the partition 
is removed. If the specific heat capacities at constant volume of 
the gases are cy, and cy,, find (a) the limiting temperature of the 
mixture, and (b) the change in entropy for each gas and for the 
system as a whole. 


(a) We assume that while the two gases exchange heat 
between them and perform work, the system as a whole 
does not exchange heat with its surroundings and its 
volume remains constant (V = V, + V2). In other words 
for the whole system the process of mixing is adiabatic, 
with no work being done. As such, AU = 0. The total 
internal energy before mixing is Cy, T, + Cv, T and 
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after the mixing it is Cv, T + Cy,T (here we have made 
use of the approximation U = CyT+ constant and 
then neglected the constant). Since AU = 0 it follows 
that 

Cy, Ti + Cr, To = Cr, T + Cy, T 


which yields 


Cy Ce 
T = 1 T. | 2 
Or tOu a Ont Oh 





Te (5.31) 


Considering that Cy, = Cy, M1 and Cy, = Cy, M2 the 
above equation becomes 
Cy, Mı 


1 2 
Cy Mı + Cy, M2 Cy, Mı + Cy, M2 


Cvp Mo 





(b) From equation (5.8) and the additive property of 
entropy it follows that for each gas the total entropy 
change is the sum of the entropy change due to 
the temperature change of the gas alone and the 
entropy change due to the volume change of the gas 
alone, i.e. 


AS, = ASır + ASiy 
AS», = ASor + ASoy 





where the lower limits for AS,7,ASjy,ASj7, and 
ASoy are given by (recall equation 5.8) 


T 
ASır = Cy, n — 
Sir Cy, Jr T, 


T 
ASər = Cy, In — 
Sor Cy, n T, 


Vi + Vo 


1 


ASiy = Cy, In ( 


y= 
Vi + V: 
) = mR, In ye 


1 


and 








y2—1 
y: 
Asv = Ov, m( 245) = Mm Rə mt 2 


2 2 


where R, and R, are the specific gas constants of the 
two gases. For the system as a whole the lower limit 
for the total entropy change is 


AS = AS, + AS» 
or 


AS = (ASir + ASor) + (ASiy + ASoy) = ASr + ASy 
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(5.2) 


(5.3) 


or 


T T 
AS = levi In — + cym In | 


Tı To 
(5.32) 
+ m,R, eet + mR, matk P 
Vi Vo 


In equation (5.32) the second term on the right-hand 
side (AS) is necessarily greater than zero. The first 
term (A Sr) is either positive or zero when T; = T (see 
problem 5.13). As such for the whole system AS > 0 
which according to (5.12) indicates irreversible process. 
This proves that mixing is irreversible. 


Calculate the change in air pressure if the specific entropy 
decreases by 0.05 J g-'K"' and the air temperature decreases 
by 5%. 


Assuming that the process is reversible we have 


dT V 
dT „dp 


Since we do not have any information on the mass of the 
air involved in the process we will write the above per unit 


mass: dT d 
Pp 
ds = C — a 


or 


dp _c,dT ds 
p RT R 





or considering the air as dry 


in Ps — 1005, (0.957: _ (50) 
pi 287 T, 287 





or 


ps = 0.994 P, 


The pressure will decrease by a small amount of about 0.6%. 


In example 4.3, what is the change in entropy in each of the 
three transformations? 


(a) The first transformation is an isothermal reversible 
compression. Thus, 
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AS = [?- 5 [60-3 = nein snr ak 
V; Pf 


(b) The second transformation consists of two branches. 
The first one is a reversible adiabatic transforma- 
tion. Thus, the entropy change in the first branch is 
zero. The second branch is an isobaric reversible com- 
pression. As such the entropy change of the second 
transformation is (recall that since T; and Ty; lie on 
the same isothermal T, = T; = T) 


as= f Pool TG i =e 


= Cp ln eu 


= —nR* ln Pr 
Pi 


= n R* ln z 
Pf 





(c) The third transformation also consists of two branches: 
one reversible isochoric and one reversible isobaric. The 
total change in entropy is the sum of the changes in 
entropy in the two branches 


-{% +f[ 8 
T Ja T 


dT 
T w + [ GF 


+ T; 
+ Cy Ing 


=Cy ln 








Tin! 


During this transformation p;V; = nR*T, and p;V; = 
NR*T py! - Thus 


Tn! = Pir, 


If we substitute the above value of Tm in the previous 


equation we obtain 


T : 
AS = Cin Ona 

Pi T, py 
Since T; = T; = T it follows that 


AS = Cy In f +C, n Č 
Pi Pf 
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or ý: 
AS = nR* ln —. 
Pf 
Thus, in all transformations the change in entropy is the 
same. 


Problems 


(5.1) 


(5.2) 


(5.7) 


(5.8) 


Starting from the first law show that 6Q/T is an exact 
differential. 
Derive the so-called Maxwell relations for reversible pro- 


5-8, 
D-8. 
ORO) 
D-0, 


Find the lower limit for the change in specific entropy of dry 
air at a given temperature and pressure which is allowed to 
expand freely into an insulated chamber to assume twice its 
original volume. (200 J K~' kg~*) 

What is the form of isenthalpic lines in a (p, V) diagram? 
Dry air at T=0°C and p=1 atmosphere is compressed isen- 
tropically until its pressure becomes 10 atmospheres. Find 
the final temperature. (527.5 K) 

In a Carnot cycle the area for each transformation corre- 
sponds to the following work — (a) isothermal expansion: 
31165J, (b) adiabatic expansion: 21517J, (c) isothermal 
compression: 24282 J, (d) adiabatic compression: 21517 J. 
Calculate (1) the heat provided, per cycle, to the gas by 
the warm source, and (2) the efficiency of the cycle. (6883 J, 
0.22) 

One mole of ideal gas occupies 5 liters at 273 K. If it expands 
into a vacuum until its volume becomes 20 liters, what is the 
lower limit for the change in entropy and the upper limit for 
the change in Gibb’s function? (11.5 J K+, —3146 J) 

A fixed volume 2V is divided into two equal volumes V 
separated by a diaphragm. Each volume contains the same 
amount of the same gas at the same temperature T. Using 
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Boltzmann’s expression for entropy, show that after the 
diaphragm is removed the entropy change is zero. How 
will you approach this problem if the temperature in the 
two volumes is different? Use the approximation In N! = 
NinN-N. 

(5.9) In a thermodynamic process a parcel of dry air is lifted 

in the atmosphere such that its pressure decreases from 

1000 to 800mb while its temperature remains constant. 

Calculate the change in the specific entropy of the parcel. 

(64 Jkg~* K~*) 

During a process a parcel of dry air descends from 900 

to 950mb and its specific entropy decreases by 30Jkg~! 

K~. If its initial temperature is 273K what is (a) its 

final temperature and (b) its final potential temperature. 

(269 K, 273 K) 

(5.11) Show that for flows on a constant pressure level the relative 
change (i.e. percent change of original value) in temperature 
is equal to the relative change in potential temperature. 

(5.12) Cooler air moves over a warmer surface. What will happen 
to the surface pressure? (Hint: the condition of problem 5.11 
does not apply here.) 

(5.13) Show that the first term on the right-hand side of equation 

(5.32) is greater than zero except when T, = T, when it is 

zero. (Hint: define b = Cy, /Cy, and z = T,/T, and find the 

first and second derivatives of f(x) = AS7/Cy,.) 

Try to explain to your friends why the second law of ther- 

modynamics implies a positive arrow of time (i.e. time goes 

only forwards). 

(5.15) Consider two isolated systems existing very far from each 
other. The same irreversible transformation occurs in both 
of them. Will the change in entropy be the same in both of 
them? 

(5.16) In an isentropic process the specific volume of a sample of dry 
air increases from 300cm*g~* to 500 cm? g~t. If the initial 
temperature was 300K, what are the final values of T and 
p? (244.5 K, 1404 mb) 

(5.17) Prove geometrically (in a p,V diagram) equation (5.2). 


(5.10 


m= 


(5.14 


ar; 


CHAPTER SIX 





Water and its 
transformations 


The fundamental equations derived in the previous chapter can only 
be applied to closed systems (i.e. to systems that do not exchange 
mass) that are homogeneous (i.e. they involve just one phase). In 
such cases we do not need to specify how thermodynamic functions 
depend on the composition of the system. Only two independent 
variables (T and p, or p and V, or p and T) need to be known. 
Since the total mass (m) remains constant, if we know the values 
of the extensive variables per unit mass we can extend the equa- 
tions to any mass by multiplying by m or by n (the number of 
moles). 

A heterogeneous system involves more than a single phase. In this 
case we are concerned with the conditions of internal equilibrium 
between the phases. Even if the heterogeneous system is assumed to 
be (as a whole) a closed system the phases constitute homogeneous 
but open “subsystems” which can exchange mass between them. In 
this case the fundamental equations must be modified to include 
extra terms to account for the mass exchanges. These extra terms 
involve a function u called the chemical potential, y = u(T, p). 
We will not go into the details of defining u; we will only accept 
that in the case of open systems something else must be included 
to account for the heterogeneity of the system. In this book we 
are concerned with a heterogeneous system that involves dry air 
(No, O2, COz, Ar) and water, with the water existing in vapor 
and possibly one of the condensed phases (water or ice). Then one 
component of the system is the dry air (which is assumed to remain 
always unchanged and in a gaseous state) and the other component 
is the water substance which can exist in two phases. Thus dry air 
is a closed system whereas the two phases are two open systems. 
The system composed of dry air, water vapor, liquid water, and ice 
is also of interest but it is quite unstable and as a result it cannot 
exist in equilibrium. 
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6 WATER AND ITS TRANSFORMATIONS 


6.1 Thermodynamic properties of water 


In this chapter we will put aside dry air and we will concentrate 
on the one-component heterogeneous system “water” comprised of 
vapor and one of the condensed phases. Water vapor, like dry air, 
can be treated to a good approximation as an ideal gas. Accordingly, 
water vapor obeys the equation of state and if it existed only by 
itself its states would be determined by two independent variables. 
When, however, water vapor coexists with liquid water or ice things 
become a little complicated (or less complicated depending on the 
point of view). In this case, the mixture does not constitute an ideal 
gas and the equations applicable to ideal gases do not apply. One 
phase requires two independent variables (p and T) and the other 
phase two more (p’ and T’). Then, equilibrium between the phases 
will require that p = p' and T = T’. However, since the mass of 
each open subsystem does not remain constant another criterion for 
equilibrium must be considered. This criterion is given by u = p’ 
and expresses equilibrium in the terms responsible for the mass 
exchange. Accordingly, for equilibrium between the two phases the 
following three constraints must be satisfied: 


p=p 
T=T' 
=w. 


The first two equations reduce the number of independent variables 
to two (say p and T). Then, the third equation w = p(p,T) = 
u'(p',T') = p'(p, T) reduces the number of independent variables by 
one more. Since u(p, T) = u (p, T) it follows that p = f(T). Thus, 
if we fix the temperature at which the phases exist in equilibrium 
the value of pressure also becomes fixed. This defines curves along 
which equilibrium between two phases can exist. We conclude that 
in the case of a one-component system involving two phases in 
equilibrium the number of independent variables (otherwise known 
as degrees of freedom) is one, not four. 

Let us now consider the case of a one-component system involving 
all its phases (vapor, liquid, solid). Following the same arguments 
as above the following constraints must be satisfied for equilibrium 


p=p = p" 
T = T’ = T” 
=K =i". 


Then the first two equations reduce the number of independent 
variables from six to two. The first equality of the third equation 
reduces them by one more and the second equality by yet one more. 


Figure 6.1 
Phase-transition equilibria 
for water. 
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Thus, we remain with zero independent variables, meaning that all 
values are fixed, which implies that coexistence of all phases at equi- 
librium requires just a point, called the triple point, in the p, V, T 
state space. It follows that the greater the number of phases the 
fewer the degrees of freedom. The following formula suggested by 
Gibbs generalizes the above as it gives the number of independent 
state variables, N, for a heterogeneous system involving C different 
and non-reactive components in a total of P phases: 


N=C+2-P. (6.1) 


Figure 6.1 shows the triple point for water as well as the curves 
p = f(T) representing the phase transition equilibria. We call 
the three equilibria vaporization (gas + liquid), fusion (liquid + 
ice), and sublimation (gas + solid). Along the curves for vaporiza- 
tion and sublimation, vapor is in equilibrium with water and ice, 
respectively. As such these curves provide the equilibrium vapor 
pressure for water and ice. Vapor in equilibrium with a condensed 
phase is often called saturated and the corresponding equilibrium 
vapor pressure is called saturation vapor pressure. There is no real 
distinction between the two terms. Note that the extension of the 
vaporization curve to temperatures below the triple point corre- 
sponds to supercooled water, with vapor pressure over supercooled 
water being greater than that over ice. This is a metastable equi- 
librium where supercooled water and water vapor coexist. Here the 
system of supercooled water and water vapor may be stable with 
respect to small changes in temperature and pressure but intro- 
duction of ice in the system makes it unstable and as a result water 
freezes. The triple point corresponds to p, = 6.11 mb and T = 273K 
(more accurately 273.16 K). Because of the different densities of 
liquid water, ice, and water vapor, the specific volumes for liquid 
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water, dy, ice, a;, and water vapor, ay, at the triple point are: 


dy = 1.000 x 107? m? kg" 
a; = 1.091 x 107? m*kg™* 
dy = 206 m? kg. 


An interesting feature of Figure 6.1 is that the vaporization curve 
ends at point C where the temperature (To) is 374°C and the pres- 
sure (pc) is approximately 2.21 x 10° mb. Beyond this critical point 
there is no line to separate liquid from vapor. Otherwise stated, 
beyond point C there is no discontinuity between the liquid and 
the gaseous phase. This indicates that beyond point C we cannot 
differentiate between the liquid and the vapor phase. In order to 
understand the significance of critical point C let us consider point 
A representing the temperature and the pressure of a sample of 
water vapor below critical point C. If we decrease the volume while 
keeping the temperature constant the pressure increases until it 
becomes equal to the pressure at point A, where it will condense. 
Similarly, if we keep the pressure constant and cool the gas until its 
temperature becomes equal to the temperature corresponding to 
point A», the gas will condense. Obviously the gas will condense if 
we perform any intermediate process described by line AA. Thus, 
there are three ways by which the gas can become liquid. We either 
have to cool it at constant pressure or compress it at constant 
temperature or simultaneously compress it and cool it. All three 
ways require crossing the equilibrium line. The fact that there is 
an end to the vaporization curve means that above Te we cannot 
liquefy a gas by compressing it at a constant temperature and that 
above po we cannot liquefy a gas by cooling it at a constant pres- 
sure. Because in our atmosphere temperature and pressure are well 
below the critical point, water vapor can condense. For the other 
gases in the atmosphere (N2, O2, Ar) the critical temperatures are 
very low (—147 °C, —119°C, —122°C, respectively) and so they do 
not condense. 

Let us now consider phase changes along isotherms in the (p, V) 
domain. Such a domain (shown in Figure 6.2) is called an Amagat- 
Andrews diagram. We start with a sample of water vapor at a state 
corresponding to point A in Figure 6.1 (i.e. at a temperature T, and 
pressure pı greater than the triple point temperature T, and pres- 
sure p,. If we compress the vapor isothermally the pressure increases 
until we reach point A. At this point liquid water and water vapor 
coexist in equilibrium. This means that some water vapor has con- 
densed to form liquid water and that liquid and vapor are in equilib- 
rium. Now what will happen if we further compress the vapor? Since 
we are on the equilibrium curve the pressure of the vapor depends 


Figure 6.2 
An Amagat—Andrews 
diagram. 
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only on temperature. As such, since the transformation is isother- 
mal, as long as liquid and vapor coexist a further decrease in volume 
cannot change the pressure. Do not get confused here. The ideal gas 
law does not apply to the mixture. Thus, if we keep on compress- 
ing the vapor we will observe that more vapor condenses to form 
liquid water (with liquid and vapor being always in equilibrium) 
until we reach point B where all vapor is condensed. This part will 
be represented in the (p, V) diagram by the horizontal (constant 
p) stretch A,B. A further (small) compression causes the liquid 
water’s pressure to increase rapidly (owing to the low compress- 
ibility of liquids). This part is indicated by the segment BB,. The 
complete isotherm is depicted in Figure 6.2 by the line AA, BB,. 
For higher temperatures we observe that the horizontal stretch of 
the isotherms decreases until it reduces to a point. This is the point 
C corresponding to the critical state Te, pc, Vo. The isotherms 
above the critical temperature are monotonic decreasing functions 
with no discontinuity. They go over into equilateral hyperbolas 
which is how ideal gases behave. For temperatures below T, the 
horizontal stretch of the isotherms increases. A solid line in Fig- 
ure 6.2 connects the beginning and end points of the horizontal 
stretches. This line together with the critical isotherm and the 
isotherm corresponding to the triple point temperature (T,) par- 
tition the diagram into six regions: vapor, gas, liquid water, ice, 
liquid water plus vapor, and ice plus vapor. Here the distinction 
between gas and vapor reflects only the fact that above To no 
two-phase discontinuity for condensation exists. 


6.2 Equilibrium phase transformations — latent heat 


In the case of a homogeneous system undergoing an isobaric 
transformation that does not result in a phase change the heat 
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exchanged is proportional to the temperature change (the 
proportionality constant is Cp). In the case of a heterogeneous sys- 
tem involving two phases we say that when the two phases are in 
equilibrium a fixed temperature implies a fixed pressure, meaning 
that isobaric transformations are also isothermal transformations. 
Thus, even though the masses of the two phases change, the tem- 
perature of the system does not (recall the horizontal stretches in 
Figure 6.2). In such cases, the amount of heat exchanged depends 
only on mass changes (which modify the internal energy of each sub- 
system) and on the work done owing to possible volume changes 
(the volume can change because the mixture does not obey the ideal 
gas law). By definition, the latent heat, L, of a transformation is 
the heat absorbed (or given away) by the system during an isobaric 
phase transition 
L = 6Q, = constant. 


Recall that enthalpy is defined as H = U + pV. Differentiation of 
this equation gives dH = dU +pdV +V dp. Since the change of phase 
takes place at a constant pressure it follows that dH = dU + pdV. 
The right-hand side of this equation is according to the first law 
equal to 6Q. It follows that 


L=dH. (6.2) 


Thus, the latent heat of a phase change is the change in enthalpy 
during the transformation. The latent heats of vaporization (liquid 
< vapor), fusion (ice © liquid), and sublimation (ice vapor) 
are denoted by Ly, Lẹ, and Ls, respectively. Note that the above 
latent heats are positive when during the transformation heat is 
absorbed (liquid — vapor, ice — liquid, ice — vapor) and negative 
when heat is released (vapor — liquid, liquid — ice, vapor —> ice). 
L,, Lẹ, and L, are also called the enthalpy of vaporization, fusion, 
and sublimation, respectively. According to the above, 


L, = H, ~ A, = U, BA Uw + Paw (Vo = Vw) 
L; = Hy — H; = Uy — U; + pwi(Va — Vi) (6.3) 
Ls a H, — H; = U, — U; + pv (Ve — Vi) 


where here p stands for the equilibrium pressure between two 
phases. At the triple point pyy = Pvi = Pwi = 6.11 mb. It follows 
that at the triple point the specific latent heats l,,l¢,and l, satisfy 
the relationship l, = le + ly. Obviously, the latent heats depend on 
temperature (see Table A.3). Nevertheless, in the range of temper- 
atures observed in the troposphere they do not vary significantly. 
Because of that it is often assumed that the latent heats are inde- 
pendent of temperature. It is easy to show that the first law in 
the case of an isobaric change of phase can be written as dU = 
L — pdV. For fusion dV ~ 0 and for vaporization and sublimation 


Figure 6.3 

The isotherm T; of Figure 
6.2. For any given change 
C — C” the change in 
volume is dV. This causes 
an amount dm of liquid 
water to evaporate. 
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dV = Vyapor— Viiquid water or ice % Vvapor- Thus, the first law reduces to 


dU = L for fusion 
dU = L—mR,T for vaporization or sublimation (6.4) 


where R, is the gas constant for water vapor. 


6.3 The Clausius—Clapeyron (C-C) equation 


Before we proceed and to avoid confusion with too many subscripts 
we need to make a change in our notation. From now on we shall 
represent the water vapor pressure as e. Furthermore, during phase 
transitions the equilibrium (saturation) vapor pressure over water 
will be denoted as e,,, and equilibrium vapor pressure over ice will 
be denoted as esi. 

For an ideal gas we saw how the equation of state p = f(V,T) 
relates changes of pressure, temperature, and volume. In the case 
of a heterogeneous system involving two phases we saw that p = 
es = f(T). Does this indicate that a relationship similar to the 
equation of state for ideal gases exists for heterogeneous systems? 
Let us consider the transformation described by isothermal T, in 
Figure 6.2 and focus our attention on segment A,B. Going from 
A, to B a change of phase takes place (vapor — liquid). As we 
discussed before, between A, and B vapor is in equilibrium with 
liquid water and p = esw = f(T). Because of that any other state 
variable or function is a function of temperature only. Let us assume 
that at point C (Figure 6.3) the mass of the liquid phase is m,, and 
the mass of the vapor phase is m,. Then the total volume and 
internal energy of the system at point C are 


V=m,ay, + my,ay 





U = MyUw + My Uy 








Y 





86 


6 WATER AND ITS TRANSFORMATIONS 


where dy, ay, Uw, Uy are the specific volumes and internal energies 
of the two phases. Now assume that the system moves from C to 
C”. This change corresponds to a change in volume dV and causes 
an amount dm of liquid to evaporate. The volume at point C” will 
be 

V +dvV = (mwy — dm)a, + (m, + dm)ay 


or using the previous equation 
dV = (a, — a,,)dm. (6.5) 


Similarly 
dU = (uy — uy)dm. (6.6) 


Recall from equation (6.3) that during this change 
Uy — Uw + €sw(Ay — dy) = |. (6.7) 
Dividing equation (6.6) by (6.5) yields 





a _ wu 
dV ay — Gy 
or using equation (6.7) 
dU ly 





a= ee Caw: (6.8) 


Here is as good a place as any to restate that while general defini- 
tions such as the first law (equation (4.5)), the second law (equation 
(5.12)), enthalpy (H = U +pV), specific heat (equation (4.14)), etc. 
are valid for any system, variations or expressions derived using the 
ideal gas law (for example, equations (4.22) and (4.18)) are not valid 
for any system but ideal gases. According to problem 6.1, for any 


system 
OU Op 
sear) ER ee ee 
(a)r (ar), -? i 


Also, according to problem 6.2, for an ideal gas (OU/OV)r = 0 
(Joule’s law). However, water in equilibrium with its vapor is not 
an ideal gas and U and esw are functions of temperature only. Then 
since the transformation is isothermal we can write that 


au) _ du 
wv), dV 


esw) _ desw 
Tj, ge 


Accordingly, equation (6.9) becomes 


dU desw 
4V = T dr "Cau 


and 
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Combining this equation with equation (6.8) yields 





desw l, 
T Sw = Sw 
dT z Ay — Aw . 
or 
desw by 





= iy l 
dT T (ay — aw) (on) 
Equation (6.10) is called the Clausius—-Clapeyron (C-C) equation. 


In its general form, 

des l 

= 6.11 

dT TAa’ (PAL 
it relates the equilibrium pressure between two phases to the tem- 
perature of the heterogeneous system. Here l is the specific latent 
heat (or specific enthalpy) of the change of phase and Aa is the 
difference in the specific volumes between the two phases at tem- 
perature T. It is for heterogeneous systems what the equation of 
state is for ideal gases. 








6.4 Approximations and consequences of the C-C 
equation 


e Temperature dependence of enthalpy of vaporization 


By differentiating equation (6.3) we have 


OL, = OH, OH, at JO 

oT OT OT 
where Cpv and C,, are the heat capacities at constant pressure 
for vapor and liquid water, respectively. Since Ly is the enthalpy 
difference between the vapor and liquid water phases at equilibrium, 
it depends only on temperature. As such the above equation can 
be written as 





dL, 

dT 
For a wide range of temperatures (—20°C to 30°C), C,, and Cpw 
vary very little (about 1%). Because of that we can consider them 
independent of temperature and integrate equation (6.12) to get 





= Cp — Cow- (6.12) 


L, = Lyo a (Coy a Cow) (T ~~ To) 


or 
ly = lvo + (cp er To) (6.13) 


where lo is the specific latent heat of vaporization at the ref- 
erence state To. For the reference state To = 273K with ho = 
2.5 x 10° The cp = 1850 Jkg™™ Kt, and cpw = 4218 Jkg™t K™}, 


P 
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equation (6.13) provides a very good approximation of l, over the 
range from —20°C to 30°C. 
We can repeat the above for sublimation and fusion to arrive at 


dL, 





aT Cov — Opi 
dL 
or 
l = lso a (Cov = Cri) (LT = To) (6.14) 








le = leo (Cow Cpi) (T To). 


Recall that the specific heat for constant volume is cy = ôq/dT. 
Since water vapor is considered an ideal gas (recall Joule’s law), it 
follows that 

_ oq duy 

= d dT? 

For liquid water (not an ideal gas) 





Cyy 


_ ôq _ duw | day, 
“w= ae aT are 


Since a,, varies very little with T the above equation reduces to 





S duw 
CVw = dT ʻ 


Similarly, for ice cy; = du;/dT. Now since h = u+ pa it follows that 
for water vapor cpv = dh, /dT = du,/dT + Ry. Therefore for water 
vapor Cp # Cyy. Liquid water and ice do not behave as ideal gases. 
Therefore pa # RT. In this case 





dh, du, daw dp duw 
= = t F aw x 
w = aT at PaT dt dT 
(because dp = 0 and da,,/dT ~ 0). Thus, cpw © Cyw = Cw. Similarly 
Cpi © Cvi = G. Using the values of Cp, Cvv,Cw, and q at T = 


0°C (1850, 1390, 4218, and 2106 Jkg~* K-t, respectively) and using 
equations (6.13) and (6.14) we find that 





dl 
— = —2368Jkg-1 K7! 
dT 8 
dl 

© x —256 Jkg 1 K7! 
dT 8 
dls LAN 
— y 2112Jkg™! Kt. 
dT Tkg 


These rates are small compared with the values of l,,1,,and lẹ at 
0°C (2.510%, 2.834x 10°, and 0.334x10°J kg~' K-t, respectively). 
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We can thus conclude that all latent heats vary very little with 
temperature. As such the above equations justify approximating 
l,,/;, and lẹ as constants. An improvement to this approximation is 
given by equations (6.13) and (6.14) which express them as linear 
functions of T. 


e Temperature dependence of equilibrium (saturation) vapor 
pressure 


In the case of vaporization, a, >> aw and C-C equation can be 


approximated as 

desw ly 
= : 6.15 
dT Ta, ( ) 
By combining equations (6.15) and (6.13) and considering vapor as 
an ideal gas (i.e. €swa@y = R,T) we obtain 








k desy -lso Ge Cr) Io. Cow ey 











en dT R,T? R,T 
or 
B AT de. lo + (Cpw 7 Cov) Lo [ dI — Cpw — Cp f dT 
Lanse 7 R, Pee R, m T 
or 





In Cage ko (Coat py) To ( 1 1 ) Cow — Cpv l T 


€so = Ry To T Ry š To f 


Considering as the reference state the triple point where Tọ = 
0°C, es = OA mb, lo = 2.5 x 10êJkg™t, cp = 1850Jkg™' 
K-t, and cpw = 4218 Jkg~* K~! the above equation reduces to 





e 1 1 T 
l sw a E Les l S 
n 2 6808 (z +) 5.09 In T, 


or 
6808 
esw = 6.11 exp (5349 re 5.0917 ) , (6.16) 


The above equation provides the relation between e,, and T. 
Another way to establish such a relation would be to integrate equa- 
tion (6.15) assuming that l, is constant and independent of T. This 
would result in an approximation to the above equation given by 


Csw l, ly 


€so ~ RyTo RT 





ln 


or 


5417 
esw = 6.11 exp (19.88 = =) (6.17) 


In both equations (6.16) and (6.17) esw is in mb and T is in Kelvin. 
In problem 6.3 you are asked to show that equations (6.16) and 
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(6.17) are nearly identical for temperatures in the range —20°C to 
30°C. 

The corresponding relation to equations (6.16) and (6.17) for subli- 
mation can be derived following the same procedure (here again 
a, > a) with I, = 2.834 x 10°Jke~? and c = 2106Jkg™ 
K~!. They are: 


2 
esi = 6.11 exp (2516 — 0293 — 0.555 In r) (6.18) 
and 
142 
esi = 6.11 exp (22.49 — =) ; (6.19) 


The C-C equation and its simplified versions describe the thermo- 
dynamic state of water vapor and liquid water (or water vapor and 
ice) when the two phases are in equilibrium. When the system’s 
pressure is greater than the equilibrium pressure, equilibrium is 
approached by condensation (sublimation) of water vapor. This 
reduces the amount of vapor in the system thereby reducing the 
vapor pressure. If the system’s pressure is lower than the equilib- 
rium pressure, equilibrium is approached by water (ice) evaporating 
(sublimating). Because of the particular exponential dependence of 
esw and e, on T in equations (6.17) and (6.19) it follows that for 
higher temperatures esw is greater than for lower temperatures. As 
such more water vapor can exist at higher temperatures. For exam- 
ple, for T = 300K from equation (6.17) we have esw = 36 mb. For 
T = 273°C, esw = 6mb. This is a substantial difference. Note that 
this does necessarily imply that warmer air can hold more water 
vapor than colder air as it is often wrongly stated. All the above 
arguments and formulas were derived in the absence of air! We only 
considered two phases with water evaporating (or vapor condens- 
ing) into a vacuum. If you think about it this is consistent with 
Dalton’s law which would indicate that in a mixture of air and 
water vapor the partial pressures are independent of the mixture. 
Therefore, keep in mind that when we say that “air at a given tem- 
perature is saturated” we do not imply that the air holds as much 
vapor as it can hold. We strictly mean that at that temperature 
the amount of vapor that can exist is maximum regardless of the 
presence of air. 

The reason for the relation between equilibrium vapor pressure 
and temperature arises from purely kinetic reasons. At equilibrium 
evaporation is zero (strictly speaking it is the net evaporation = 
evaporation minus condensation that is zero as molecules inside 
water have as much chance to evaporate as vapor molecules have 
to condense). At higher temperatures the molecules inside the liq- 
uid acquire greater speeds and their chance to escape increases. As 
a result (net) evaporation increases. 
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The exponential dependence of esw (or esi) on T has important con- 
sequences for our climate system. Tropical warm water can transfer 
more water vapor into the atmosphere than colder extratropical 
water. Most of this water vapor is then precipitated as a result of 
organized convection in the tropics thereby releasing large amounts 
of heat. This heat can then be converted to work and generate 
kinetic energy which helps maintain the general circulation against 
frictional dissipation. 


e Changes in the melting and boiling points for water 


At boiling point (T = 100°C), l, = 2.26 x 10°Jkg™*, a, = 
1.673 më kg~*, and aw =0.00104 më ke~*. Then from equation (6.15) 


it follows that 
desw 


= 3621 Pa K7! 
IT 36 a 





or 
desw 


dT 


Accordingly, a decrease in pressure by 0.03575 atm (36.2 mb) lowers 
the boiling point by 1°C. Note that if we evaluate the approximate 
relationship (6.16) at T, = 373K and T, = 372K, we will find 
that esw(Tı) 1003 mb and esw(T2) ~ 969 mb, or Ae, ~ —34 mb. 
If we use instead equation (6.17), we will find that e.,(Ti) ~% 
1237 mb, esw(T2) ~ 1190 mb, and Ae,,, ~ —47 mb. The differences 
are due to the fact that in deriving the approximate expressions 
we have assumed that the specific heats and latent heats are either 
constant or linear functions of T. While in the range of temper- 
ature relevant to weather these approximations may be good, for 
temperatures outside this range the differences may be significant. 
At melting point (T = 0°C), l = 0.334 x 10°Jkg™’, a, = 
1.00013 x 107? m?kg~*, and a; = 1.0907 x 107? m?kg~*. Substi- 
tuting these values into equation (6.15) yields 


= 0.03575 atm K~?. 








dpwi 
w=-1 PaK“! 
IT 3 503 800 Pa 
or F 
Pwi Ea —1 
ar ~ 133.3atm kK. 


This result indicates that lowering the melting point by just 1°C will 
require a tremendous increase in pressure. This is actually rather 
anomalous as in most cases the melting point increases with increas- 
ing pressure. The strange behavior of water and ice is due to the 
fact that ice is less dense than water whereas in most cases the solid 
is denser than the liquid. The fact that the melting point of ice is 
lowered by increasing pressure is very important to geophysics as 
it explains the motion of glaciers. When a mass of ice encounters 
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a rock on the glacier bed the high pressure of the ice against the 
rock lowers the melting point of ice causing the ice to melt and 
then to flow around the rock. Once behind the rock the pressure 
is restored and ice freezes again. In this way ice propagates around 
obstacles. 


Examples 


(6.1) 


Calculate the work done, the amount of heat absorbed, and the 
change in the internal energy during (a) melting of | g of ice to 
water at a temperature of 0°C and under a constant pressure 
of | atmosphere, and (b) evaporation of | g of water to water 
vapor at a temperature of 100°C and under a constant pressure 
of | atmosphere. 


(a) When the effect of absorbing an amount Q at a con- 
stant pressure is a change in the physical state of the system 
(which also occurs at constant temperature), then Q is pro- 
portional to the mass that undergoes the change. In our case 
Q = ml; where I; is the specific latent heat of fusion of 0°C. 
For ice lp = 79.7calg~'. Thus, 


Q = lem = 79.7 cal. 


During melting the volume changes. The corresponding work 
done is 


water water 
W = j pdV = pf dV = DP(Vivater a Vice). 


ice 


The density of ice at 0°C is 0.917 gcm~*, so 1g occupies a 
volume of 1.09 cm?. The density of water at 0°C is 1 gcm™?. 
Accordingly, 1g of water occupies a volume of 1.0cm?. It 


follows that 
W = 1013 x 10? x (1 — 1.09) x 107° J 


or 
W = —0.00912 J 


or 


W = —0.0022 cal. 


The work done is negative because the volume of the result- 
ing water is smaller than the volume of ice. From the above 
we then find that 


AU = Q — W = 79.7022 cal. 
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(6.2) 


(b) In this case at 100°C, 1, = 539cal g~t. Thus, 
Q = l,m = 539 cal. 


Water at 100°C has a density of 0.958 g cm~? which for 1g 
gives a volume of 1.04cm?. The volume of 1 g of water vapor 
at 100°C is 1673cm?. As such 


W = p(Veapor — Vwater) = 1013 x 10? x (1673 — 1.04) x 107° 
= 169.37 J = 40.47 cal. 


It then follows that 
AU = 498.53 cal. 


Note that in both cases the internal energy increases even 
though the temperature remains constant. This is due to 
the fact that the degrees of freedom increase. The degrees 
of freedom of the molecules in H2O is much greater in the 
liquid phase than in the ice phase and much greater in water 
vapor than liquid water. 


Can raindrops grow from cloud droplets through condensation 
alone? 


In order to address this question we have to make some rea- 
sonable assumptions. Cloud droplets form over nuclei and 
begin to grow according to the available water vapor. If the 
vapor pressure is greater than the equilibrium vapor pres- 
sure droplets grow by condensation. So we need to get an 
idea of how much vapor is available and how many droplets 
are going to compete for this vapor. Let us assume that 
the concentration of nuclei is 100cm~* which corresponds 
to rather clean air. Now we need to estimate the amount 
of vapor that will be available during the process of grow- 
ing by condensation. We can get a gross estimate of this 
amount by assuming that at the level where clouds form (say 
1000 m) the air is saturated and the surface temperature is 
about 30°C. If this air ascends adiabatically, the tempera- 
ture decreases. As the temperature decreases the equilibrium 
(saturation) vapor pressure decreases. This means that not 
as much vapor is needed for equilibrium or that the vapor 
pressure becomes greater than the equilibrium vapor pres- 
sure. Thus, the droplets begin to grow. If the air ascends to 
the maximum level (say the top of the troposphere at about 
16 km), then we can estimate the equilibrium vapor pressure 
at that level by assuming that no water vapor condenses 
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(6.3) 


using the equation 


e 5417 
ln -= = 19.83 — —— 
Eo AASA 
where T = 30 — (16 x 9.8) = —126.8 °C = 146.2 K. It follows 
that at this level e,,, = 0. At cloud base where the temper- 
ature is % 20.2°C, esw & 23.7 mb. As such in this scenario 
the equivalent of 23.7mb in water vapor will be available 
to the droplets for growth. How much water vapor is that? 
From the ideal gas law we find that for V = 1m? 

esw V 


m = iy = 0.0175 kg 








Assuming that a typical raindrop has a radius of 0.5mm it 
follows that the mass of a raindrop is Myaindrop = V/p ~% 
0.00052 g. Given the total amount of water vapor available 
this translates to ~ 33650 raindrops. However, in 1 m? of air 
there are 108 “competitors” with equal rights. Simply, none 
of them will grow to be a raindrop by condensation only. 
Other processes must also take place (collision, coalescence, 
for example). 


Consider a sample of 2 moles of supercooled water (i.e. liquid 
water existing at environmental temperatures less than 0 °C). 
Assume that the temperature of the surroundings is —10 °C 
and that the supercooled water freezes. In this process the 
latent heat of fusion (freezing) is lost to the surroundings and 
finally ice and surroundings return to the original temperature. 
Calculate AU, AH, and AS for both supercooled water and 
surroundings. 


Two moles of water equal 36g. If we neglect volume and 
pressure variations during phase changes between water and 
ice it follows from the first law and the definition of enthalpy 
that for both supercooled water and surroundings 


AH = AU =Q. 


There is a correct and a wrong way to calculate Q and AS 
during this process. If we calculate the amount of heat given 
away from the supercooled water as it freezes at —10°C, we 
find that Qw = —ml;(—10°C) = —0.036 x 0.312 x 10°J = 
—11232J. Then, if we use equation (5.7) we will find that 
AS, = f 6Q/T = —42.7J K~'. The amount of heat released 
is gained by the surroundings. Thus, Qs = 11 232 and AS; = 
42.7 JK~*. It follows that the change in entropy of the whole 
system (surroundings + supercooled water) is zero. But there 
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was a change in the system, which, according to the sec- 
ond law, should increase the entropy of the system. The 
reason for this violation is that we have indirectly assumed 
that the spontaneous freezing is a reversible process, which 
is not. 


A correct way would be the following: Start with the super- 
cooled water at —10°C and warm it reversibly to 0°C (step 
1), freeze it reversibly at 0°C (step 2) and then cool it 
reversibly to —10°C (step 3). In this scenario we have 


Step | 
Qi = MepwAT = 1518.5 J 
= dT 
AS, = MCpw— = 5.67J K! 
7 T 
Step 2 
Qz = —ml;(0°C) = —12013.2 J 
AS, = Q2 = —44.0 JK 
To 
Step 3 
Qs = mci AT = —758.2 J 
Ty T 
AS3 = eae = —2.83 JK" 
Tə T 


In the above calculations we have considered that T} = 
263K. Ty = 273K, Cpw = cw = 4218Jkg™ K}, cp = 
¢ = 2106Jkg™ K-t, and l(0°C) = 0.3337 x 10° Jkg™! 
K-t. It follows that the amount of heat released by the 
supercooled water is 


Qw = Q1 + Q2 + Q3 = —11 252.9 J, 
and the change in entropy of the supercooled water is 
AS, = AS, + AS; + AS; = —41.16JK™. 


The surroundings (Ts = —10°C) gain the amount that the 
freezing of supercooled water releases. Thus, 








AU; = AH; = Qs = —Qy = 11 252.9 J. 


At the same time their change of entropy is 


Qs Qs 5 
= —— = — = . . 
ASs Te = Fe = 42 795K 
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It follows that the change in the entropy of the whole system 
(surroundings + supercooled water) is 


AS = AS; + AS, = 1.63JK™! > 0. 


Problems 


(6.1) 


(6.2) 


(6.3) 


(6.4) 


(6.5) 


(6.6) 


(6.7) 


Prove that for any system 


(a), =7 (az) 

eS Se eel 

ƏV jy oT) y 

(Hint: express dU using equation (2.1) and substitute into 
the first law. Then consider the definition of dS = ôQ/T and 
the fact that dS is an exact differential, perform the partial 
differentiations and collect terms.) 

Using the equation derived in problem 6.1, show that for an 
ideal gas U is only a function of temperature and does not 
depend on V (Joule’s law). 

By plotting esw against T show that in the range —20°C to 
30°C equations (6.16) and (6.17) are virtually identical. If 
we were very picky about choosing one over the other which 
would you choose? 

Inside an open container is 500g of ice with temperature 
T = —20°C. Heat is then provided to the ice at a rate of 
100calmin™* for 700 min. Plot (a) the relationship between 
temperature (in °C) and time (in seconds), (b) the relation- 
ship between the heat absorbed by the ice and time, (c) 
the relationship between the heat absorbed and tempera- 
ture. Neglect the heat capacity of the container and assume 
that ls = constant = 79.7calg-! and c = constant = 
0.503 calg-1 K+. 

A closed box of volume 2m° and temperature 120°C con- 
tains 1 m? of water saturated with vapor. Calculate (a) the 
mass of water vapor, (b) the mass of water vapor that must 
slowly escape in order for the temperature inside the box to 
decrease to 100°C while maintaining saturated conditions. 
(1.43 kg, 0.71 kg) 

An open container of volume 2 liters contains dry air and 
some quantity of water at standard conditions. After equi- 
librium is achieved the container is closed and heated to 
a temperature of 100°C. What is the least amount of 
water that must have been added if after the heating the 
air is saturated? What is the final pressure inside the box? 
Use equation (6.16). (1.16 g, 2.35 atmospheres) 

A substance obeys the equation of state pV1? = 10°71. 
A measurement of its internal capacity inside a container 


3 


Figure 6.4 
Illustration for problem 
6.10. 


PROBLEMS 97 


(6.8) 


(6.9) 


(6.10 


m= 


(6.11) 








| Ice | 


a a a 





having a constant volume of 100 liters shows that under 
these conditions the thermal capacity is constant and equal 
to 0.1calK~'. Express the energy and the entropy of the 
system as function of T and V. (Hint: Is the substance an 
ideal gas?) 

How would you have to alter your assumptions in example 
6.2 in order to make cloud droplets grow to raindrops by 
condensation only? 

Calculate the change in entropy when 5kg of water at 10°C 
is raised to 100°C and then converted into steam at that 
temperature. At 100°C the latent heat of vaporization of 
water is 2.253 x 10°Jkg~*. In the above range of tempera- 
tures, the specific heat capacity under constant volume for 
water is nearly independent of temperature with a value of 
about 4.18 x 10? Jkg~* K1. (> 35972 J K~*) 

A block of ice rests as shown in Figure 6.4. A wire on which 
we attach a weight surrounds the ice. Explain why the wire 
will go through the ice block without the ice breaking into 
two pieces. 

The warm source of a Carnot cycle has a temperature of 
100°C while the cold source consists of melting ice. When 
the cycle operates for 1h it is observed that 1 ton of ice 
has melted. Calculate (1) the heat absorbed in the cycle 
from the warm source (2) the heat given away in the cycle 
to the cold source. The latent heat of fusion at 100°C is 
0.3337 10° Ike *. (333.7 x 10° J, 244.2 x 10° J) 


CHAPTER SEVEN 





Moist air 


Our atmosphere is basically a two-component system. One compo- 
nent is dry air and the other is water existing in vapor and possibly 
one of the condensed phases (liquid water or ice). According to Dal- 
ton’s law, in a mixture of ideal gases each gas can be assumed to 
behave as if the other gases were absent. As such, in a mixture of dry 
air, water vapor, and a condensed phase, the “water” system (water 
vapor + condensed phase) can be treated as being independent 
of the dry air. In this case, the concepts developed in the previ- 
ous chapter (for the one-component heterogeneous system “water” ) 
are valid for the two-component heterogeneous system “dry air + 
water”. We will call the system consisting of dry air and water 
vapor “moist air” and it can be unsaturated or saturated with water 
vapor. Since liquid water is absent, moist air is a two-component 
system with one phase present. As such, according to equation (6.1) 
we need three state variables to specify the system’s state. Usually 
these variables are taken to be pressure, temperature, and a new 
variable (to be defined soon) called mizing ratio. If the condensed 
phase is present and in equilibrium with the vapor phase, then two 
variables are needed (typically temperature and pressure). 

For clarity in our notation we will use the subscripts d, w, v to indi- 
cate dry air, liquid water, and water vapor, respectively. The only 
exception will be for the vapor pressure, which we will denote sim- 
ply as e. Variables with no subscripts will correspond to a mixture 
of dry air and water vapor. 

The following table gives various systems, their type, and the 
variables needed to specify their state (see equation 6.1). 
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7 MOIST AIR 
System Type Variables needed to 
specify system’s state 
Dry air and Heterogeneous— 3 (p, T mixing ratio) 
water vapor two components, 
one phase 
Dry air, water vapor, Heterogeneous — 
and liquid water two components, 
two phases 2 (p, T) 
Dry air, water vapor, Heterogeneous — 
liquid water, and ice two components, 
three phases 1 (T) 
Liquid water and Heterogeneous — 
water vapor one component, 
two phases 1(T) 
Liquid water, water Heterogeneous — 
vapor, and ice one component, 
three phases 0 
Water vapor Homogeneous — 
one component, 
one phase 2 (p, T) 
Dry air Homogeneous — 
one component, 
one phase 2 (p, T) 








7.1 Measures and description of moist air 


7.1.1 Humidity variables 


In a sample of moist air, dry air and water vapor have the same 
temperature T and occupy the same volume V. Thus, for water 


vapor 


or 


eV =m R,T 


e=p,R,T 


where e is the vapor pressure, py is the density of water vapor, and 
R, is the specific gas constant for water vapor. Since the molecu- 
lar weight of water M, = 18.01 gmol7’, the value of R, is equal 
to R*/M, = 461.5Jke~*K~-?. Since R* = RaMa it follows that 


RaM, = R,M,. Thus, 


7.1 MEASURES AND DESCRIPTION OF MOIST AIR 101 


Ra M, 
= — = — = 0.622. 7.1 
ee, (7.1) 
We define the specific humidity, q, and mixing ratio, w, as 
py M 
q = — = — 
p m 
and es 
yi ass (7.2) 
Pa Ma 


where m = mg+m, is the total mass of the mixture (py is also called 
the absolute humidity). Since p, = eM,/R*T, pa = paMa/R*T, 
and pg = p — e it follows that 








e 
w=e ; 7.3 
5 (7.3) 
At saturation the saturation mixing ratio is 
es 
Ws = € 7.4 
=e (7.4) 


where e, is either the equilibrium (saturation) vapor pressure with 
respect to liquid water (esw) or with respect to ice (esi). In general 
€s,€ < p which reduces the above two equations to 


ween, w YES. (7.5) 
Pp p 
Since 1/q = (ma +my,)/my = (ma/m,) + 1 it follows that 
1 1 
-= — + 1 
q w 
or 4 a 
= =—_., 7.6 
tira Eie (7.6) 


In the atmosphere both w and q are very small (w, q < 1). For this 
reason we can always assume that w 7 q. 
The relative humidity, r, is defined as 


My 





T= 
Mys 
where m, is the mass of water vapor in a sample of moist air of 
volume V and m,, is the mass of water vapor the sample would 
have had if it were saturated. Because of the ideal gas law we can 
write r as 
r= 


e 
= (7.7) 
Using (7.5) it follows that approximately 


w 
re. 


Ws 
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7.1.2 Mean molecular weight of moist air and other quantities 


According to equation (3.13) we have that the mean molecular 
weight of moist air is 


= Ma F My 
M = Md 4 my * 
Ma Mv 





We can manipulate the above equation by writing it as 


1 Ma | My 1 
M Ma M, Ma + My 
































or 
Le sa 1 q Ma 
M ~ Ma Ma T My KE M, 
or 
1 1 Ma q q Ty /Ma 
M Ma Ma + My M,/Ma 
or 
1 1 1 G+”) 
M ~ Mal+w E 
or 
1 1 1 q 
M Ma 1l+w E 
Since 
l _q 
ltw w 
and 
nn 
1-q’ 
then 
1 1 (1 4 ) 
M Ma q € 
or 
1 1 1 
aaa 
or i i 
= = — (1 + 0.61g). 
77 M, | q) 
Then for moist air the equation of state is 
R* R* 








= Roca T e Ss 1 +0.61q)T 
pa moist M a, gg q)T, 


d 
= Ra(1 + 0.619)T. (7.8) 


The above equation defines the virtual temperature, 


Tan = (1 + 0.61q)T, (7.9) 
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which can be interpreted as the temperature of dry air having the 
same values of p and a as the moist air. Otherwise stated, virtual 
temperature is the temperature that air of a given pressure and 
volume (or density) would have if the air were completely free of 
water vapor. Since q in reality is always greater than zero, Tūirt 
is always greater than T. The above equation also gives the gas 
constant of the mixture: 


R = (1 + 0.61q) Ra. (7.10) 


Similarly, we can define the specific heat capacities of moist air. 
Suppose that a sample of such air receives, at a constant pressure, 
an amount of heat 6Q which increases its temperature by dT. Some 
of this amount (Qa) is absorbed by the dry air and some (6Q,) by 
the water vapor. In this case we can write that 








ÒQ = 6Qa + 6Qy 
or 
ÔQ = madga + My dqv 
or py in 
ôq = 3 S x 
a ma + m, da + Ma nee 
or 


ôq = (1 — q)ôqa + qoqv. 


Recalling the definition of specific heat capacities we then have 


o= gat + 
or 

Cp = (1 — q)cpa + q Cpv 
or 

pE 
or 
Cp = Cpa (1 + 0.87q) ~ Cpa (1 + 0.87w). (7.11) 

Similarly 


cy = (1 — q)cva + qevy 
which leads to 
cy = cval 1 + 0.97q) © cva(1 + 0.97w). (7.12) 
From equations (7.11) and (7.12) we also have 


= Cp _ Cpa (1 +0.87q) 
T Cy o Cva (1 T 0.97q)’ 
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which because q < 1 can be approximated with 
y & yall + 0.87q)(1 — 0.97q) 
or after ignoring the second-order term 


7 Yall — 0.19), (7.13) 
where ya = 1.4. Similarly, 
y—1 R  Ra(1+0.61q) 
y Cp Cpa (1 + 0.879) 
x ka(1 + 0.61q)(1 — 0.87q) 
= ka(1 — 0.26q), (7.14) 





k= 





where ka = 0.286. Note that similar expressions in the presence of 
liquid water cannot be derived because liquid water + vapor does 
not constitute an ideal gas. 


7.2 Processes in the atmosphere 


In this section we will discuss scenarios that are relevant to basic 
processes in the atmosphere. To make the discussion flow effec- 
tively we will start with processes happening at a constant pressure 
level (for example, at the surface). After that we will proceed with 
processes involving ascent, first of unsaturated air and then of sat- 
urated air. Finally, we will discuss processes that involve mixing of 
two air masses in the horizontal and in the vertical. The laws and 
formulas describing each process will be derived analytically. Subse- 
quently, each process will be described visually in a (p, T) diagram. 
As we will see later, better diagrams can be designed. However, for 
now the simple and familiar (p,T) diagram will be used for illus- 
trating the various processes. In the following presentation when a 
condensed phase is present, we will assume that it is liquid water. 
The same theory holds when the condensed phase is ice as long as 
the appropriate substitutions are made (for example esi for esw, ls 
for ly, c for cw, etc.). 


7.2.1 Isobaric cooling — dew and frost temperatures 


Consider a parcel of unsaturated moist air and assume it is a 
closed system. In this case q and w remain constant. If this parcel 
begins to cool at a constant pressure the vapor pressure e will also 
remain constant as long as condensation does not occur (recall that 
e = wp/(w + €). However, the equilibrium vapor pressure will not 
remain constant. As we saw in the last chapter both esw and esi are 
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strongly dependent on T, and as T decreases so does es. Since es 
decreases it follows that the relative humidity increases and satura- 
tion is approached. The temperature at which saturation is reached 
is called the dew point temperature, Taew, if saturation is reached 
with respect to liquid water or the frost point temperature, T;, if 
saturation is reached with respect to ice. 

These temperatures derive their name from the fact that if the par- 
cel cools isobarically below these temperatures droplets (or ice) will 
form. However, keep in mind that for condensation (or sublima- 
tion) to happen at temperatures just below saturation point either 
condensation nuclei or solid surfaces must be present. Spontaneous 
condensation in the absence of nuclei or surfaces needs extremely 
high relative humidities — supersaturation — which cannot easily 
occur in the atmosphere. 

Isobaric cooling can happen naturally because of radiation cooling 
of the surface or the air itself. In this case the ground can lose 
enough heat to cool below the temperature at which atmospheric 
water vapor condenses. If this happens at temperatures greater than 
0°C water vapor in the layer very close to the surface condenses 
on the surface and becomes visible as dew. If it happens at tem- 
peratures less than 0°C the product is frost. If the cooling is very 
strong or if the air itself cools by radiation, the layers of air above 
the surface may also cool below their saturation point leading to 
the formation of radiation fog. Alternatively, isobaric cooling can be 
achieved when air moves in the horizontal (p ~ constant) and trav- 
els over colder regions. Then saturation may occur as the air cools 
by heat conduction to the surface. This leads to advection fogs. 
At Tuew vapor pressure becomes the equilibrium (saturation) vapor 
pressure (and to a very good approximation the mixing ratio 
becomes the saturation mixing ratio). As such if e is the partial 
vapor pressure at T, the dew point must satisfy the equation 


€sw(Taew) = €. (7.15) 


This is shown graphically in Figure 7.1 where point B is assumed 
to have a temperature T and vapor pressure e. To approach 
saturation isobarically we must draw a line from B perpendicu- 
lar to the pressure axis. Point C corresponds to a temperature 
where the water vapor pressure is its equilibrium vapor pressure. 
This point thus defines Ty... In this diagram AT = e,,,(T) and 
BT = CTiew = € = €sw(Taew). Obviously, since the equilibrium 
vapor pressure curve is known the relative humidity and dew point 
temperature of any point in the diagram can easily be found 
(r = TB/TA, Taw = OT — BC). If only the temperature is 
known then we need to know T).,, in order to estimate the rela- 
tive humidity and the absolute humidity. From the figure it follows 
that neither relative humidity nor the dew point spread (or dew 
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Figure 7.1 

Graphical definition of the 
dew point temperature, 
Taew- Clearly, 

€sw(Taew) = e(T). 


Figure 7.2 

Depending on where a 
given state is in a (p, T) 
diagram, Taew or T's may 
or may not be defined. 
For example, for state P 
only Taew can be defined. 
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0 Tis T T (°C) 


point depression), T — T4ew, are direct measures of the actual water 
vapor amount. Both describe how far we are from the equilibrium 
curve but similar values for r or T — Ty. correspond to different e 
and thus to different p,(absolute humidity). Note that, depending 
on where a given point in the (e, T) diagram is, Tyew (or T;) may or 
may not be defined. For example, in Figure 7.2, given point P only 
Tew Can be defined. For point P’, Tiew = Tg = triple point tem- 
perature. For point P” and in the presence of supercooled water, 
both T; and Ty, can be defined with Tj’ > Tikw- If supercooled 
water is not present, 77, cannot be defined. 

It can be shown that the dew point temperature of a rising parcel 
decreases. If we denote as Tyew,o and Tuew,z the dew point temper- 
atures at the surface and at some higher level z, respectively, we 
have that 

esw(Taew,0) = €9 = Ry pyTo 


Caw (Ldew,z) = €; >= Ry pyT,. 


As a parcel rises and expands, both its temperature and vapor 
density decrease (because the vapor is distributed over a larger 


p (mb) 











Y 
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volume). Thus eọ > e, and esw(Taew,o) > Csw(Tuew,z). Then from 
equation (6.17), it follows that Taew,z < Tuew,o- 


Estimating Taew 


Since e = @sy(Tuew) the C-C equation can be written as 








de — he 
dtc Beate 
or 
desw z l,dT 
ew RTZ 


Considering 1, œ~ constant and integrating from T4ew (where 
esw(Taew) = €) to T (where esw(T) = esw) we obtain 


Esw ie T= Tuew 
In = 
e Ry TTiew 
l = bs T= Tiew 
Se, Ry, TTaew 


T — Taew = —R,TTaew lnr/l. (7.16) 





or 


or 


Equation (7.16) provides Ty. from temperature and relative 
humidity. Using l, © 2.501x 10° Jkg~? and R, = 461.51 J kg™' K7}, 
the above equation becomes 


T — Taew = —1.845 x 10~4TTaew In r. 
Following the same procedure we can obtain (for ls = 2.8345 J kg~*) 
T — T; = —1.628 x 10-*TT; In r. 


Estimating the dew point temperature can be visualized in the fol- 
lowing (p, T) diagram (Figure 7.3) where pressure decreases along 
the y-axis. 

In this diagram a couple of dry adiabatic (solid) and a couple of 
saturation mixing ratio (broken) lines are shown. The dry adiabats 
can be calculated from Poisson’s equation 


In 0 = In T + kg In 1000 — kg In p 


by setting In 0 = constant: 


1 
In p= — ln T + constant. 
ka 
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Figure 7.3 

Illustration of the 
estimation procedure of 
the dew point 
temperature. 
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p (mb) 





1000 


Tis 02 0, > 05 T 


In a (p,T) coordinate system, this is not a straight line. The sat- 
uration mixing ratio lines can be calculated from equations (7.5) 
and (6.17) 








€€sw (I 
yw ol) 
Wsw 
a 6.11 5417 
pa = E exp(19.83 — ) 
sw T 
= 6.11 5417 
In p ~ In =E + 19.83 
Wsw T 


or for some Wsw = constant 
B 
hpx A+ — 
A T 


where A = In 6.1le/Wsw + 19.83 and B = —5417. In a (p, T) dia- 
gram this is not a straight line either. For this reason both the 
dry adiabats and the saturation mixing ratio lines are depicted as 
nonlinear curves. Since the purpose here is to illustrate the process 
graphically, these curves are realistic but by no means exact. Let 
us now assume that an air parcel is at a pressure level p and has 
a temperature T. Let us further assume that its mixing ratio is 
10gkg~*. The parcel is not saturated because at this (p, T) state 
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the saturation mixing ratio is 15gkg~’. In order to estimate the 
dew point, we then have to move, as the arrow shows, to the left 
(lower temperatures) (but remain at the same pressure level) until 
we “meet” the 10 gkg™* saturation mixing ratio line. At that tem- 
perature, the mixing ratio of the parcel is equal to the saturation 
mixing ratio and thus the parcel becomes saturated. That temper- 
ature is then the parcel’s dew point temperature. Note that if from 
(p, T) we follow the dry adiabat 6, down to 1000 mb we obtain the 
potential temperature of the parcel, which will be greater than its 
temperature, T. 


7.2.2 Adiabatic isobaric processes — wet-bulb temperatures 


Let us now suppose that our system is unsaturated moist air and 
liquid water and that the mixture has a temperature T and a com- 
bined water mixing ratio w, = (my + Mmw)/Ma. Since this is an 
unstable situation the system will tend to approach equilibrium. 
This will be achieved by water evaporating. If we assume that the 
system is closed and the process is taking place adiabatically (i.e. 
the process is also isenthalpic), then the liquid water must draw 
heat for the evaporation from the unsaturated moist air. Since no 
other processes are taking place, this decrease in temperature must 
be related to the change in the mixing ratio of moist air due to 
evaporating liquid water. Note that, since at any time during this 
process the system differs slightly from saturation, the process is 
spontaneous and irreversible. If there is enough liquid water, sat- 
uration will be reached. The temperature at saturation is called 
the isobaric wet-bulb temperature, Tw. An important atmospheric 
process described by this type of transformation is cooling by evap- 
oration of rain at a given level. The system consists of a certain air 
mass plus water that evaporates into it from the rain falling through 
it. If instead of water we had ice present, then we would define the 
isobaric ice-bulb temperature, Ti. Note that Tw is not the same as 
Taew Since the two cooling processes are different. Now, if we follow 
the reverse procedure, we could define the temperature which air 
would achieve if all its vapor were to condense in an adiabatic iso- 
baric process and with the produced water being withdrawn. This 
limiting temperature is called the isobaric equivalent temperature, 
Ta. Note, however, that the definition of this temperature is not a 
“legal” one. The reason is that we assume that the reverse of an 
irreversible process can happen, which is not possible. 


Derivation of the wet-bulb temperature 


Consider a closed heterogeneous system consisting of dry air, water 
vapor, and liquid water. If this system undergoes a change, then 
the change in enthalpy is 
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AH = H; — H; 


where the subscripts 7 and f refer to the initial and final states, 
respectively. We can write the above equation in terms of the three 
components 


AH = (Haş + Hss + Hwg) — (Hai + Hvi + Hwi) 
or 


AH n (mahas + My phy f F Mwy trys) g (Mahai + Myihyi oe Mwihwi). 
(7.17) 


If we recall the definition of l, (ly = h, — hw) and consider that the 
total amount of water Mm, = my + Mw is conserved (i.e. Myvi + Mwi = 
Myf + MwF = m), then the above equation becomes 


AH = Malhas = hai) + M (hwf = hwi) + My fly f = Myilyi- (7.18) 


If we now assume that in the range of meteorological interest Cpa 
and c,, are nearly independent of temperature, then we can write 
(recall equation (4.18)) that 


hag — hai = Cpa(T =) 
hf = hwi = Cw (Ty = T;). 


In addition, if T, is not greatly different from T; we can assume 
that 


lyi aaa lys = ly. 
Then equation (7.18) becomes 


AH = MaCpa(Ty a T;) + MeCw (Ty oo T;) + (my — Myi)ly 





or 
mM Myf — Myvi 
S a ae ae 
or 
AH/Ma = (Cpa + Weew) (T; — Ti) + (wy — wily 
or 
AH/ma = (Cpa + Wcw) AT + Aw. (7.19) 





Note that w, (the total water mixing ratio) is a different con- 
stant for each system and varies only with the concentration of 
liquid water. As such it does not depend on temperature. Since we 
have also assumed that Cpa; €w; ly # f(T) we can write the above 
equation in differential form: 


dH /ma = (Cpa + Wrew)dT + lydw. (7.20) 
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Equation (7.20) is a general expression for the change in enthalpy of 
a closed heterogeneous system of dry air, water vapor, and liquid 
water during any process. If we restrict ourselves to isenthalpic 
processes only, then equation (7.20) becomes 


(Cpa + Wcw) (T; = T;) = (wy = wily. (7.21) 


Equation (7.21) indicates that in an adiabatic isobaric process, 
where moist air cools from evaporating liquid water, the change 
in temperature relates to the change in the mixing ratio. It also 
tells us that T; depends on w, with Ty becoming minimum when 
w, > 0. Since m, = my + Mw cannot be zero (if it were neither 
the mixing ratio nor temperature would change), this condition 
requires that ma — oo. Of course this cannot be true but to a good 
approximation mg >> m,. We define this limiting temperature as 
the wet-bulb temperature, Tw, which is determined by setting in 
equation (7.21) T; = T, T; = Tw, W; = W, Wf = Wsw, and w, = 0: 


Cpa(T ~~ Fy) = (Wsw _ wl, 


or 
ly ly 
Ty + —UWew = T + >w (7.22) 
Cpa Cpa 
where Wsw is the saturation mixing ratio over liquid water at 
temperature T,,. As expected, since Wsw > w, Ty < T. 


Estimation of the isobaric equivalent temperature 


In order to estimate this temperature (Tyi), the only thing we have 
to do is to set in equation (7.21) w; = w, = w (the mixing ratio of 
the moist air), ws = 0, T; = T, and T; = Ta. Then we obtain 
pepe A = (7.23) 
Cpa + WCw 
Obviously Tei > T. Also, since for typical values of w (around 
10 gkg~*) ly /(Cpa + wcw) ~ 2500 it follows that Tea > Tyirt. Neither 
the wet bulb nor the isobaric equivalent temperatures is estimated 
graphically. They are calculated from equations (7.22) and (7.23), 
respectively. 


Relation between wet-bulb and dew point temperatures 


If we employ the general approximation w ~ ee/p we can write 
equation (7.21) as 


€ 
cp(T, — Ty) = a Sth 
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where cp = Cpa + Wcw. Setting T; = T, T; = Tw, and ef = egy (Ty) 
yields 
_ Per 


€i = esw(Tw) l E 


(T-T). 


From the definition of dew point temperature we have that e; = 
esw(Taew). Then, the above equation becomes 


pe 
ésa Las) = esw(Tw) = reas oe T): 


Since T — Ty > 0 it follows that esw(Taew) < esw(Tw). Then, from 
equation (6.17) it follows that 
Tiew < Tyi 


Are you keeping up with all these temperatures? Up to now we have 
discussed five different temperatures: T, temperature; Tyirt, virtual 
temperature; Tyew, dew point temperature; Tw, isobaric wet-bulb 
temperature; and Tai, isobaric equivalent temperature. They are 
related as follows: 


Tee XS RLT KD < Ter (7.24) 


Brace yourself. There are more! 


7.2.3 Adiabatic expansion (or compression) of unsaturated 
moist air 


Recall Poisson’s equation for an adiabatic expansion or compression 
of an ideal gas from p, T to p', T” 


IN k 
rool’ 
p 


which for p’ = 1000 mb defined the potential temperature, 0, 


6=T (=) f (7.25) 


Poisson’s equation is valid for any ideal gas. As such it is valid for 
moist air provided that the appropriate k is used and that con- 
densation does not take place. We can then define the potential 
temperature of unsaturated moist air as 


1 ka(1—0.26q) 
Om =T (=) , (7.26) 
p 


7.2 PROCESSES IN THE ATMOSPHERE 113 


If we substitute T from equation (7.25) applied to dry air into 
equation (7.26) we obtain (assuming that pa ~ p) 


—kqg0.26q 
0. = Os (=) 
p 


—0.07q 
On = 04 (=) , (7.27) 


or 


Pp 


As we have already discussed, in our atmosphere q < 1 which 
indicates that to a very good approximation 


(AN A Og. 


As in the case of virtual temperature, here again we can define the 
virtual potential temperature, Oy, by substituting Tire for T in 
equation (7.25): 


1 ka 
~) , (7.28) 


Ovirt PE Tyirt (= 
P 


This will be the temperature dry air will acquire if it were to expand 
or compress from a level (Tirts, p) to the 1000 mb level. Since Tyi; > 
T it follows that byin > 0. 

Note that even though q and w are very small (of the order of 
0.01) we should not jump to the conclusion that all variables or 
parameters of unsaturated moist air are for all practical purposes 
the same as those of dry air. If we consider a typical value of q = 
w = 0.01 we find that differences in virtual temperature of the 
order of 2-3°C are not uncommon. Such differences are often not 
negligible. On the other hand, as we saw above, differences in moist 
potential temperature are indeed negligible (of the order of 0.1°C). 
Consequently, from now on the symbol 0 will be used to indicate 
the potential temperature of both dry and unsaturated moist air. 
Similarly to the dry adiabatic lapse rate we may define here the 
moist (unsaturated) adiabatic lapse rate as 


g g 





= = 
Cp  Cpa(1 + 0.87w) 
or 
Pye ~ Iy(1 — 0.87w) (7.29) 
a TOS $ oy et 


7.2.4 Reaching saturation by adiabatic ascent 


As an unsaturated moist air parcel rises adiabatically, it cools. Dur- 
ing the ascent (and as long as no condensation occurs) its q or w 
remains the same but the equilibrium vapor pressure decreases. At 
the same time the vapor pressure e also decreases. This can be 
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easily seen by recalling that e = wp/(w + €). As w and e€ remain 
constant e/p remains constant. Since during the ascent p decreases, 
e must decrease if e/p must remain constant. Accordingly, when 
unsaturated moist air ascends adiabatically, its relative humidity 
(r = e/e,) increases only if e decreases at a slower rate than es. 
If this condition is not met clouds will not form in ascent but they 
could still form in descent. However, in our experience, clouds never 
do form as a result of descending air, so this condition is apparently 
satisfied. But what is the mathematics and physics behind it? 

If we log-differentiate equation (7.7) we get 


dlnr=dlne—d In ew. (7.30) 
For moist air, Poisson’s equation tells us that Tp 7 = constant. 


Since e/p remains constant, it follows that Tex = constant. Then, 
again by log-differentiating, we obtain 





4 
du T= dhe. (7.31) 
J 
Combining (7.30) and (7.31) with the C-C equation yields 
dl Yipo ap (7.32) 
n r= —— dln T — ; ; 
y—1 R,T? 


The first term on the right-hand side of the above equation accounts 
for the change in relative humidity due to the decrease in e while 
the second accounts for the change in relative humidity due to 
the temperature decrease. As one term is positive and the other 
negative, the net result can be either positive or negative. This 
indicates that adiabatic expansion can lead to an increase or to a 
decrease in r, which in turn implies that the opposite (i.e. adiabatic 
compression) can also lead to a decrease or an increase in r. This 
opens the road to bizarre scenarios where clouds can actually form 
on descent. If we rewrite equation (7.32) as 


de č r y lg 
dr T\g-1 RT)’ 


S ag 

y-1° RT 
then dr < 0 when dT < 0 and dr > 0 when dT > 0, which would 
indicate that clouds can form on descent. Assuming the typical 
values for y ~ 1.4, l, = 2.5 x 10ê J kg™*, and R, = 461.5 Jkg™* K7! 
we find that inequality (7.33) is satisfied only when T > 1550K. 
Obviously, such a situation does not exist on our planet and clouds 
can only form on ascent. However, all it takes is an exotic planet 
whose atmosphere has a condensable gas with smaller l, and all 
sorts of things may happen. 





we see that if 





(7.33) 
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Now that we have settled this issue, we can expect that ascend- 
ing moist unsaturated air will eventually reach a relative humidity 
of 100%. Immediately after that condensation will take place in 
order to sustain equilibrium conditions. The level at which satu- 
ration of ascending unsaturated moist parcels is achieved is called 
the lifting condensation level (LCL). The temperature at this level, 
Tic, is called the saturation temperature. By integrating equa- 
tion (7.32) from an initial state T, T4ew,r to the final state Tror = 
Taew,LCL> TLCOL = 1 we obtain 


Y Tic ly 1 1 
Inr = l . 7.34 
Tal T tile T oe) 








Since r = e/e,(T) and e = esw(Taew) we have from the C-C 
equation 





ly ly 
€sw(Taew) = 6.11 exp (sa ET ) 
v v £ dew 


and 








L ly 
esw(T) = 6.11 exp Ga ae 7) l 


Dividing the above two equations and taking logs yields 
ly T— Tiew 
ln r= 
Ry T Tew 
Accordingly, equation (7.34) becomes 


T a Tiew Trc 1 1 
———— =AIl . 
TT soy DE i (= 7) oa) 


aoon 


Equation (7.35) can be solved numerically to provide the saturation 
temperature Tucu from T and Ty... If you would rather not solve 
this equation you may use the approximation provided by Bolton 
(1980) 








where 


1 
+55, (7.36) 


Trcr = 1 lnr 
T—55 2840 








where T is in K, r = w/Wew, Wsw = €[Csw/(P — Csw)], esw = 
6.11 exp(19.83 — 5417/T), and T — Taew = —RvTTacw In r/ly. 

A parcel that rises from an initial level where the temperature is T 
to a level where the temperature is Tyo, has expanded. Therefore, 
the vapor density at LCL is smaller than that at the initial level. 
From the ideal gas law it then follows that e(Tici) < e(T). Because 
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Trot = Taew,ucL and e(T) = esw(Taew) this inequality is equivalent 
to €sw(Tict) < €sw(Taew). Then from equation (6.17) it follows that 
Tror < Taew. Thus, relation (7.24) can be extended to include yet 
one more temperature: 


Trcr < hiv < Tw <T < T, < Ta. (7.37) 


Estimating the height of LCL 


The height of LCL for a given parcel (zcu) depends only on its 
temperature and relative humidity. At zo, the parcel’s tempera- 
ture equals its dew point: Tuc, = Taew,.ct- AS the parcel is lifted 
from a reference level z its temperature decreases according to 


T(z) = To — Taz — zo) (7.38) 


where T is the temperature of the parcel of zọ and we have assumed 
that I, œ% TIa. Similarly, its dew point temperature decreases 
according to 

Tew (Z) = Taew,0 T Taew(2 gi 2). (7.39) 


Here both I4 and Iyew are considered constant. From equations 
(7.38) and (7.39) we have that 


Trcr = To — Ta(Zucr = Zo) 
Taew, LCL = Taew,0 wi Daew(ZLCL T zo) 
or since Trcr = Taew,LCL» 


To E Taew,0 
I ~ Loew 





ZLCL — 20 = 


(7.40) 


In the above equation we know everything except for yew. In order 
to estimate [yew we proceed as follows. The equation that defines 
dew point is e = esw(Taew). Differentiation with respect to z gives 


desw(Taew) = de 
dz dz 
or 


desw (Tiew ) dT jew de 
= —. .41 
dTaew dz dz ee) 





Now we know that e = wp/(w + €). If we assume that w remains 
constant until condensation occurs, we have 


de w dp 


dz wtedz 





or 


de edp — €sw(Taew) dp 
= = P A2 
dz pdz p dz (ee) 
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Combining equations (7.41) and (7.42) gives 
1 desw(Taew) dTaew _ 1 dp 


Esw (Tiew) dTaew dz p dz ` 





After considering the C-C equation we can reduce the above to 


lL, dTiew ldp 


RT? dz p dz (oaa) 





Since the ascent is assumed to be adiabatic, we have from Poisson’s 
equation that 


Tars ye 
y 


or 


dT y-1 dp 
= i 44 
Tdz y pdz (een 





From equations (7.43) and (7.44) it follows that 


Le AT sew = aA dT 
R,T? 


dew 





dz y— 1Tdz 


~ aaa Oe) 


or assuming y © Yq 





Ll, AT sew = g 
RT2w dz Ral 


which then gives 


dTaew pee g Ties 


d; ew — EE ‘ 
s dz el, T 





The value of the right-hand side in the atmosphere ranges from 
about 1.7 to 1.9°Ckm™. Therefore we can assume that on the 
average Taew decreases with height at the rate of 


Pie = 1.8°Ckm™*. 


By dropping the subscript 0 in equation (7.40) it follows that the 
height of LCL is given by the approximate relationship 


T— Tiew 
8 7 


Equations (7.35) or (7.36) and (7.45) give the temperature and 
height (in km) of the level where clouds begin to form from the 
temperature and dew point temperature at some reference level z. 
Reaching saturation via adiabatic ascent and the estimation of Trc 


(7.45) 


ZLCL T Z X 
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Figure 7.4 

Illustration of the 
estimation procedure of 
the lifting condensation 
level and saturation 
temperature. 
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and the height of LCL can be visualized in Figure 7.4. The setup in 
this figure is the same as in Figure 7.3. Since the parcel is initially 
unsaturated, as it ascents, it follows the dry adiabat passing from 
its initial state (p, T) (recall that 0m œ% ĝa). As it rises (follow the 
arrow) it cools and approaches the Wsw = 10gkg~* line, which it 
intersects at a pressure level (height) pict (Zucx). At that point, the 
temperature is Tyc,. As the figure illustrates, Tyo, < Taew- 


7.2.5 Saturated ascent 


Once saturation is reached, further ascent results in a further 
increase in relative humidity, which means that the vapor pressure 
becomes greater than the equilibrium vapor pressure. At this point, 
the system will return to equilibrium after the “extra” water vapor 
condenses and water droplets (or ice depending on the temperature) 
form on condensation nuclei. From this point on two extreme pos- 
sibilities may be considered. (1) The condensation products remain 
in the parcel. This is a reversible process because if we reverse 
the ascent the products will evaporate. It is also adiabatic as we 
assume that no heat is exchanged between the parcel and the envi- 
ronment. Since it is reversible and adiabatic it is also isentropic. (2) 
All condensation products fall out and the parcel consists always 
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of dry air plus saturated water vapor. This makes the parcel an 
open system. Obviously, this process is neither reversible nor adi- 
abatic and thus it is not isentropic. We will call this process a 
pseudo-adiabatic process. In reality these two extreme cases may 
never exist as some condensation products are likely to remain sus- 
pended in the parcel. However, because of the conditions imposed 
the extreme cases are easier to treat. Nevertheless, the study of the 
extreme cases offers useful insights and approximations for the real 
cases in between. 


Reversible saturated adiabatic process 


Since the process is assumed to be reversible, the heterogeneous 
system of dry air, water vapor, and liquid water (or ice if it is too 
cold) must be at equilibrium at all times. In this case the total 
entropy of the system, S, is conserved: 


S = Sa + Sy + Sw = constant. 


If we divide all the terms in the above equation by the mass of dry 
air we arrive at the following expression 


m. 
= Sq + Sy 4 Si 
Ma Ma Ma 





Recalling that the parcel is saturated (w = Wsw) and the definition 
of the total water mixing ratio (w, = (m, + Mw)/Ma), the above 
equation can be written as 


S 
— = Sa + WewSv + (W; — Wsw)Sw- (7.46) 
Ma 


For a reversible process 


jeez Oa Ove a (37) ar+ (5) av] + Zav: 
V T 





T T T T |\ oT OV 
(7.47) 
In addition, 
Os OSs 
= = ae . A 
dS (Se) ar+ (55) av (7.48) 


From equations (7.47) and (7.48) it follows that 


(ar), =r (ar), a 


(S). r (S): (7.50) 





and 
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For liquid water and water vapor in equilibrium at temperature T 
and pressure p equation (7.50) gives 

Se miSa. grate 
V, z Vw 7 V, Fs Vy 


T 





or 


T(S, — Sw) =U, — Uw + p(V, — Va). 
If we use equation (6.3) the above equation becomes 
T(S, — Sy) = H, — Hy = L, 


or 
T(sy — Sw) = ly. 


Substituting s, from the above equation into equation (7.46) results 
in 





= Sa + UWiSw 4 ; (7.51) 
Now recall that 
Sa = Cpa ln T — Ra ln pa + constant [p = pa + esw(T)] 
and, from equation (5.9), 
Sw = Cy In T + constant. 
Accordingly, equation (7.51) becomes 


S ly Sw 
— = (Cpa + Wicw) In T — Ra In pa + = 
Ma 





+ constant. (7.52) 


Since here S and mg are conserved, it follows that 


l 
(Cpa + wcw) In T — Ra ln pa + S 


= constant’. (7.53) 
As previously, here we can assume that w, Cpa, Cw # f(T). However, 
Wsw through equations (7.5) and (6.17) is a function of T. If we 
assume that l, is a function of T (and that may be a good idea 
here as the parcel may rise to very cold temperatures), then we can 
write equation (7.53) in differential form as follows: 


ly sw 
(Cpa + Weew)d In T — Rad In pat+d (=) =0. (7.54) 


Equations (7.53) and (7.54) describe reversible saturated adiabatic 
processes. Note again that here T and p correspond to saturation. 
If we now define 0’ as 


Ra/(cpatwtew) 
=T (=) (7.55) 
Pa 
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it follows that 
(Coa + WeCw)d In O’ = (Cpa + Wicw)d ln T — Rad ln pg. (7.56) 


Combining (7.56) and (7.54) gives 


(Cpa + WrCw)d In 0" = —d (=) 


or 





dln 0’ = al by Wsw | 


Cod + WiCw)T 
Integrating the above equation results in 


1 Wsw 


0’ exp | ——————_ | = constant. 7.57 
laren wee 
The above equation defines the family of curves that describe 
reversible saturated adiabatic processes. We can evaluate equation 
(7.57) at a state where all the vapor has condensed (w,,, = 0). This 
defines a new temperature, the equivalent potential temperature 6, 
given by 
Ly Wsw 


b = 6 ——— 
oe lo + wcw) T 


| = constant. (7.58) 


Although 0. is defined for saturated air it can also be defined for 
unsaturated air at temperature T, pressure p, and mixing ratio w 
by considering it to be lifted to the LCL. At LCL we have 


: 1000 Ra/(cpatweew) 
= Tuc (=) 
Pa 


(equation (7.55)) and thus 


exp l. (Tuc) w 
(Cpa + WeCw)Tiucr 








(7.59) 


d 


1000 Ra/(cpat+weew) 
be = Ttc ( ) 
Pp 


Note that in the above two expressions pa = Pror — €sw(Tict)- 
Another way to express 6, for unsaturated air at T, p, w is given by 
Emanuel (1994): 


Ra/(cpatwtew) 
6. =T (=) i p7 wRv/(cpatweew) 
Pa 


x exp |lyw/(cpa + Wiew)T] (7.60) 





where r is the relative humidity. Here py = p — e(T). When the 
initial mixing ratio is zero, both (7.59) and (7.60) give (as expected) 
e = 0, the potential temperature. 
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Recall that the moist potential temperature 0,,,(7.26) is 


1 Ra/cpa(1—0.26w) 
=F (=>) l 
Pp 


Comparison between 6,,(% 0) and 0’ indicates that 6’ is only 
approximately equal to 0. Given the definition of 6’, it follows from 
(7.60) that the equivalent potential temperature is approximately 
equal to the potential temperature a parcel would have if it were 
lifted to very low pressure so that all its water vapor condensed 
(w = 0). 


Pseudo-adiabatic processes 


The description of reversible saturated adiabatic processes depends 
on the value of w,. The problem with this dependence is that 
although Wsw can be determined at a given state, (p, T), w, can- 
not. The liquid water mixing ratio can take on any value and as 
such at a given (p,T) point there may exist an infinite number of 
reversible saturated adiabats. If we were to describe such a pro- 
cess in the (p, T) domain this would present a great inconvenience. 
This problem can be dealt with by defining a new process called 
a pseudo-adiabatic process. In such a process we do away with liq- 
uid water by assuming that it is removed as soon as it is produced. 
Obviously, our system is now an open system and this process is not 
reversible. However, we can think of the whole process as a two-step 
process. During the first step we have a reversible saturated adia- 
batic process, with condensation of a mass of water (dm,,). During 
the second step the water produced is leaving the system. For the 
first step the entropy may be defined as previously but with the 
liquid water term omitted. Accordingly, we set in equation (7.52) 
Wt = Wsw and we obtain 


ly Wsw 


T 





Ma 


D = (CpatWewCw) ln T—Ra In pet ( ) constant. (7.61) 


Since in this stage S and mg are conserved we can write the above 
equation in differential form 


ly Wsw 





A[(Cpa + WswCw) ln T] — Rad In pa + d ( ) =0. (7.62) 
After water forms it is removed immediately, thereby decreasing the 
entropy of the system. However, this process has no effect on the 
values of T and p. As such equation (7.62) describes the variation 
of p and T in a pseudo-adiabatic process. Equation (7.62) is very 
similar to equation (7.54) but here Wsw depends on temperature 
whereas w, in equation (7.54) does not. This makes equation (7.62) 
very difficult to solve analytically. However, numerical techniques 
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can be employed to define, similarly to the case of the equivalent 
potential temperature, a pseudo-equivalent potential temperature, 
Oep- As is shown in Bolton (1980), within 0.3°C, fep is given by 


1000 \ ®2850-0-28w) 3 
bep =T (=>) exp fea + 0.81w) ( ee 2.54) 
Pp Tror 





(7.63) 


where (T, p, w) is any state of the parcel saturated or not (w here 
is dimensionless and T, Tycy are in K). 

The pseudo-equivalent potential temperature can be interpreted 
as the actual temperature achieved by the parcel under the fol- 
lowing thermodynamic processes: (a) pseudo-adiabatic expansion 
to zero (or to a very low) pressure level at which we can assume 
that all water vapor has condensed and fallen out, and (b) sub- 
sequent dry adiabatic descent to 1000 mb. Note that according to 
the definition of 0. (the condensation products do not leave the 
parcel) no similar meaning can be attached to 0e. Nevertheless, 
the differences in the temperature of a parcel in a reversible sat- 
urated adiabatic and in a pseudo-adiabatic process are not great 
and are often neglected as processes in the atmosphere lie some- 
where in between true reversible saturated adiabatic processes and 
pseudo-adiabatic processes. This can be easily verified if we assume 
that Wsw © Ws K Cpa (ie. Cpa + WswCw © Cpa) and pa > ews (i.e. 
p = pa). In this case equations (7.54) and (7.62) become identical: 
Cpad In T — Rad In p + d (wswly/T) = 0. 

Graphically a saturated ascent can be described with the following 
diagram (Figure 7.5). At LCL the parcel is saturated. Now, as it 
continues to rise, it follows a path (depicted by the broken line) 
dictated by the saturated adiabat passing through (prc, Tucu) (in 
practice by the pseudo-adiabat, i.e. from equation (7.62) or (7.63)). 
At some very low pressure level (say, 200 mb), we may assume that 
all vapor has condensed and the condensation products have fallen 
out. At that limit, the pseudo-adiabat approaches a dry adiabat 
(solid line). This dry adiabat corresponds to ep (in other words if 
followed down to 1000 mb, it gives the pseudo-equivalent potential 
temperature). 


7.2.6 A few more temperatures 


The dry adiabatic descent from zero pressure level to the 1000 mb 
level defined the pseudo-equivalent potential temperature, Oep. At 
the original (p,T) level it defines the pseudo-equivalent tempera- 
ture, Tep (Figure 7.6). Since the descent is dry adiabatic we have 
that 
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Figure 7.5 

Illustration of the 
estimation procedure of 
the pseudo-equivalent 
potential temperature. 
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p (mb) 


200 






(Pict, Ta) 





1000 
0 Pep T 
p \ři 
Ta = b. (2) . 
p = Oev | T000 
Since T = 0 (p/1000)** it follows that 
bop 
Toet- (7.64) 


During the process by which Tep is defined, much more water vapor 
condenses than it would have if it condensed during an adiabatic 
isobaric process at an initial level below LCL. Therefore (recall 
the definition of Tai), if all that vapor were to condense at the 
initial level it would result in a temperature greater than Tai. Thus, 
To > Tei, which makes Tep the highest temperature. 

Many more temperatures can be defined but I will spare you this! 
Just keep in mind that at a given state, if the air is unsaturated, a 
dry adiabat can be followed to 1000 mb to define a corresponding 
potential temperature. As such from the state (Ti, p) the isobaric 
equivalent potential temperature, @,;, can be defined. Because Tep > 
Ta it follows that bep > bei. It also follows that bep > @. 

Similarly, if the air is saturated, a pseudo-adiabat can be fol- 
lowed to either the original pressure level or to the 1000 mb level 
to define yet more temperatures. For example, and recalling that 
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Figure 7.6 p (mb) 
Illustration of the 

estimation procedure of 

Top, Twp, and Owp- 200 






(Pictur Tic) 


1000 


wp ep 


Ty is attained at saturation in an adiabatic process, following 
the pseudo-adiabat corresponding to point (Tcu, pict) to (T,p) 
level defines the pseudo-wet-bulb temperature, Twp, and continuing 
down to 1000 mb defines the pseudo-wet-bulb potential tempera- 
ture wp (see Figure 7.6). Obviously, since both ep and wp are 
defined by the same curve they must be related and conserved 
in a pseudo-adiabatic process. Oy, can be estimated from equa- 
tion (7.63) by setting T = Tic, = Owp and w = w’ = u,(p = 
1000 mb, T = wp), ie. 





Bap = Ourp EXP bwa + 0.81w') (= z 2.54) . (7.65) 
wp 

Similarly, following the pseudo-adiabat passing through point 
(Tw,p) to 1000mb defines the isobaric wet-bulb potential 
temperature 0w. Here we should note that Twp and Tẹ (or Owp 
and 6,,) are not exactly the same. T, and 6, correspond to a 
true saturated adiabatic process whereas Twp and wp correspond 
to a pseudo-adiabatic process. Since water remains around in a 
true saturated adiabatic process it turns out that Twp < Tw and 
Dis Oss 
It is interesting to note here that if at the initial point (p, T) (see 
Figure 7.6), the amount of vapor content (w) is higher, then the dry 
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adiabatic ascent will last a shorter time and therefore the LCL will 
be found at higher pressure and temperature. The pseudo-adiabat 
passing through this point will be of a greater temperature than 
before. This means that higher wp (or Oep) values would correspond 
to higher moisture (in accordance with 7.63). 

Here is then the final breakdown of all the temperatures and 
potential temperatures and their relationships: 


Trons Lee < Tp = de OE ik = Ta <en 


(7.66) 
Ozi < bw < ox Om) < Ovirt < boi < Dep: 


7.2.7 Saturated adiabatic lapse rate 


Since condensation of water vapor releases heat during a reversible 
saturated adiabatic ascent, the cooling of the ascending air slows 
down. As such the corresponding lapse rate should be smaller than 
that of the dry or moist unsaturated adiabatic ascent. 

We could derive an expression for this saturated lapse rate, I}, by 
starting with equation (7.54) 


dT dpa L, Wew = 
Cop Ma +a( T )=0 





where Cp = Cpa + WiCw- Multiplying all terms of the above equation 
by T/cp gives 
R T 1 ly sw 

a” d(p — egy) + —Ad(lyWew) — aT =0, 


CpPa Cp Cp 





dT 


where we have considered that p = pg+e,,,. The second term, using 
the C-C equation in the form de, /dT = eswly/ R,T? and the fact 
that Raesw/Rvpa = Wsw, can be written as 


R T mE N 
gpi nar. 
CpPa Cpl’ 





Substituting this into the previous equation, dividing all terms 
by dz, assuming that dp œ~ dpa, and using the hydrostatic 
approximation, leads to 

dT g l d 


= L,Wew)- 7.67 
dz Cp Cp a Wsw) en 





Note that since in the above equation c, # Cpa, the term —g/c, is 
not exactly the dry adiabatic lapse rate. By performing the differ- 
entiation on the right-hand side we can write equation (7.67) as 
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aT g Wsw dl, ly dWsw 











dz Cp Cp AZ Cy dz 
or 
qT g Wsw dl, dT ly dWew 
dz Gs Cp dT dz Cy dz 
or 
T l 
E E ee (7.68) 
Cp dT) dz Cp Cp dz 


Using equation (6.12) the above equation becomes 


h Wsw (Cw _ 2] dT = g ly dWsw 


dz Cp Cp dz 





Cp 


Typical values for Cy,Cpv,Cp, and Wsw are 4220Jkg™*K~-}, 
1850 J kg™* K~!, 1050 J kg™* K-t, and 0.01, respectively. Such values 
make the term Wey (Cw —Cpv) /Cp of the order of 0.02. As such, it can 
be neglected to arrive at the following approximation 

dT g ly dWew 


= (7.69) 


dz Cp Cp dz 





Now all we have to do is to find an expression for dwsw/dz. Log- 
differentiating equation (7.5) leads to 


1 dWsw = 1 dp 1 desw 





Wsw dz pdz Caw dz 


1 dp 1 desw dT 
pdz esw dT dz 








where p is the total pressure of the parcel or the surroundings. If 
we again use the hydrostatic approximation and the C-C equation 
the above equation becomes 


l, dT 
RT? de 


1 dWsw g 


ww dz RT ' 





(7.70) 


Combining equations (7.70) and (7.69) leads to 


dT g 1l+lhWww/RT 
r. = = 7.71 
i dz Cy 1+ wey /CpRyT? en) 





where Cp = Cpa + WiCw and R is the gas constant for the surround- 
ings. In other thermodynamics books you may find some slight 
differences in the formula giving I. Most of the differences are due 
to certain assumptions and approximations. For typical values of 
saturation mixing ratio and temperature in the lower troposphere, 
T, ranges around 5°C km™+, which is half the dry adiabatic rate. 
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7.3 Other processes of interest 
7.3.1 Adiabatic isobaric mixing 


Let us consider two air masses of temperatures T} and Tọ and 
humidities qı and q2, which are mixing at the same pressure. This 
corresponds to a horizontal mixing. If we assume the mixing to be 
an adiabatic process, then it is also isenthalpic. If we assume that 
no condensation takes place, then we can write for the change in 
enthalpy of the whole system 


where 


Ah, = Cp: AT = cp (T = T) 
Aho = Cp AT = Cp(T 3 Tə) 


and T is the final temperature of the mixture. Considering equation 
(7.11), equation (7.72) becomes 


M1Cpa(1 + 0.87q1)(T a Tı) + M2Cpa(1 + 0.87q2)(T — Tə) =0. 
Solving for T we get 


= (7m Ty + m2T>) + 0.87 (miq T; + M2q2T2) , (7.73) 
m+0.87m, 





In the above expression m = mı + Mə is the total mass and 
My = M141 + M2q2 is the total water vapor mass. Since there is no 
condensation, m, remains constant. If the final specific humidity is 
q, then the water vapor mass is qm. It follows that 

+ 

ga C (7.74) 

m 
i.e. the final q is the weighted average of qı and q2. Using equation 
(7.74) we can rewrite equation (7.73) as 


= (mT, + mT) + 0.87(m191T) + M2q2T>) 





; 7.75 
m(1 + 0.87q) ee) 
If we neglect all vapor terms the above equation becomes 
pe mht mth (7.76) 
m 


indicating that, approximately, the temperature of the mixture is 
also the weighted average of the initial temperatures. One can 


Figure 7.7 

When two initially 
unsaturated air masses 
(Ay at rı and Ag at r2) 
mix isobarically, the final 
product (A) may be 
saturated. 
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show that the potential temperature and the vapor pressure of the 
mixture will also, approximately, be the weighted averages of the 
initial values (see problem 7.1): 


M10, + M202 


Ox (7.77) 


m 


Mı€1 + M2€2 


en 


(7.78) 


m 


Equations (7.76) and (7.78) indicate that in a vapor pressure- 
temperature diagram (Figure 7.7) if the two masses correspond 
to points A, and A, then the mixture will correspond to the point 
A which will lie on the straight line connecting A, and As. From 
the similarity of triangles AA A; and A2AjA it follows that point 
A will be found at a distance such that A,;A/AA2 = AA‘/A2A, = 
(Tı — T)/(L — T2) = m/m, (here equation (7.76) was used). The 
thin solid lines in Figure 7.3 represent the equal relative humidity 
curves corresponding to points A, A2, and A. Because of their cur- 
vature it can easily be seen that the final relative humidity will be 
greater than the weighted average. This can easily be visualized if 
we consider the simple case where the air masses have the same rel- 
ative humidity, i.e. points A, and A, lie on the same equal relative 
humidity curve (Figure 7.8). Then point A will lie to the left of 
the curve indicating that the final relative humidity is greater than 
the weighted average, (mr + mr)/2m = r. We conclude that iso- 
baric adiabatic mixing increases the relative humidity. This is a very 
interesting result because it implies that two air masses that are ini- 
tially unsaturated (but close to saturation) may result in a mixture 
that is supersaturated. In this case condensation will occur. This 
mechanism can produce mixing fogs and other important processes 
in the atmosphere. 


Fanar 


e 


e2 
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Figure 7.8 

When two saturated 

(r = 1) air masses (Ai 
and A.) mix isobarically 
the final product (A) is 
supersaturated. 
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7.3.2 Vertical mixing 


Let us now look at vertical mixing occurring because of turbulent 
and/or convective processes. Such a process is much more com- 
plicated than the previous one because now the vertical variation 
of T and p come into the picture. However, we can simplify the 
treatment of vertical mixing by (1) considering two isolated air 
masses M, and mz at pressure levels p, and pə (pı > p2) and tem- 
peratures T} and Tp, which are taken to a pressure level p where 
they mix isobarically, and (2) redistributing the mixture in the 
layer Ap = pı — p2. The excursions to level p can be considered as 
adiabatic ascent and descent (during which the specific humidities 
qı and qz remain constant) resulting in new temperatures for the 


air masses: 
k 
Pi 


k 
T: =T, (4) 
P2 


Once at level p the air masses mix isobarically. From what we 
discussed in the previous section it follows that in this case 


mT} + MTs 
m 
Mıqı + M242 


qa (7.79) 


m 


Ts 


and 





6 x mıı a a M24 


m 


where 6, and 62 are the initial potential temperature values, 


which remain constant during the adiabatic ascent and descent. 
The subsequent redistribution of the mixture in the layer Ap 
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will preserve 0 assuming that it consists of adiabatic ascents and 
descents. It follows that when the layer is well mixed, 0 will not 
vary with height. We can generalize equation (7.79) by considering 
mixing of n air masses and employing the general expression for 
the weighted average. For example, 


pee 
m 
_ Jo adm 
= 4 
m 
” wd 
pam = 
m 


Using the hydrostatic approximation we can express the above 
equations as 
Aime Sged 
qt = z = 
Jo pdz Pı — P2 





where zx is 0, or q, or w. 
The temperature variation with pressure can then be expressed 
according to 

raa Ey 7.80 

(iw) c ee) 

If in the initial distribution (before mixing) there is a level 
(p,T < T), where r = 7, then this level is defined as the mixing 
condensation level. 


7.3.3 Freezing inside a cloud 


As air rises and cools and clouds begin to form, it is often the case 
that the developing cloud would reach the freezing level. From that 
point on, freezing of cloud droplets (or rain drops) would occur 
and water vapor and ice will co-exist. Does this freezing affect the 
further development of the cloud? We can get an answer to this 
question by making the following simple assumptions: 


1) Initially the air is saturated with respect to water at 
temperature T. 

2) As water freezes latent heat is released into the environment. 

3) Once water freezes, water vapor is supersaturated with respect 
to ice. Depotition is then initiated releasing more latent heat. 
This continues until the vapor pressure becomes equal to the 
saturation vapor pressure over ice. 

4) The whole process occurs at a constant pressure level and it is 
adiabatic. Since it is adiabatic and isobaric, it is also isenthalpic. 
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Accordingly, we can divide this process into three “treatable” steps: 


Step 1: 


Step 2: 


Step 3: 


Water freezes at a constant temperature T. In this case we 
have a change of phase at constant pressure and the change 
in enthalpy is (see Chapter 6) 





AH, F —6@Qp=constant = —L; = —liMw 


where Mw is the amount of water that freezes. The minus sign 
is due to the fact that heat is released. Dividing the above 
equation by the mass of dry air yields 

AH, 


Ma 





= -hw (T) (7.81) 


Vapor condenses on ice at a constant temperature T until the 
water vapor pressure become the saturation vapor pressure 
over ice at a temperature T’ > T. Like in Step 1, the change 
in enthalpy during this step will be 


AH, = =l (Mw = mi) 
or 


AH, 


Ma 


= —I,[wy(T) — w;(T’)| (7.82) 


The whole system is warmed from T to T’. The change in 
enthalpy during this term will be (recall 4.18) 





AH; = C AT 
or 
AH; 
ma TAT (7.83) 


where it is assumed that C,/ma œ% Cpa. Since the whole 
process is isenthalpic, it follows that 


lew (T) = lafwa (T) — wi(T)] + cpa AT = 0 
or 


AT = lpww(T) + 1,[ww(T) = w(T’)| 





Cpa 
or using (7.5) 
celsesw(T) + €ls[€sw(L) — esi(T’)] 


PCpa 


AT = 
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In the above equation and in the range of typical conditions in 
our atmosphere, the right-hand side is positive. As such the effect 
of freezing is to increase the temperature inside the developing 
cloud. Typically, this temperature increase is of the order of 3K, a 
significant warming during cloud development. 


Examples 


(7.1) 


The initial state of a parcel of air is p = 1000 mb, T = 30°C, w = 
14 gkg"t. Find (a) T, (b) Tirts (c) Taews (d) Ts (e) Tis (f) Tic, 
(8) zxcrs (h) Peps (i) Toeps (I) wps (K) Avice: 


(a) 


0 UD! 
ae 


From equation (6.17) we estimate that esw(T) = 
43.1mb. From equation (7.3) we then find that 
e(T) ~ 22.0mb. Thus, r ¥ 0.51, Wsw © 28eke7". 


Tan = (1 +0.61q)T = f +0.61 (= )] T = 305.55K. 


l1+w 
From equation (7.16) we obtain Tyew = 291.75 K. 
From equation (7.22) we have 


Ty =T + eig — Wsw) 
Cpa 
where l, (T = 303K) = 2.43 x 10° Jkg~*. Now recall 
that wey in the above formula is the saturation mix- 
ing ratio at Tw. Making use of the approximations in 
equations (7.5) and (6.17) we have 





(eg E E T (1983 z | 
CpaP Tw 

Numerical solution of the above equation gives Tw ~% 
295.4K. 
From equation (7.23) for cw = 4218 Jkg~* we find that 
Ta ~ 335K. 
From equation (7.36) it follows that Trcr ~ 289K. 
From equation (7.45) and assuming that the initial 
level is z = 0 we find that zcr = 1.4km. 
Using equation (7.63) with T = 303K, p = 1000mb, 
w = 0.014gke"*, Tuor = 289K, we obtain that Oep = 
345K. 
From equation (7.64) it follows (since T = 0) that 
Top = Oop = 345K. 
Considering that w’ = w, (1000mb, T = 6,,) = 
(e/p)6.11 exp[19.83 — (5417/0wp)], numerical solution 
of equation (7.65) yields 0,,, + 295 K. 
From equation (7.28) and since p = 1000 mb it follows 
that Ovire = Tvis = 305.55 K. 
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(7.2) 


Note that all estimated temperatures and potential temper- 
atures are consistent with equation (7.66). 


The temperature of saturated air at 1000 mb is 10 °C as it begins 
to drop through radiation loss. As a result, some water vapor 
condenses and radiation fog forms. During the formation of fog, 
12x 10° J kg", of heat is lost to the surroundings. Find the final 
temperature and the decrease in vapor pressure. 


The process is non-adiabatic and thus non-isenthalpic. 
From equation (7.20) we have 


dh = c,dT + l,dw 


where Cp = Cpa + WtCw and we have assumed that dH/ma ~ 
dh. During the process the air remains in a state of satura- 
tion, thus w corresponds to saturation. Also, the process 
is isobaric, so dh = ôq and we can write the above 
equation as 


ôq = C,dT + ly dwew- (7.84) 


Using the approximate relationship Wsw œ% €€sw/p, we find 
that dwsw © (€/p)desy or, using the C-C equation, 
€ Ly esw 


dw = pal: (7.85) 





By substituting (7.82) into (7.81) we obtain 


2 
ôq = (o + ely Cw ) dT. 





pR,T? 


In the above equation, c, is to a good approximation 
independent of temperature. Also, for small temperature 
changes, l, is nearly independent of temperature. The sat- 
uration vapor pressure, however, is dependent on T (via 
equation (6.17)). This makes the integration of the above 
equation very difficult. One way to overcome this difficulty is 
to assume that esw/T? is nearly constant for small tempera- 
ture differences and then treat the differentials as differences. 
In this case, for T = 10°C, esw = 12.16 mb, (10°C) = 
2.4774 x 10° Jkg™*, ôq = —12 x 10° Jkg™’, and c, © Cpa = 
1005 Jke~* K7}, we find that AT = —5.3°C. Thus, the final 
temperature is 4.7°C. Note that this result can be improved 
if we do not assume that cp © Cpa. However, this will require 
the value of w, which is not known. We could of course 
assume that w, © Wsw(10°C) and consider cp % Cpa + WswCw- 
In this case, we find that the difference in the result is 
small (~ 0.066°C). In order to find the decrease in vapor 


Figure 7.9 

Diagram for example 7.3. 
When two initially 
saturated air masses mix 
isobarically, the final 
product is supersaturated 
(T, e). In order for the 
mixture to reach the 
equilibrium state (T”, e’), 
the excess vapor will 
condense. The amount of 
vapor that condenses is 
proportional to the 
difference e — e’. 
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pressure we combine equation (7.81) with dwsw ~ (€/p)desw 
to arrive at 


l 
ôq = cp dT + a, 
p 


Then incorporating the C-C equation gives 


De 26h 
5q = (== an =e) desw. 
ly Cow p 





Here again if we treat the differentials as differences and use 
the values mentioned above we find that Ae,,, = —4.326 mb. 


Two equal air masses, both at 1000mb, mix adiabatically. 
Their initial temperatures and mixing ratios are T} = 293.8 K, 
w, = 16.3gkg"*, Ta = 266.4K, and wz = 1.3gkg~*. First 
show that the mixture is supersaturated. Then find the final 
temperature, the final mixing ratio, the amount of liquid water 
content per unit mass of air, and the amount of liquid water per 
unit volume of air. 


As the two air masses mix isobarically, the initial tempera- 
ture of the mixture will be (equation (7.75)) T = 280.2K. 
Since qı = 16.0386gkg~' and q2 = 1.2983 gkg™' we find 
that, for the mixture, q = 8.6685gkg~* (equation (7.74)), 
w = q/(1 — q) = 8.7443gkg"*, and e = wp/(e+w) = 
13.86 mb (using the approximate equations (7.76) and (7.78) 
results in T = 280.7K and e = 14.1 mb). 

For T = 280.2 K equation (6.17) yields esw = 10.05 mb. This 
result indicates that initially the mixture is supersaturated 
(e > esw). Thus, immediately after mixing the excess vapor 
will condense so that the mixture reaches the equilibrium 
state (T’,e’). This step, graphically shown in Figure 7.9, 
will result in a decrease in vapor pressure and an increase 
in temperature (remember the mixing is adiabatic, so the 


(7, e) 


(T, e) e =eylT') 
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(7.4) 


heat gain due to condensation remains in the mixture). The 
increase in temperature is denoted by the line connecting 
points (T, e) and (T", e) and the decrease in water vapor pres- 
sure by the line connecting points (T’, e’) and (T’,e). Since 
the process from (T, e) to (T”, e’) is isobaric and adiabatic it 
is isenthalpic. As such, equation (7.21) applies: 


(Cpa + Weew)(L — T’) = (w — w)ly 


or 
(Coa + wcw) (T — T’) = (e — e) 


Slo 


or 





(cpa + eS = oan exp (1988 a) | = 
where, assuming that the mixture is a closed system, 
w, = w = 8.7443¢ke~t. Numerical solution of the above 
equation, with cw = 4218Jkg~'K~! and l, = 2.4774 x 
10° Jkg™*, gives T’ ~ 282.8K. This corresponds to e’ = 
€sw(I”) = 12mb and hence to w = 7.55 gkg™". 

The amount of liquid water produced per unit mass of air 
is given by the difference q- q ~ w—w' ~ ṣ(e— e)a 
1.16gkg~*. This amount must be equal to the decrease of 
water vapor per unit mass, —dm,, which is proportional to 
the decrease in water vapor pressure, —desw. From the ideal 
gas law we know that the mass of water vapor per unit vol- 
ume at temperature T is e,,/R,T and its variation with 
temperature is 


desw/dT = desw/RyT — (esw/ R,T? )AdT. 


Using the C-C equation we can verify that the second term 
on the right-hand side of the above equation is much smaller 
than the first term. As such, the amount of water vapor 
lost is approximately equal to desw/ R,T. It follows that the 
amount of liquid water per volume is equal to desw/ R,T, or 


equal to 
e— e 


x 144gm. 
RT aoe 





The following problem (Salby, 1996) is a classic example, demon- 
strating the effect of the downslope wind called Chinook wind 
(Chinook in American Indian means “snow-eater”). Moist air 
from the Pacific is forced to rise over the continental divide. On 
the west side, the surface lies at 8300 mb where the temperature 
is 293 K and the mixing ratio is 15 g kg™*. If the summit lies at 
600 mb and assuming that any condensation products fall out, 
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find the surface temperature and relative humidity at a place on 
the leeward side at 830 mb. 


For T = 293K equation (6.17) gives esw = 23.4mb. From 
equation (7.4) we then find that the saturation mixing ratio 
at state (T, p, w): (293K, 800mb, 15gkg~*) is 18.74gkg*. 
Since w < ws, the parcel is initially unsaturated. As the 
air begins to rise over the mountain it cools and its relative 
humidity increases. From equation (7.36) we find that, for 
r = 15/18.74 = 0.8, Ticu = 288.7 K. 

The ascent to the level where the air first becomes saturated 
is described by Poisson’s equation 


_. (Tre \"" 
PLCL = P T 


where k = ka(1 — 0.26 w) = 0.2849. It follows that prcur ~ 
760mb. At this point the air is saturated. As it continues 
to rise condensation takes place and all condensation prod- 
ucts fall out. As such, the ascent above puc is described 
by a pseudo-adiabatic process. Therefore, 6,, from 760 to 
600 mb is conserved. From equation (7.63) we find, for T = 
Trcr = 288.7 K, p = prou = 760 mb, and w = wew(Trot) = 
15gkg™*, that 0 ~ 359K. Evaluating the same equation 
at p’ = 600 mb gives 








0.285(1—0.28w’) 
359 = T’ (=5) 
p' 
3376 
x exp wa + 0.81w") ( — 2.54) (7.86) 
Tuc 


where T’, p’, and w’ is the state of the air at the summit. In 
the above equation we know p’ but T’ and w’ are unknowns. 
However, w’ can be expressed as a function of T’ alone by 
combining equations (7.5) and (6.17): 


(T) 611 5417 
w! = Wey(T’) & É ( De: E exp (1983 - ) 
Pp P T 








Thus, in effect, equation (7.83) has only one unknown, T”. 
Solving this equation numerically gives 7” ~ 279K. It fol- 
lows that at the summit w’ ~ 9.6gkg~'. From that point 
on the air descends adiabatically with the mixing ratio 
remaining constant (9.6 g kg~*). From Poisson’s equation we 
then find that at 830mb, T ~ 306K. At this temperature 
Wsw & 38.5gkg~* and at p = 830mb the relative humid- 
ity is r = 0.25. Thus the air on the leeward side is warmer 
and drier than that on the west side. Note that we assume 
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(7.5) 


that the parcel from 760 to 600 mb remains saturated with 
respect to liquid water. This is a safe assumption because 
Ticu = 288.7K and the vertical distance between 760 and 
600mb in a standard atmosphere is about 1.5km. So, we 
would expect that the temperature at the summit would be 
above 0°C. 


What is the maximum rate of precipitation when saturated air 
rises? 


Consider a saturated parcel of air rising. For simplicity, we 
assume that this parcel is a unit mass. As the parcel rises 
water vapor condenses. If we assume that all condensed 
products (rain or snow) fall out we can estimate the upper 
limit for the rate of precipitation. In order to achieve this we 
consider equation (7.54) that describes reversible saturated 
adiabatic processes, which we simplify by assuming that l, = 
constant and Cpa + WiCw © Cpa- 





cpad In T — Ryd In pat bd (=) -~0 


or 


dWsw Wsw 








Cpad In T — Rad In pa + ly T T dT =0 
or by multiplying all terms by T 
ly Wew 
CpadT — T'Rad In pg + lydwsw — dT =0 





or 





l 
(us a) dT ad In pa + lydwsy = 0 


or using the hydrostatic approximation 








ly SW 
(ea 7 Jar + gdz + dwew = 0. 


For typical conditions in the atmosphere l,wey/T is very 
small compared to cpa. We can thus neglect it to arrive at 


CpadT + gdz + lyd = 0 


or 
Ll dWsw dT 


dz =e Ty RY 





or 
ly dwyy aT g 





T 
Cpa dz dz Cpa 
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Now, recall that T4 = a and that — = T,. Thus, the 


dz 
above equation becomes 








lL, dWsw 
Cpa dz A 
or 
dwsw — CpalTa — Ts) 
dz = 
or 


dWsw dt = Cpa(La a I.) 
dt dz ly 





where a is the rate of ascent or speed of the parcel, U. It 


follows that 
dWwsw Cpa(La — £3) 
dt l tU; 





The above equation gives the mass of condensed water per 
unit mass of air and unit time. If we multiply it by the mass 
of air, m, we can obtain an expression for the mass of water 
condensed and precipitated per unit time 


pe Coal Fa) ym 
or 
P= ae yyy 
or 
TSR 
P= waa T) yan, 


ly 


where A is the area and Az the thickness of the rising 
volume. Thus, per unit area and unit time 


Coal = Is) 


P= 
ly 


UpAz. 

In the above equation, P is given (in MKS units) in 
kgm~’s~?. In practice, the amount of condensed water and 
precipitation per unit area is expressed by the depth of the 
collected water in millimeters (mm). Considering that the 
density is about 10?kgm~*, then 1kgm~? is equivalent to 
1mm of precipitated water. If, in addition, we consider that 
1h = 3600s, then we can express the above equation in 
mmh!, which is the unit for the rate of precipitation used 
in practice 


(Ta -T;) 


P = 360024 UpAz. 


v 
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Problems 


(7.1) 


(7.2) 


(7.3) 


(7.4) 


(7.5) 


(7.6) 


(7.7) 


(7.8) 


Show that when two air masses mix adiabatically and isobar- 
ically without condensation taking place, the final potential 
temperature and final vapor pressure are the weighted 
averages of the initial values (make necessary assumptions). 
Show that 
PLCE __ p 
esw(TroL) Cael dew) 


Then use equation (6.17) to derive a relationship between 
Trcr, Taew, P, and Prow- 

An air mass has a temperature of 20°C at 970mb pres- 
sure with a mixing ratio of 5gkg~*. After some time this 
mass has acquired a temperature of 10°C and a pressure 
of 750mb. Assume that no condensation or mixing between 
the air mass and the environment takes place and calcu- 
late the initial and final values of vapor pressure, relative 
humidity, and dew point temperature. (7.74mb, 5.98 mb, 
0.33, 0.49, 276.5 K, 272.8 K) 

Show that evaporation of dm grams of water into ma grams 
of dry air (under constant temperature) requires an amount 
of heat absorbed given by 


ÔQ = mal, dw. 


A Chinook wind blowing at 800mb has a temperature of 
38°C and a mixing ratio of 4gkg™*. Is this the same as the 
air at 1000 mb on the windward side of the mountains having 
a temperature of 294.5 K and a mixing ratio of 10gkg~* or 
is it the same as the air at 800 mb with 278K and 5 gkg™*? 
(It is the same as the air at 1000 mb) 

Show that if an isothermal layer of temperature To con- 
tained between isobars p, and p mixes vertically, the final 
temperature distribution is given by 


T(p) = [To/(1 — k)] (pi * — p3") p” /(P1 — p2). 

What will be the lowest possible temperature to which very 
dry air at T = 40°C can be cooled by evaporation? Consider 
p = 1000 mb and l, constant and equal to 2.4 x 10° Jkg™'. 
(~288K) 

When it is cold outside our breath creates mixing clouds. Is 
it possible to create a mixing cloud when the temperature of 
our breath is less than the air temperature? Assume that the 
temperature of breath is 30°C. (Hint: Plot esw as a function 
of T using equation (6.17) and make assumptions about this 
curve for high temperatures. Then use the theory discussed 
in section 7.3.1.) 
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(7.9) 


(7.10 


= 


(7.11) 


(7.12 


ai 


(7.13) 


(7.14 


Na 


(7.15) 


On a summer day in Milwaukee the maximum temperature 
is predicted to reach 35°C with dew point around 28°C. 
What will be the lowest temperature to which the air can 
be cooled by evaporation in this case? Does the answer make 
sense? Elaborate on the difference between this answer and 
the answer in problem 7.7. Consider l, constant and equal 
to 2.4 x 10° Jkg~*. (~303 K) 

Moist air moves over land where the temperature is 20°C 
and the relative humidity is 70%. Through contact with the 
ground the air cools. At what temperature will fog form? 
Consider l, constant and equal to 2.45 x 10° Jkg~*. (287.5 K) 
The main cabin in an airplane is maintained at a pressure of 
900 mb and at a temperature of 25°C. If it suddenly decom- 
presses adiabatically to a pressure of 500 mb, what should 
the relative humidity be to avoid cloud formation in the 
cabin? (6.29%) 

Outside air has a temperature of —15°C and a relative 
humidity of 0.6. If the air indoors is 25°C, what is the rela- 
tive humidity inside the room (assuming that the air inside 
is only heated, not humidified)? If the room has a volume 
of 100m? how much water vapor must be supplied so that 
the relative humidity rises to 50% and what will then be the 
total mixing ratio? If the temperature change due to evap- 
orating water is neglected, what amount of heat must be 
added for this to happen? Assume l, constant and equal to 
2.5 x 10°Jke~*. (3.6%, 1.08 kg, 8.97 gkg~*, 2.7 x 10° J) 
Which of the following quantities are conserved in (a) 
reversible adiabatic unsaturated transformation, (b) rever- 
sible adiabatic saturated transformation, (c) adiabatic 
unsaturated isobaric transformation, (d) adiabatic satu- 
rated isobaric transformation? 

W, Q, Wey Its €, Csws Wsw, T, PLOL, Licr; Tw, 9, Ovirt, Oe, Ow, S, h. 
A refrigerator of volume 2m? is initially filled with air at 
temperature 303K and with a relative humidity of 50%. 
At what temperature condensation will form on the walls? 
Assuming that the desired final temperature is 275 K, how 
much water vapor must condense to achieve the final tem- 
perature? What is the total amount of heat that will be 
given away to the surroundings during the whole process? 
Assume l, = 2.5 x 10°Jkg™*, cy = 4218 Jkg™ Kt, pa = 
1.293 kg m~’. (~ 19°C, 19.8 g, 117368 J) 

If the temperature inside a room is 25°C and the tem- 
perature outside is —10°C calculate the maximum rela- 
tive humidity that can be accommodated inside the room 
without windows fogging when (a) the interior part of the 
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(7.16 


(7.17 


(7.18 


(7.19 


(7.20 


pe 


a 


) 


) 


Sinise 


windows is not thermally isolated (i.e. heat is transferred 
freely) from the exterior part, and (b) it is thermally iso- 
lated. Assume l, = 2.5 x 10°Jkg~', cw = 4218 Jke~'K-}. 
(0.32, 1.0) 

If the wet-bulb temperature of a parcel at 1000 mb is 15.6°C 
and its mixing ratio is 6gkg~’, find the relative humidity 
of the parcel after it is lifted dry adiabatically to 900 mb. 
Assume l, = 2.47 x 10° Jkg™*. (36.14%) 

Raindrops at 10°C evaporate into air of temperature 20°C. 
If the saturated mixing ratio of air at 10°C is 8g ke‘, what 
is the mixing ratio of the environmental air? Make assump- 
tions and consider l, = 2.47 x 10° Jkg™*. (3.93 gkg~*) 

If water vapor comprises 1% of the volume of air (i.e. 
it accounts for 1% of the molecules) what is the virtual 
temperature correction? (Tuir, = 1.0038 T) 

If all the water vapor in a column of air extending from the 
surface to the troposphere were to condense and fall to 
the surface, it would occupy the lower dw of the height of 
the column. In this case, we can write that 


putty = f pydz 
0 


where pw and py are the densities of liquid water and water 
vapor, respectively. Show that the upper limit for d,, is given 
by 





EA Toeso 
Kai Pwly LD 


where I’ is the lapse rate, which is assumed to be constant, 
Ty is the surface temperature, and eso is the saturation vapor 
pressure with respect to water at To. Get an estimate of dw 
by assuming realistic values for Tọ and I’. What do you 
think of the estimated dw? Is it realistic? [dw is often called 
precipitable water.] 

Show that the mixing ratio of a parcel of air is given exactly 
by 





TrWs 


= gee 


(7.21) In Meteorology Book I, Aristotle says: “Both frost and dew 


are found when the sky is clear and there is no wind.” What 
do you think? Is he right or wrong? Explain why. 


CHAPTER EIGHT 





Vertical stability in the 
atmosphere 


Over large scales the atmosphere is very nearly hydrostatic, which 
means that the pressure gradient force balances the force of gravity. 
Because the net force in the vertical is zero, the atmosphere over 
those scales exhibits slow vertical motions of constant speed, i.e. 
there is no upward acceleration of air. However, on smaller scales 
hydrostatic equilibrium may be invalidated. In this case, acceler- 
ated motion leads to convection which embraces many phenomena 
observed in the atmosphere, from the structure of planetary bound- 
ary layers to the dynamics of hurricanes. In this chapter we will 
investigate the conditions that allow such accelerating motions. 
In particular, we will examine the fate of a parcel at equilibrium 
with its environment when it is subjected to a small perturbation. 
We will make the following assumptions: (1) the environment. is 
in hydrostatic equilibrium; (2) the parcel does not mix with its 
surroundings; (3) the parcel’s movement does not disturb the envi- 
ronment; (4) the process is adiabatic; and (5) at a given level the 
pressure of the environment and the pressure of the parcel are equal. 
In problem 8.1 you will be asked to elaborate on the validity of the 
above assumptions. 


8.1 The equation of motion for a parcel 


Because the environment is in hydrostatic equilibrium the following 
equation holds: 

dp _ 

ie Pg: 
For the parcel, the above equation is not applicable because as the 
parcel is displaced (upwards or downwards) it has some acceleration 
(d?z/dt?). Assuming that gravity and the pressure gradient force 
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are acting on the parcel, Newton’s second law dictates that per unit 
volume (primes identify the parcel) 


pee 2 dp 
Pae 99 a@ 





or , 
A 

dz 
Because of assumption (5) above, dp'/dz = dp/dz. Thus, the above 
equation can be written as 


as (9 
Z=-g-—a’(-= 
a 
£ a —a 
z= 9 
a 


s= (Z£). (8.1) 


The right-hand side gives the force per unit mass acting on the 
parcel due to the combination of gravity and pressure gradient and 
is called the buoyancy of the parcel. Using the definition of virtual 
temperature and assumption (5), we can write the ideal gas law for 
both environment and parcel as: 


Z=>-g- 


or 





or 





p= PRaTvirt 
and 


p= p RaT iir 


(8.2) 


The reason for adopting these expressions rather than p = pRT and 
p = pP'R'T is that this way we deal with Ty, and 7’... rather than 
R, R’,T, and T’. This makes the analysis more straightforward. 
Combining equations (8.1) and (8.2) yields 


TY, = Ty ir 
i= g (Aa). (8.3) 
virt 


What we are interested in here is to investigate the effect on the 
motion, described by equation (8.3), of small displacements (z < 1) 
from an original equilibrium level. If, for simplicity, we take this 
level to be the z = O level where the temperature is Tyist and 
express both Tir and T’, in terms of Taylor’s series we get 


irt 








dT int 1 OT int 2 
Tr at Ree ek ae N 
virt virt,0 de vA 5 T2 
and m ‘a 
fer = Tyirt,o virt z virt 2 











dz 2 dz? 
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If we now neglect terms of order higher than the first and define the 
environmental virtual temperature lapse rate as —dTyirt/dz = Dvint 
and the parcel’s virtual temperature lapse rate as —dT’,.,/dz = T! 
we arrive at 


irt 


a Ivins — Toirt)? 
e ee as OR 8.4 
i Tyirt,o = Doirt ( ) 
In the above equation we can manipulate the term 1/(Tyirto—Lvirt2) 
to arrive at 
1 1 1 
Tyirt,o 3 Iire? E Tyirt,o 1 — vinz 


Tvirt,0 


1 a) Dyirt2 
1+ because <1. 
Tertio Tyirt,o virt ,0 





Q 








Then equation (8.4) becomes 


as g Dirt Z 1 
= 1 Dir a I i 
í Tyirt,o ( z Tyirt,o ) ( i viet) 








or by eliminating terms involving z? 


g (Doire — T? 


T virt 
virt,0 


\Z 


ž == 





or 
EEE i 


virt 
virt,0 





a Dyirt) z=0. (8.5) 


8.2 Stability analysis and conditions 


The solution of the differential equation (8.5) depends on the 
constants. Three possibilities exist. 


(1) Tir = Toii > 0 
In this case equation (8.5) takes the form ž + A?z = 0 and has 
the solution 


z(t) = A sin At + B cos At 


where the oscillatory components are characterized by 





E y T Pi — Pan) > 0. 
This is the so-called Brunt-Väisälä frequency. Since we 
assumed that the initial level is z = 0, it follows that in this 
case B = 0 and as such z(t) = A sin At. This indicates that 
the parcel will oscillate in time about its original position with 
a period T = 27/A. This represents a stable case as the parcel 
does not leave the original level. 
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(2) Leire = Tvis < 0 
In this case equation (8.5) takes the form ž — A?z = 0 and has 
the solution 


z(t) = Ae™ + Be* 


where now 





\= y Fo Da- Din) > 0. 
Since at t = 0, z(0) = 0 it follows that A+ B = 0. This indicates 
that A = —B # 0 (the possibility of A = B = 0 leads to 
the trivial solution z(t) = 0 which we ignore). Since A # 0 
it then follows that as t > oo, a parcel’s displacement grows 
exponentially. Note that since for t > œo, dz/dt = Are, the 
parcel’s motion is an accelerating motion. This is the unstable 
situation where the parcel leaves the original level and never 


returns. 

(3) Lit — virs = 0 
In this case equation (8.5) becomes % = 0 and has the linear 
solution 


z(t) = At+B 


indicating that the displacement grows linearly with time. The 
parcel’s motion is now one of a constant speed (dz/dt = A). 
This is the neutral case where the parcel leaves the original 
level and never returns. Thus, the only difference between neu- 
tral and unstable situations is whether or not the motion is an 
accelerating one. 


The physical meaning of the above is that as long as the parcel’s 
lapse rate is greater than that of the environment and the parcel is 
forced upwards (downwards), it will become colder (warmer) than 
the environment and it will sink (rise) back to the original level. 
If the parcel’s lapse rate is less than that of the environment and 
the parcel is forced upwards (downwards) it will become warmer 
(colder) than the environment and it will continue to rise (sink) 
thus moving away from the original level. The difference between 
the two lapse rates is the net force acting on the parcel. If the net 
force is zero the motion is one of a constant speed. Otherwise the 
net force acts in the direction of the initial impulse if I” pe —Lvire < 0 
or in the opposite direction if I¥,.,— Ivins > 0. 
Let us now examine in detail what happens when a parcel is dis- 
placed in a layer, which is characterized by a virtual temperature 
lapse rate Iir. From equation (7.9) we have that for an unsatu- 
rated parcel T’ = (1 + 0.61w’)T’, where w’ is the initial mixing 
ratio, which because of assumption (2) remains constant. Thus, 

AT ine n dI" 

a (1+ 0.61w’) re 
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or 
PaA 


virt 


= (1 +0.61w") Ia 
or using equation (7.29) 


T! = (1 + 0.61w’)(1 — 0.87w') T; 


virt 


or 
T! a œ (1 — 0.26w") T). 


virt 


The term 0.26w’ in the above relation is rather small and therefore 
we may neglect it. Thus, for unsaturated parcels 
T! 


virt 


aly. (8.6) 
Note that for the environment, w is not constant but a function of 


z. As such, for the environment we have 


dTor dT d 
= = (1+ 0.6lw) + (IT — 
z 


dz dz 
or 
d 
Pine = (1+0.61w)P — 0.61. (8.7) 


The second term on the right-hand side of equation (8.7) is not 
negligible and therefore we cannot arrive at a relation similar to 
equation (8.6) for Ivin and I (the environmental lapse rate). This is 
why we have to use the virtual temperature in our stability analysis 
rather than the actual temperature. If Iyi were approximately 
equal to I’, then we could use T rather than Tirs in our stability 
analysis. Because this is not the case, it is more appropriate to use 
virtual temperatures. 

We can then state the condition for stability for an unsaturated 
parcel in a layer with J\;,4 as follows: 


if Iun >a layer is unstable, 
if Ive =a layer is neutral, (8.8) 
if Ivun <Ia layer is stable. 


For saturated parcels the mixing ratio decreases with height. 
Proceeding as above we have 





Tea = (14 0.61w')T" 
or 
aT". dT’ dw’ 
=(1 61w’ 617’ — 
dz (1 + 0.61w") —— + 0.6 T: 
or 
dw’ 


rr, =(1+0.61w')P, — 0.617’. 
(1 +0.61w) T, — 0.617’ 


virt 
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Figure 8.1 

Relative position of 
fundamental lines and 
regions of stability. 
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In this case the second term of the right-hand side is much smaller 
than the first term. Thus, we can approximate the above equation as 
I 


virt 


wo Las 


Accordingly, the condition for stability for a saturated parcel in a 
layer with yirs can be stated as 


if Ive > I, layer is unstable, 
if Lyin =I, layer is neutral, (8.9) 
if Daire < I, layer is stable. 


Since T4 = 9.8°Ckm™* and T, ~ 5°Ckm™t, conditions (8.8) and 
(8.9) can be combined in one as follows: 


if Darn > Ta layer is absolutely unstable, 
if I; < Isin < Ta layer is conditionally unstable, (8.10) 
if Dan <I layer is absolutely stable. 


The word “absolutely” indicates that the stability criterion holds 
for any type (saturated or unsaturated) of parcel. The term “condi- 
tionally unstable” means that the layer is stable for displacement of 
unsaturated parcels and unstable for saturated parcels (Figure 8.1). 
Now recall that for unsaturated environments equation (7.28) 
applies: 

1000 ) i 


Oyirt = Tyirt (= 
P 


Logarithmic differentiation of the above equation gives 


1 dO virt = 1 AT yixt ka dp 


vine dz Tyirt dz Pp dz 
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or 

1 dO virt 1 ka p 

= rT virt g 
virt dz Tyirt Pp RaTyirt 

or 

1 dO ir = I virt +4 1 g 

Ovirt dz T Tiiri Tyirt Cpa 
or 

1 dOyir 1 
t= (La — Dain) - (8.11) 


virt dz Tyirt 











Combining (8.8) with (8.11) yields an alternative way to express 
stability conditions for unsaturated parcels: 
if vine >0 layer is stable, 
dz 
Ovir 1 : 
if dO vi = 0 layer is neutral, (8.12) 
dz 
dôi 3 
if = <0 layer is unstable. 
Z 


For a saturated parcel similar conditions can be applied if we substi- 
tute Ovire by Oe (or in practice Oep) which is invariant along saturated 
adiabats (or pseudoadiabats) 








dé 
if 2 >0 layer is stable, 
dz 
0 
if Ws =0_ layer is neutral, (8.13) 
dz 
0 
if as <0 layer is unstable. 





dz 


The above set of criteria applies to the stability as a parcel rises 
in a motionless layer. If this is not true (for example, when entire 
layers are lifted or lowered), then the stability of the layer may be 
affected. Here is why it happens. 

Let us consider a stable layer in hydrostatic equilibrium whose bot- 
tom has a pressure pı and whose top has a pressure pə. Let us 
further assume that the difference in pressure, Ap, between top 
and bottom remains constant during the lifting or sinking. If the 
whole process is an unsaturated adiabatic process, then equation 
(8.11) applies. Since in this process Oyj; is conserved then the dif- 
ference in Ovit between top and bottom is conserved. It follows 


that 
1 dO ire = 1 dO yixt 


vine dz 7 Dirt dp ( p9) ~ Tyirt 





(La Tair) 
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Figure 8.2 
Graphical representation 
of convective instability. 
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or 
1 dO virt _ Ra 


Ovirt dp 7 P g 





(Ta — Toirt) = constant 


or 
Ralla — Dvir) = constant x g x p. 


During lifting pressure decreases, and according to the above equa- 
tion Iyi 4 Ia. Since the layer is initially stable this tendency will 
cause the layer’s lapse rate to approach the region of instability 
(recall Figure 8.1). During sinking the opposite is true. As such, 
lifting of a layer tends always to destabilize the layer and sinking 
of a layer tends always to stabilize the layer. 

If the lifting process is strong enough to cause the layer to saturate, 
then the picture is quite different. In this case the stability changes 
depend on the way saturation is reached, which in turn depends on 
the vertical distribution of moisture and the vertical structure of 
temperature among other factors. 

For example, let us consider an initially stable unsaturated layer 
at level 1 in which the temperature profile is BT and in which 
dé./dp < 0 (i.e. the bottom of the layer has more moisture than 
the top of the layer). As Figure 8.2 illustrates when the layer is 
lifted, the bottom of the layer (B) and the top of the layer (T) 


p(mb) 






Level 2 


0; b2 > 0; T 


Figure 8.3 
Graphical representation 
of convective stability. 
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will cool at the dry adiabatic rate (solid lines). However, the bot- 
tom will reach saturation sooner than the top. Thus, at level 2 the 
bottom will start cooling at the moist adiabatic rate (broken line) 
whereas the top will continue cooling at the dry adiabatic rate. 
Because of that, at position 3, independently of whether or not the 
top has reached saturation, the temperature profile will be rotated 
to the left in relation to the initial profile at level 1. This coun- 
terclockwise rotation causes an initially stable situation to become 
less stable (see Figure 8.1). Thus, when d6./dp < 0 the initially 
stable unsaturated layer is potentially or convectively unstable. 

If in the same layer d0./dz > 0 (i.e. the top has more moisture than 
the bottom), then the top will reach saturation first and it will begin 
to cool at the moist rate while the bottom will continue cooling at 
the dry rate. In this case the initial temperature profile BT turns 
clockwise (Figure 8.3), which means that the initially stable layer 
will tend to become even more stable. Thus, when d6./dz > 0 
the layer is potentially or convectively stable. If in the same layer 
dé. /dz = 0 (i.e. all levels of the layer have the same moisture), then 
both top and bottom will reach saturation along the same moist 
adiabat. In this case the initially stable layer is called potentiallly 
or convectively neutral. 


p(mb) 





Level 3 






Level 2 
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If the whole layer becomes saturated, the stability criteria for the 
layer translate to the following 


























dé. , s ; 

if q >0 the saturated layer will be stable with 

@ respect to saturated parcel processes, 
. dhe i p f 
E 0 the saturated layer will be neutral with 

E respect to saturated parcel processes, 
.- dh : = : 
iog <0 the saturated layer will be unstable with 

f respect to saturated parcel processes. (8.14) 


where d6,/dz refers to the initial profile of 6, (or in practice Oep) 
in the layer, which initially is not saturated. In the atmosphere, 
a potentially unstable layer tends to form cumulus-type clouds 
and often convective precipitation. On the other hand, a poten- 
tially stable layer produces stratiform clouds with light (if any) 
precipitation. 


8.3 Other factors affecting stability 


In our initial assumptions we accepted that the parcel does not mix 
with the environment, that the process is adiabatic, and that the 
environment is not disturbed (i.e. there are no compensating ver- 
tical motions in the environment as a parcel or a layer rises in it). 
When these assumptions are violated they affect the temperature 
lapse rate of the parcel and so the stability criteria derived previ- 
ously do not apply. A detailed treatment of these subjects is beyond 
the scope of this book. Interested readers can consult Iribarne and 
Godson (1973). 


Examples 


(8.1) The table below provides pressure, temperature, and dew point 
temperature measurements at various pressure levels. (1) Inves- 
tigate the stability conditions of parcels rising in each layer. (2) If 
each layer were lifted en masse until it became saturated would 
it become convectively stable, neutral, or unstable? (Assume that 
ly is independent of T and equal to 2.45 x 10°Jkg'.) 


First, we need to establish whether or not at a given level 
the parcels are saturated. The temperature is at all levels 
greater than the dew point temperature except at 850 mb. 
Accordingly, for the stability of the parcels, conditions 
(8.12) apply for all layers except for the layer 850-800 mb. 


EXAMPLES 











P T Taew 
(mb) (°C) (°C) 
1000 30 23 

950 27 21 

900 23 20 

850 20 20 

800 18 10 

750 15 5 

700 10 2 

650 5 0 
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For the layer 850-800 mb, conditions (8.13) apply. For the 
convective stability of each layer conditions (8.14) apply. 
Thus, in order to answer the questions we need to calcu- 
late Ovit and Oep at each level. In order to estimate Oyirt 
and ep (recall equations (7.28) and (7.63)) we need to 
find at each level w and Trcr. The mixing ratio can be 
estimated by combining equations (7.3), (7.7), (7.16), and 
(6.17). Tuc, can be estimated from equation (7.36). After 
plenty of calculations we arrive at the following table; here e, 
is in mb. 








Taew 


Tyirt 


vir 


Oep 


€s 


Trcr 





1000 
Layer 1 

950 
Layer 2 

900 
Layer 3 

850 
Layer 4 

800 
Layer 5 

750 
Layer 6 

700 
Layer 7 

650 


303 


300 


296 


293 


291 


288 


283 


278 


296 


294 


293 


293 


283 


278 


275 


273 


0.0182 


0.0169 


0.0166 


0.0176 


0.0099 


0.0149 


0.0063 


0.0058 


306.4 


303.0 


299.0 


296.1 


292.8 


290.6 


284.1 


279.0 


306.4 


307.5 


308.1 


310.2 


312.1 


315.5 


314.6 


315.6 


357.4 


355.1 


355.0 


360.2 


340.1 


361.7 


333.3 


332.9 


0.66 


0.70 


0.83 


1.00 


0.61 


0.51 


0.58 


0.70 


43.1 


36.0 


28.2 


23.4 


20.6 


17.0 


12.2 


8.6 


294.3 


292.6 


292.3 


293.0 


281.3 


276.0 


273.4 


272.0 








(1) According to the results in the above table and condi- 
tions (8.12) and (8.13), we find that for parcels rising 
in each of the layers, layer 1 is stable (d6,;./dz > 0), 
layer 2 is stable (dO,i4/dz > 0), layer 3 is stable 
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(8.2) 


(dôve /dz > 0), layer 4 is unstable (d0.,/dz < 0), 
layer 5 is stable (dO,i4/dz > 0), layer 6 is unstable 
(dOyir/dz < 0), and layer 7 is stable (d6,;../dz > 0). 


(2) Similarly, we find that layer 1 is convectively unstable 
(d0.p/dz < 0), layer 2 is convectively neutral (d6.,/ 
dz ~ 0), layer 3 is convectively stable (d0.,/dz > 0), 
layer 4 is convectively unstable, (d0.,/dz < 0), layer 
5 is convectively stable (d6,,/dz > 0), layer 6 is 
convectively unstable (d60.,/dz < 0), and layer 7 is 
approximately convectively neutral (d6.,/dz ~ 0). 


In an atmospheric layer the virtual temperature is constant and 
equal to 10 °C. If an initial fluctuation causes the parcel to rise 
adiabatically, calculate the energy per unit mass that has to be 
given to the parcel in order for it to rise | km above the bottom 
of the layer before it begins to sink. 


When a parcel rises in the atmosphere, a certain amount of 
work is performed by or against the buoyancy force depend- 
ing on whether the motion is along or against the direction 
of the buoyancy force. If the buoyancy force is directed 
downwards (negative buoyancy) a certain amount of work is 
done against buoyancy and if the buoyancy force is directed 
upwards (positive buoyancy) a certain amount of work is 
done by the buoyancy force. In the above problem since the 
layer is isothermal the layer is stable. Thus, the buoyancy on 
the parcel is negative and a parcel forced to rise will have to 
return to the initial level. But before it begins its return, it 
will reach a maximum height which depends on how strong 
the initial impulse was. Recall from Chapter 4 that 








OW = Fdz = madz = mdz. 


Thus, the work done, W, when a parcel is forced to rise from 
a level 7 to a level f is 


f 
wf mdz. (8.15) 


Using equation (8.3), the above equation becomes 


W= [om a |! Se (8.16) 


w=omf TOig f dz 
Toirt z) i 





or 
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or, assuming that the ascent is dry adiabatic, 


w=om l Ivo — sa: gm(f — i). 


Toirt 





Since Tuire(2) = Tuirt,o and assuming that the bottom of the 
layer is at z = 0 (i.e. i = 0m and f = 1000 m), we find that 
the work done per unit mass is equal to —169.7 J kg™t, i.e. 
work is done against the buoyancy force. Since the ascent is 
assumed to be adiabatic, 6Q = 0. Then from the first law 
we obtain 


du = 169.7 Jkg™'. 


This is the amount of energy per unit mass that must be 
given to the parcel for the above process to take place. 


Note that using the hydrostatic approximation, equation 
(8.16) can be expressed as 


ir = Tyi 
W-= A z (Ze vint :) dp 
Tyirt 
Pp = PRaT virt 


f d 
WaR fg Tad E 
i p 





or using 


or 


W = -Rm f (r art — Lvirt)d In p. (8.17) 


In a (T — Inp) diagram the above equation is proportional 
to the area enclosed between the vertical profile of the 
environment and the vertical profile of the parcel’s virtual 
temperature. Equation (8.16) or (8.17) gives the maximum 
work done by the buoyancy and when divided by the total 
mass defines (remember we are always dealing with adia- 
batic processes, i.e. ÂQ = 0) what is called the convective 
available potential energy (CAPE) and the convective inhi- 
bition energy (CINE). More on CAPE and CINE and their 
applications will be discussed in the next chapter. 


Problems 


(8.1) Elaborate on the five assumptions made in derivation of the 
stability conditions of a parcel that is subject to a small 
displacement from equilibrium. 
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(8.2) 


(8.5) 


If the virtual temperature profile in a layer is Tyii(z) = 
Tyirt,0a/(@ + z), where a > 0 and Tyit.o > ala, derive the 
vertical profile of the virtual potential temperature. What 
is the condition for stability for unsaturated parcels in this 
layer? (Unstable for z < ze, stable for z > z., and neutral 


for z = ze, where ze = «/(aTyirto/La) — a) 


If the profile of virtual potential temperature in a layer is 





what is the condition for stability in the layer for unsaturated 
parcels? (Stable when z > a, In(a;/az), unstable when z < 
a, In(a,/az2)) 

An atmospheric layer is isothermal. A dry parcel is sub- 
jected to an upward displacement and begins to oscillate 
about its original level. Plot the oscillation period as a 
function of the temperature of the layer. What do you 
observe? 

If we define the geopotential ¢ according to dọ = gdz with 
¢(0) = 0 show that under hydrostatic conditions 


Az= (= In 2) Titis 
g P2 


where Az is the thickness of a layer bounded by p, and 
p2 (pı > p2) and Tin is the mean virtual temperature 
of the layer. Using the above relationship show that if 
a layer is lifted en masse while the mean virtual tem- 
perature remains constant, then the upper level pressure 
change, dpə, is related to the lower level change, dp, via the 
equation 
dp. _ dps 


Pı P2 


Show that for a constant virtual temperature lapse rate 
atmosphere (Tyire(Z) = Toirt — Dyirt2) 


Lire Z 


9/Ralvirt 
a] 


P = Po h = 
Show that for an atmosphere where Tyirt(2) = Tyirt,o 


—g9z/RaTvirt,o 


P = Poe 


From the following data determine the stability and convec- 
tive stability of each layer. 
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(8.9) 


(8.10) 


(8.11) 








p(mb) Tyire(°C) — Taew(°C) 





1000 30.0 22.0 
950 25.0 21.0 
900 18.5 18.0 
850 16.5 15.0 
800 20.0 10.0 
750 10.0 5.0 
700 —5.0 —10.0 
650 —10.0 —15.0 
600 —20.0 —30.0 








In an unstable layer of air extending from the ground to the 
pressure level of 900 mb, the virtual temperature decreases 
with height at a rate of 25°Ckm"*. A parcel of air at the 
surface is given an initial velocity of 1ms~. If the virtual 
temperature at the surface is 7°C, and assuming that the 
parcel is and remains dry, find its position and velocity after 
1min. (81m, 2.1ms~*) 

Sea-breezes and lake-breezes are established when the air 
over the land becomes warmer than the air over the water. 
The warmer air over land rises thus destroying the pressure 
stratification and resulting in horizontal pressure gradients 
between land and water. The height of the breeze is the 
altitude where the horizontal pressure gradient vanishes. 
Assuming that the vertical virtual temperature profiles over 
land and water are isothermal and the same, show that this 
height is given by 





Ra ln(pw/p) 
1 


1 
g( Tiy T 


virt,1 


h= 








) 


where the subscripts | and w denote the conditions over land 
and water, respectively. 

Consider problems 8.6, 8.7, and 8.10. How would the results 
change if these problems were stated in terms of the tem- 
perature profiles rather than virtual temperature profiles? 


CHAPTER NINE 





Thermodynamic diagrams 


Chapters 1 to 6 laid down the fundamental physical and mathemat- 
ical concepts pertaining to thermodynamics. While we always kept 
the discussion close to our atmosphere, it was not until Chapter 7 
that, in a more applied mood, we presented how these concepts 
are applied to yield quantities useful to atmospheric processes. 
However, even if we understand all the mathematics, we still need 
an efficient way to present and visualize thermodynamic processes 
in the atmosphere. Thermodynamic diagrams can do that. Up to 
this point we made an effort to visualize the thermodynamic pro- 
cess using a (p,V) or a (p,T) diagram. However, such diagrams, 
while simple, may not be very convenient when it comes to their 
utilization. 

Since the purpose of a diagram is efficiently and clearly to display 
processes and estimate thermodynamic quantities, the following are 
very desirable in a thermodynamic diagram: (1) for every cyclic 
process the area should be proportional to work done or energy 
(area-equivalent transformations), (2) lines should be straight (easy 
to use), and (3) the angle between adiabats and isotherms should 
be as large as possible (easy to distinguish). The (p,V) diagram 
satisfies the first condition (pda = dw), but the angle between 
isotherms and adiabats is not very large (Figure 4.5(a)). Because 
of this, while it is used for illustration purposes, this diagram is not 
used in practice. 


9.1 Conditions for area-equivalent transformations 


When we are constructing a new diagram, in effect, we go from the 
x = a,y = —p domain to a new domain characterized by two new 
coordinates, say u and w. These two coordinate systems are shown 
in Figure 9.1. If we want the first condition to be satisfied then an 
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Figure 9.1 
Area-equivalent 
transformations. 
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v 
cv 


—_ -> 
area dA = | dx x dy | in the x,y domain is mapped onto an area 
— — 
dA’ = | du x dw | in the u, w domain, such that dA = ¢dA’ (x here 
is the cross product between two vectors and ¢ > 0 is a constant). 
Since any pair of points x and y are transformed to u and w we 
have that x = x(u,w) and y = y(u, w). We can thus write 

















des z du +2 dw 
> Oy > Oy 
Then 
Ox Oy Oxdy\|> => 
_ (Oz oy _ Ox oy 1 
dA (Za sect) du x dw}, (9.1) 
where we have used the facts that 
du x du = dw x dw= 0 
and 
du x dw = — dw x du. 
In equation (9.1) | du x dw | = dA’. Thus, 
dA = ¢dA' (9.2) 
where 
dx ae 
Ou Ow 
C= J — 
dy Oy 
Ou Ow 


is the Jacobian of the coordinate transformation. It follows that 
the condition for area-equivalence is that the Jacobian should be a 
constant (if J = 1, then this is an equal-area transformation). The 
above approach tells us whether or not the two new coordinates 
provide an area-equivalent transformation, but both u and w must 
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be specified. If only one is specified, then the following approach 
provides the other one. Since dA = ¢dA’ it follows that 


f ude = $ -pda =¢ f wdu 


foda + ¢wdu) = 0. 


or 


The above condition is true if pda + Çwdu = dz is an exact 
differential. Thus, z = f (a,u) and 


dz(a,u) = ($) da+ ($) du. 


It follows that the sufficient conditions for an area-equivalent 


transformation are 
_ (02 
E= Oa} , 


Oz 
a (a) 


From the above two equations we have 


and 


Op Oz 

ðu), ~~ ðaðu 
¢ Ow = O?z 

da), aðu’ 


ON oe 


then the areas will be equivalent (again here if ¢ = 1 the areas will 
be equal). From equation (9.3) if we specify u we can determine 
what w should be in order for the transformation to be an area- 
equivalent transformation. 





and 





Thus, if 


9.2 Examples of thermodynamic diagrams 
9.2.1 The tephigram 


If we suppose that u = T then 


(si), = (ar), 
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Using the ideal gas law, the above leads to 


a 
du), a 


From equation (9.3) we then have 


Ow R 
S 


¢ (=) da = Pia 


Oa a 


or 


or after integrating 
w= Rin a+ f(T) (9.4) 


where f(T) is the constant of integration which we can choose. By 
manipulating Poisson’s equation we can write 


1 RTN" 
6 \1000a 


1 
Ina =; (In ð In T) +1In T +1n R—1n 1000 


or 








or 
Rilna=c, ln 64 f'(Z). 


If we combine the above equation with equation (9.4) and choose 
f(T) = —f'(T) and ¢ = c,, we arrive at an area-equivalent diagram 
with coordinates 


w=ln0 (9.5) 
u=T. 


Obviously, in this diagram the dry adiabats are straight lines and 
perpendicular to the isotherms (Figure 9.2). The equation for the 
isobars can be derived from Poisson’s equation for p = constant: 


In 6 = In T + constant. 


In a coordinate system (In 0, T) the above equation describes loga- 
rithmic curves. As expected, the work done on a unit mass of air in 
a reversible cycle (—q) is proportional to the cycle: —q = — ¢ Tds = 
=, $ Td ln 6 = —cCpAtepn Where Atepn is the area of the cycle in the 
tephigram* (positive if described counterclockwise). If the air is dry, 
C= "Cpa 


* The word “tephigram” originated from the coordinates of this diagram, which 


are T and entropy. Entropy was originally denoted as ¢ (phi), but is now denoted 
by S. 


Figure 9.2 

Isobars, isotherms, and 
dry adiabats on a 
tephigram. Also shown is 
the ws = 10gkg~+ 
saturation mixing ratio 
line (solid bold line) and 
the pseudo-adiabat 

Oe = 40°C (broken line). 
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9.2.2 The emagram 


Let us assume again that u = T. By taking the logarithm of the 
equation of state we get 


Ina=—-Inp4+ln R+I1nT. 
Combining this equation with equation (9.4) yields 
w=-C€C'Rlnp+¢C [Rn R+Rln T+ fT). 


If we now choose f(T) = —RlIn R— Rin T and ¢ = R, we arrive 
at the area-equivalent diagram with coordinates 
w=-—lInp (9.6) 
u=T. 
In this diagram (which is called an emagram — Energy per unit 
MAss diaGRAM) the isobars are parallel straight lines and per- 


pendicular to the isotherms. The dry adiabats are again derived 
from Poisson’s equation 


ln 0 = ln T + kg In 1000 — kg In p 
or 


1 
—ln p= E ln T + constant. 
d 
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Figure 9.3 400 
Same as Figure 9.2 but for 
an emagram. 
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Temperature (°C) 


Since the coordinates are In p and T, the above equation is not a 
straight line. In the emagram the dry adiabats are slightly curved 
lines (Figure 9.3). Their angle with the isotherms is about 45°. 
Here, again, the work done on a unit mass of air in a reversible 
cycle (— ¢ pda) is proportional to the area of the cycle: 


= f pda = f [—d(pa) + adp] 
= f |- RdT + adp] 
= - $ RaT + $ adp 
= fod 
= f rr 


= -R $ Ta~ m p) 
= —RAen 





where Aem is the area in the emagram enclosed by the cycle, which 
is positive if described counterclockwise. If the air is dry, R = Rg. 
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9.2.3 The skew emagram (skew (T-ln p) diagram) 


Let us now assume that u = — ln p. Then 


haaa 
6 da ee ðln p), 
Ow 
sn 


¢ (==) da = —pda. 
ln p 





or 


or 
Oa 
Integrating for p = constant yields 
Cw = —pa + F (ln p) 


or 
¢w = —RT + F(l p). 


We may choose F'(In p) = —€1n p where € is a constant. Then 
Cw = —RT—Eln p. 

If we choose ¢ = R, we arrive at a new diagram with coordinates 
w=-T-€/Rinp 
u=-—lIn p. 


Since the sign of an area in a diagram involves only the direction 
in which a cycle is carried out, the new coordinate system can be 
written as 


£ 
w=T+—2In 9.7 
pe? (9.7) 
u=-—In p. 
This represents a correction (in fact a rotation) to the emagram 


which provides larger angles between any adiabats and isotherms 
(Figure 9.4). An isotherm in this diagram has the equation 


w = constant + ln p 


or 


Ç 
w = constant — >u. 
R 


In the (u,w) domain this equation represents a straight line with a 
slope equal to £/ R. Thus, we may choose a € for which the isotherms 
are perpendicular to the dry adiabats. Note that the equation for 
dry adiabats is not linear. From Poisson’s equation we have that 
In 0=In T + ka In 1000 — kg In p or In p= (1/ka)ln T+ constant. 
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Figure 9.4 
Same as figure 9.2 but for 
a skew (T —In p) diagram. 
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While here In p is a coordinate, In T is not and so the dry adiabats 
are not exactly straight lines. Here again it is easy to show that 
the work done in a cyclic process is proportional to the area of the 
cycle. The skew emagram and tephigram are the most commonly 
used diagrams. Other diagrams have been proposed (see problems 
9.1-9.3), but their practical use is limited. A detailed skew T—In p 
diagram is supplied at the end of the book (Figure A.1). 

For a complete description of a thermodynamic process in a dia- 
gram, in addition to isotherms, isobars, and adiabats, we also need 
to define the constant mixing ratio lines and the saturated adiabats 
(see also Chapter 7). From the definition of mixing ratio we have 


Sw T 
p = esw(T) + aes ( ) 
Wsw 

or 
€€sw(T) 


Wsw 


(9.8) 





p7% 


For a constant Wsw the above equation defines a family of curves, 
which in the (T,lIn p) domain are approximately (but not quite) 


straight lines: 
6.11 5417 
pe i exp (1983 — =) 





T 


sw 


or 
11 41 
wpe E pig 
Wsw T 





Figure 9.5 

Graphical procedure to 
obtain thermodynamic 
variables, from a skew 


(T-ln p) diagram, given 


an initial state (T’, p). 
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or B 
l xA+ — 
np ty 


where A = ln (6.lle/wsw) + 19.83 and B = —5417. The above 
equation for typical values of Wsw and T is nearly a straight line. 
Since w = Wew(Tuew, p), the dew point of unsaturated air defines 
the mixing ratio w at (T, p). It follows that the mixing ratio at 
(T, p) defines the dew point temperature. 

The saturated adiabats are determined numerically from equation 
(7.63). They can be labeled in terms of the estimated fep or in terms 
of the related (through equation (7.64)) Owp- 


9.3 Graphical representation of thermodynamic 
variables in a 7'—In p diagram 


Figure 9.5 shows the graphical procedure to obtain all the ther- 
modynamic variables derived analytically in Chapter 7, for a skew 
T-ln p diagram. We start with the state (T,p) where the air is 
assumed to be unsaturated and to have a mixing ratio w. The 
intersection between the w = constant line and the temperature 
axis gives the dew point temperature. By following the dry adiabat 
passing through point (T,p) (line a,) down to 1000 mb we obtain 
the potential temperature 0. By following it up until we intersect 
the w = constant line we arrive at the lifting condensation level 
(pict), where the parcel first becomes saturated and has a temper- 
ature Tucu. Through point (Tuc, Pict) passes a pseudo-adiabat s 














200 mb 
DRY ADIABAT a, 
PSEUDO-ADIABAT s 
Pict 
<— DRY ADIABAT 
w= CONSTANT “2 
p a 
Tict Tdew 
1000 mb 
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Figure 9.6 

Sounding at Davenport, 
IA, on 21 June 1997 at 
0000 UTC. As explained 
in the text, the 
information provided by 
this diagram would not 
support the development 
of deep convection that 
day. Nevertheless, deep 
convection did occur later 
in the day and very severe 
weather resulted in parts 
of the Midwest. See text 
for details. 
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defined by 6.,. If we follow this pseudo-adiabat down to the original 
level p, we obtain the pseudo-wet-bulb temperature, Twp, and con- 
tinuing down to the 1000 mb level we obtain the pseudo-wet-bulb 
potential temperature, wp (this also indicates the one-to-one rela- 
tionship between 6,,, and fep hence the labeling of pseudo-adiabats 
in terms of either Oep or Owp). If we follow the pseudo-adiabat s up 
to a pressure level where we can assume that all vapor has con- 
densed and has fallen out of the parcel (typically this level is taken 
to be the 200 mb level), and then follow the dry adiabat passing 
through that point (line a2) down to the original level we obtain 
the pseudo-equivalent temperature, Tep. Continuing down to the 
1000 mb level yields the pseudo-equivalent potential temperature, 
Oep. Note that Toi, Tw, Ow and ĝe are not directly estimated in dia- 
grams. As we discussed in Chapter 7, even though they are very 
close to each other, in actuality Tei Æ Tep, Tw A Twp, Ow Æ Owp, and 
Oei £ Oep (recall equation (7.66)). 


9.3.1 Using diagrams in forecasting 
Figure 9.6 shows in a skew T-ln p diagram a sounding taken at 


Davenport, IA, on 21 June 1997 at 0000 UTC. The solid line is the 
vertical profile of environmental temperature and the heavy broken 
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line to its left is the vertical profile of dew point temperature. The 
parcels at the surface are unsaturated. The lifting condensation 
level of these parcels is about 860mb where Tcr = 20°C. Above 
the LCL the light broken line is the pseudo-adiabat passing through 
the point (Ttor, pict) and below the LCL the light broken line is the 
dry adiabat passing through point (T, p) at the surface. The pseudo- 
adiabat intersects the temperature profile at point A. Up to this 
point the parcel is colder than the environment. As such the layer 
from surface to point A is stable. This means that the buoyancy 
in this layer is negative and that work must be done against the 
buoyancy if parcels originating at the surface are to reach that 
level. The work needed is proportional to the darker shaded area 
(negative area). This work or the magnitude of this negative area 
is often called the convective inhibition energy (CINE) given by 
equation (8.17) for i = Psurface and f = pa (here and in practice we 
make no distinction between T and Tuire). 


PA 
CINE = Rf (T’ — T)dln p. 


Psurface 


In our example, Psurface % 980 mb and py, % 740 mb. Assuming that 
the average T’ — T in this layer is about —0.8°C and that R ~ Ra, 
we find that CINE ~ 64 Jkg™". If we now recall equation (8.15) we 
have that 





"A du 1 1 
CINE = i oe: vdt = 5 UA a Maurice 
where v is the vertical velocity. Assuming that va œ% 0 it follows 
that if a parcel that originated at the surface were to reach point 
A, the minimum initial impulse has to be 


Umin = V2CINE x 11.4ms7?. 
Therefore, unless an initial impulse of 11.4ms~! is applied to the 
parcels they will not make it to point A. This point is called the 
level of free convection (LFC). Above this level the parcel is warmer 
than the environment all the way up to 150 mb. The buoyancy from 
LFC to 150 mb is positive and now work is done by buoyancy. This 
work is proportional to the lighter shaded area (positive area). In 
this layer the parcels that made it to LFC are now free to rise 
and accelerate. Note that owing to the temperature inversion in 
the stratosphere, above 150 mb the temperature normally increases 
with height. As such, there is a point (B in our example) where 
the temperatures of the parcel and the environment become equal 
again. Above this level, which is referred to as the level of neutral 
buoyancy (LNB), work is done against buoyancy and the parcels 
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begin to decelerate. The magnitude of the positive area is called 
the convective available potential energy (CAPE) and is given by 

LNB 

CAPE = -r | (T’ — T)dIn p. 

LFC 
Assuming that between LFC and LNB the average difference 
between the temperature of the parcels and the environment is 
about 7°C and that R ~ Rg we can estimate that CAPE ~ 
3200 Jkg~*. As we did above, we can show that in this case the 
maximum velocity a parcel will attain is 


Umax = V2CAPE x 80ms"1. 


Both CINE and CAPE are very useful as they provide information 
on whether or not convection will occur (CINE) and on how severe 
a storm might become (CAPE). 

In our example CAPE is rather high. Thus, if the parcels make it 
to LFC, they will then accelerate greatly and deep convection will 
ensue. The problem in predicting deep convection in this particular 
example, however, is that the value of CINE is considerable. Such 
values may inhibit development of convection. In practice, situa- 
tions of high CAPE values should be monitored to see whether the 
negative area is removed by warming of the air in the lower levels, 
or whether upward acceleration is aided by a low level jet or fron- 
togenesis. In our case all of the above happened thereby allowing 
deep convection and severe weather to develop later in the day in 
the Midwest. 

Together with CAPE and CINE we may calculate the precipitable 
water from soundings. If we recall the definition of precipitable 
water in problem 7.19 and the definition of specific humidity in 
equation (7.2) we have 


1 co 
dy = — f qpdz 
Pw 0 


or using the hydrostatic approximation 


1 Psurface 
dy S= qdp. 
GPw Jo 

The above equation indicates that precipitable water can be esti- 
mated by integrating the specific humidity profile. The specific 
humidity profile as a function of pressure can be obtained from 
the vertical profiles of T and Ty using the equations presented in 
Chapter 7. 


Example 


(9.1) An air parcel at 1000mb has a temperature of 20°C and a 
mixing ratio of 10 gkg"'. The parcel is lifted to 700 mb by passing 
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over a mountain. During the ascent 80% of the water vapor 
that condenses falls out of the parcel. Find the temperature and 
potential temperature of the air when it returns to 1000 mb on 
the other side of the mountain. 


The parcel is located at the point with coordinates T = 
20°C and p = 1000mb (w = 10gkg™'). At this point, 
the saturation mixing ratio is ~ 16gkg~'. Therefore, the 
parcel is unsaturated and as it begins to rise it will fol- 
low the dry adiabat at that point (which because we are 
at 1000mb is 20°C). The point at which this dry adia- 
bat intersects the w = 10gkg™* line is at p = 910 mb and 
T = Tict ~ 12.4°C. This is the lifting condensation level 
and from that point on the parcel will follow the correspond- 
ing pseudo-adiabat (6., = 49.5 °C, which in the skew T—ln p 
diagram is labeled in terms of the related wp = 16.5°C). 
At 700mb the saturation mixing ratio is ~ 6.66gke’. 
Therefore, 10 — 6.66 = 3.34g of vapor per kg of dry air 
have condensed during the saturated ascent from 910 mb to 
700mb. Since 80% of this amount is removed by precipi- 
tation, it follows that 0.668 gkg~' of liquid water remains 
in the parcel as it begins its descent. This means that ini- 
tially the descent will be dictated by the same saturated 
adiabat until all the water evaporates, i.e. until the mixing 
ratio increases from 6.66 to 6.66 +0.668 = 7.328 gkg~*. This 
will happen at the point where the pseudo-adiabat inter- 
sects the w = 7.328 line. This is the point with coordinates 
p= 750mb and T ~ 4°C. From this point on, the descent 
is dictated by the corresponding dry adiabat (0 ~ 28°C). 
Thus, on the other side of the mountain at p = 1000mb, 
T = 0 = 28°C. 


Problems 


(9.1) 


(9.2) 


(9.3) 


(9.4) 


Show that the Refsdal diagram with coordinates u = ln T, 
w = —T ln pis an area-equivalent diagram and that neither 
its isobars nor its dry adiabats are straight lines. 

Show that the Stüve (or pseudo-adiabatic) diagram with 
coordinates u = T,w = —p*4 is not an area-equivalent dia- 
gram. Show that isobars, isotherms, and dry adiabats are all 
straight lines. 

Using the Jacobian method show that the diagram 
with coordinates u = T,w = c,In@ is an equal-area 
diagram. 

Given the sounding shown in the table, 
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(9.5) 


(9.6) 


(9.7) 


950 22.5 
900 18.0 
850 15.0 
800 16.0 
750 12.0 
700 7.0 
650 4.0 
600 —1.5 
500 —10.0 


400 —20.0 


determine using a diagram (1) the mixing ratio, (2) 
the relative humidity, (3) the potential temperature, (4) 
the pseudo-wet-bulb temperature, (5) the pseudo-wet-bulb 
potential temperature. Having determined the above, com- 
pute the wet-bulb temperature, the equivalent tempera- 
ture, the pseudo-equivalent potential temperature, and the 
pseudo-equivalent temperature. Assume that ly = 2.45 x 
10° Jkg™*, cw = 4218Jkg™*, and that the dew point at 
950 mb is 15.7°C. (12.4gkg~*, 0.68, 27.5°C, 18°C, 20°C, 
19°C, 51°C, 66°C, 60°C) 

Approximate CAPE if in the previous problem we extend 
the sounding as follows: 


p(mb) = T(°C) 


300 —35.0 
200 —55.0 


150 —60.0 


Would you expect deep convection to occur? 

Consider the following scenario (Figure 9.7). On a given day, 
on a skew T-lIn p diagram the temperature profile from the 
surface (1000mb) to 200mb has a constant lapse rate of 
7°Ckm™*. Above 200 mb the temperature is constant. If the 
parcels at the surface are saturated, what will the maximum 
vertical velocity be on that day assuming that the saturated 
adiabatic ascent can be approximated by the constant lapse 
rate of 6.5 °C km~'? The temperature at the surface is 30 °C. 
Assume R = Ra (~ 54ms~*). 

A parcel of air with an initial temperature of 15°C and dew 
point 0°C is lifted adiabatically from 1000 mb. Determine 
the lifting condensation level and temperature at that level. 
If the parcel is lifted a further 250mb what will its final 


PROBLEMS 173 


Figure 9.7 p (mb) 
Figure for problem 9.6 
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temperature be and how much liquid water is produced? 
(~ 800mb, ~% —3°C,~ —22 °C, ~ 2.7 gkg7™') 

(9.8) In the previous problem if all the condensed water falls as 
precipitation, can you roughly estimate the height of water 
collected on a rain collector of area 1m?? Here you have to 
make your own assumptions. How does this result compare 
with the result obtained using the formula in problem 7.19? 
Can you explain any differences? 

(9.9) An air parcel at 900mb has a temperature of 10°C and a 
mixing ratio of 5gkg~'. The parcel is then lifted to 700 mb 
by passing a mountain. If in this process 70% of the con- 
densed water is removed by precipitation, what will its 
temperature, potential temperature, and pseudo-wet-bulb 
potential temperature be when it returns to its initial level 
of 900mb on the other side of the mountain? (~ 14°C, 
~ 23.5°C,~ 11°C) 

(9.10) If a temperature profile in °C could be described as T(z) = 
28 — 8z for z > 1km and T = 20°C for z < 1km, and 
the mixing ratio at the surface is 16 gkg~*, estimate using a 
diagram the level of free convection. (~ 2.4km) 


CHAPTER TEN 





Beyond this book 


The previous chapters have presented the basics in atmospheric 
thermodynamics. As we know, in atmospheric sciences the ulti- 
mate goal is to predict as accurately as possible the changes in 
weather and climate. Thermodynamic processes are crucial in pre- 
dicting changes in weather patterns. For example, during cloud and 
precipitation formation vast amounts of heat are exchanged with 
the environment that affect the atmosphere at many different spa- 
tial scales. In this final chapter we will present the basic concepts 
behind predicting weather changes. This chapter is not meant to 
treat the issue thoroughly, but only to offer a glimpse of what comes 
next. 


10.1 Basic predictive equations in the atmosphere 


In the Newtonian framework the state of the system is described 
exactly by the position and velocity of all its constituents. In the 
thermodynamical framework the state is defined by the tempera- 
ture, pressure, and density of all its constituents. In a dynamical 
system such as the climate system both frameworks apply. Accord- 
ingly, a starting point in describing such a system will be to seek 
a set of equations that combine both the mechanical motion and 
thermodynamical evolution of the system. 

The fundamental equations that govern the motion and evolution 
of the atmosphere (and for that matter of the oceans and sea ice) 
are derived from the three basic conservation laws: the conserva- 
tion of momentum, the conservation of mass, and the conservation 
of energy. For the atmosphere the equation of state relates temper- 
ature, density, and pressure. In summary, from the conservation of 
momentum we derive the following set of predictive equations of 
motion (Washington and Parkinson (1986)) 
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dt Ot acoso OA a ð¢ 





dX 
u = a COs PTE 


dz 
dt 
where u,v, and w are the horizontal and vertical components of the 
motion, ¢ is the latitude, A is the longitude, a is the radius of 
the earth, f is the Coriolis force, and F is the friction force. In 
the above equations the forces that drive the motion are the local 
pressure gradients, gravity, Coriolis, and friction. 
Equation (10.3) can be approximated by the hydrostatic equation 
(by assuming that dw/dt = 0 and F, = 0) 
10 
(ae (10.4) 
poz 
which relates density to pressure. 
From the law of conservation of mass we can derive the continuity 
equation 
Op 1 o o o 
= u) + = (pv cos — — (pw 10.5 
T T 9, ) ag? o) gz P ) (10.5) 
which provides a predictive equation for density. The last two equa- 
tions in this basic formulation come from thermodynamics. The first 
law of thermodynamics (which expresses the law of conservation of 


energy) 





dT 1l1dp dQ 
Cy = i 

dt pdt dt 
where dQ/dt is the net heat gain, provides a predictive equation 
for temperature. The above equations (10.1, 10.2, 10.3 (or 10.4), 


(10.6) 
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10.5, and 10.6) make a system of five equations with six unknowns 
(u,v,w, p,p, T). The equation of state 


p = pRT (10.7) 


provides the additional equation that connects pressure, density, 
and temperature, thus resulting in a system of six equations (called 
the primitive equations system), with six unknowns (assuming of 
course that dQ/dt, Fy, and F are constants and known). Note, 
however, that this system is not a closed one because dQ/dt, Fy, and 
F, must be determined from the other variables. These terms are 
very important for climate simulations but for short-term weather 
prediction they are often ignored. 


10.2 Moisture 


The above system of predictive equations was derived without the 
inclusion of moisture. Even though a model can be derived without 
moisture in it, including moisture can dramatically improve mod- 
eling and prediction. In a manner analogous to the continuity of 
mass, the changes of moisture must be balanced by the moisture’s 
sources and sinks. An equation of continuity for water vapor mixing 
ratio can be written as 


dw å 
Y NF aL (10.8) 
p 


dt 


where M is the time rate of change of water vapor per unit volume 
due to condensation or freezing, E is the time rate of change of 
water vapor per unit mass due to evaporation from the surface and 
to horizontal and vertical diffusion of moisture occurring on scales 
below the resolution of the model, and w is the mixing ratio. Often 
the above equation is combined with equation (10.5) to obtain 


O( pw) 1 à) a ə 7 
ðt acoso By (oe) + pg Pwr cos ¢) +a, (eww) = M+pE. 
(10.9) 


When water vapor changes to water (or ice) and vice versa, heat 
is added to or removed from the atmospheric system. In this case 
dQ = d(lw,/T). 

The above general and simple equations constitute only the start- 
ing point for studying and investigating the fascinating dynamics 
of atmospheric motion. The equations have to be modified and 
adjusted in order to include the effect of many other factors (for 
example, radiative processes) and to be applied over a specific range 
of scales (mesoscale models, for example). Such studies are crucial 
in the effort to understand the climate and ecology of our planet. 
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Table A.1. Units 




















Unit 
Physical quantity MKS system egs system 
Length m cm 
Time s s 
Velocity ms * cms * 
Acceleration ms” cms? 
Mass kg g 
Density kg m~?’ gem? 
Force kg ms *(newton, N) gcms *(dyne, dyn) 
Pressure kgm~'s~?(pascal, Pa) gcm~*s~?(microbar, ubar) 
Energy kg m?s~*(joule, J) gcm? s~? (erg) 
Specific energy m’s *(Jkg~*) cm’ s~?(erg g~*) 
Relations: 
1N = 10° dyn 
1 Pa = 10 ubar 


1 bar = 10° ubar = 10° Pa 


1 atmosphere = 1013 mb = 1013 x 10? Pa = 1013 hPa 


1J = 10’ erg 
lcal = 4.185 J 
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Table A.2. Selected physical constants 








General constants 


Avogadro’s number, Na 6.022 x 107° mol7* 
Boltzmann’s constant, k 1.38 x 107 JK! 
Gravitational acceleration at sea level, g 9.807 m s7? 
Universal gas constant, R* 8.314 J K~* mol* 
Constants for dry air 
Molecular weight, Ma 28.97 g mol * 
Density, pa 1.293kgm~* (at STP)* 
Specific gas constant, Ra 287 Jkg7* KT 
Isobaric specific heat (at 273 K), Cpa 1005 Jkg~* K7* 
Isochoric specific heat (at 273 K), cva 718Jkg-*kK~* 
Ratio of specific heats, ya 1.4 
Constants for water 
Molecular weight, My 18.015 g mol* 
Density (water) 1000 kg m~°* (at STP) 
Density (ice) 917kgm7~° (at STP) 
Specific gas constant, Ry 461.51 Jkg~* K7* 
e = Ma/M, = Ra/Ry 0.622 
Isobaric specific heat (at 273 K) 1850 Jkg~* K7* 
(wapor, cp.) 
Isochoric specific heat (at 273 K) 1390 Jkg~* K7* 
(vapor, cvv) 
Specific heat capacity (at 273 K) 4218 Jkg~* K“ 
(water, Cpw © Cvw © Cw) 
Specific heat capacity (at 273 K) 2106 Jkg~* K7* 
(ice, Cpi = Cvi = Gi) 
Latent heats see Table A3 








STP: Standard temperature and pressure (273 K, 1013 mb) 
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Table A.3. Specific latent heats 








T (°C) Is (10°Skg-") i (10°Jkg™™) l (10° Jkg~*) 





—100 2.8240 
—90 2.8280 
—80 2.8320 
—70 2.8340 
—60 2.8370 
—50 2.8383 0.2035 2.6348 
—40 2.8387 0.2357 2.6030 
—30 2.8387 0.2638 2.5749 
—20 2.8383 0.2889 2.5494 
—10 2.8366 0.3119 2.5247 
0 2.8345 0.3337 2.5008 
5 2.4891 
10 2.4774 
15 2.4656 
20 2.4535 
25 2.4418 
30 2.4300 
35 2.4183 
40 2.4062 
45 2.3945 


50 2.3893 
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Figure A.I 

The skew T—In p 
diagram. Saturated 
adiabats are labeled in 
terms of Owp. 
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cloud formation 
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Dalton’s law, 20-21, 99 
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dew point temperature, 
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variation (of dew point) with 
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disorder, 68-69 
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processes, 38-40 
lapse rate, 41—42 
dry air 


composition, 22 
molecular weight, 22 


energy 
conservation of, 5 
definition of, 26 
internal, 5 
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potential, 5 
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enthalpy 
definition of, 33 
of fusion, 84 
of sublimation, 84 
of vaporization, 84 
of temperature dependence 
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entropy 
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and disorder, 68-69 

and irreversible processes, 
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and potential temperature, 
70 
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definition of, 62 
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for liquids and solids, 62 

of an isolated system, 64 

of a mixture, 72-74 

equal distribution of energy 
theorem, 9 


equations of motion in the 
atmosphere, 176 

equivalence between heat and 
work, 31 

equivalent (isobaric) potential 
temperature, 124 

equivalent potential 
temperature, 121 

equivalent (isobaric) 
temperature, 111 


equilibrium 
between vapor and liquid, 
80 
metastable, 3 
stable, 3 


unstable, 3 
vapor pressure 
for ice, 81 
for water, 81 
temperature dependence of 
phase transformation, 
84 
equilibrium state, 3 
exact differential 
definition, 7 
conditions for, 8 


131-133 
66-67 


freezing in a cloud, 
fundamental relations, 
fusion, 81 

latent heat of, 84 


gas constant 
for dry air, 22 
for moist air, 102 
specific, 18 
universal, 19 
Gibbs formula, 81 
Gibbs function, 66 


Guy-Lussac 
1st law, 13 
2>4 law, 14 


heat capacities (see thermal 
capacities) 
Helmholtz function, 66 
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hydrostatic approximation, 
42, 126, 143 
humidity 
absolute, 101 
relative, 101 
specific, 101 


ideal gas 
definition, 10 
law, 12, 17-19 

Joule’s law, 34, 96 


Kelvin’s postulate, 58, 59 
kinetic theory of heat, 9 


latent heat 
definition of, 84 
of fusion, 84 
of sublimation, 84 
of vaporization, 84 

lifting condensation level (LCL) 
definition of, 115 
estimation, 116-117 


Maxwell’s relations, 76 
mixing and irreversibility, 74 
mixing ratio, 101 
saturation, 101 
water, 109 
mixture of gases, 
moist adiabatic 
lapse rate, 113 
processes, 109, 112, 113 
moist air 
definition of, 99 
mean molecular weight, 102 
moisture in the equations of 
motion, 177 
molar heat capacity, 33 
mole definition of, 16 
molecular weight of dry air, 


20, 21 


22 
Poisson’s equations, 38-40 
potential temperature, 43—44 


of unsaturated moist air, 
112 
precipitation rate, 138 
precipitable water, 142, 170 
process 
adiabatic expansion of 
unsaturated moist air, 
112 
adiabatic isobaric, 109 
adiabatic isobaric mixing, 
128 
isenthalpic, 109 
isentropic, 119 
isobaric cooling, 104 
pseudo-adiabatic, 122-123 


saturated adiabatic, 
119-122 
saturated ascent, 118 
vertical mixing, 130 
pseudo-equivalent temperature, 
123 
equivalent potential 
temperature, 124 
wet-bulb temperature, 125 
wet-bulb potential 
temperature, 125 


raindrops, growth by 
condensation, 93 


relative humidity, 101 


saturated adiabatic lapse rate, 


126 
saturation temperature (at 
LCL), 115 


specific heat capacity 
at constant pressure, 33 
at constant volume, 33 
specific humidity, 101 
specific gas constant, 18 
stability criteria, 
for saturated parcels, 148, 
149 
for unsaturated parcels, 
147, 149 
of parcels when the layer is 
lifted, 149-152 
stability, convective, 151 
state 
equation of, 
functions, 2 
mechanical, 1 
thermodynamical, 2 
variables, 2 
sublimation, 81 
atent heat of, 84 
supercooled water, 
supersaturation, 
system 
closed, 1 
homogeneous, 
heterogeneous, 
isolated, 1 
open, 1,2 
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table of constants, 185 
table of units, 181 
temperature 
dew point, 105-109 
equivalent (isobaric), 111 
equivalent (isobaric) 
potential, 124 
equivalent potential, 121 
Kelvin scale, 15 


potential, 43—44, 112 
of unsaturated moist air, 
112 
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pseudo-equivalent, 123-124 
pseudo-equivalent potential, 
124 
pseudo-wet-bulb, 125 
pseudo-wet-bulb potential, 
125 
relationship between dew and 
wet-bulb, 112 
saturation (at LCL), 115 
thermodynamic definition of, 
67 
virtual, 102-103 
virtual potential, 113 
wet-bulb, 109-111 
wet-bulb potential, 125 
temperature lapse rate 
dry adiabatic, 41-43 
moist adiabatic, 113 
saturated adiabatic, 
126-127 
thermal (heat) capacity 
at constant pressure, 33 
for dry air, 36 
for moist air, 103-104 
at constant volume, 33 
of diatomic gases, 36 
of monatomic gases, 35 
thermodynamics 
definition of, 1 


statistical nature of, 67—70 
1st law of, 32, 37 
2°¢ law of, 62, 64 


combining the two laws, 66 
thermodynamic diagrams 

and forecasting, 168-170 

conditions for area-equivalent 
transformation, 
159-161 

emagram, 163 

graphical representation in, 
167 

skew emagram (T-ln p 
diagram), 165 


tephigram, 161, 162 
transformations 
adiabatic, 4, 38, 65 
cyclic, 4, 38 
irreversible, 3,8 
isentropic, 65, 71 
isobaric, 4, 37, 65 
isochoric, 4, 37, 65 
isothermal, 4, 37, 65 
polytropic, 40—41 


reversible, 3, 8 
triple point, 82 


units, 181 
universal gas constant R*, 
19 


vaporization, 81 
latent heat of, 84 
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variables 
dependent, 2 
extensive, 2 
independent, 2 
intensive, 2 
virtual potential temperature, 
113 


virtual temperature, 102-103 
water 
density of, 82 
mixing ratio, 109 
precipitable, 142, 170 
supercooled, 81, 94-95 
transformations of, 80-82 
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wet-bulb temperature, 
109-111 


wet-bulb potential temperature, 
125 


work 
and kinetic energy, 29 
definition of, 27 


